Data pre-processed and analyzed from chat logs and web logs by Marcus Collins.

Analysis based on R-Scripts from Jacob LaRiviere.

Experiment designs from Marcus Collins and Jacob LaRiveier.

R-Markdown and plots from Mike Wise.

Initialization of libraries

library(tidyverse,quietly=T,warn.conflicts=F)
library(lubridate,quietly=T,warn.conflicts=F)
library(scales,quietly=T,warn.conflicts=F)
library(zoo,quietly=T,warn.conflicts=F)
library(lmtest,quietly=T,warn.conflicts=F)
library(sandwich,quietly=T,warn.conflicts=F)
library(gridExtra,quietly=T,warn.conflicts = F)
library(knitr,quietly=T,warn.conflicts = F)

Set random seeds, record start time and version.

set.seed(1234)
version <-0.6
versionstring <- sprintf("Version %.1f",version)

starttime <- Sys.time()
startfmttime <- sprintf(format(starttime, "%d %b %Y - %H:%M:%S"))

print(sprintf("%s created on %s",versionstring,startfmttime))
## [1] "Version 0.6 created on 27 Mar 2017 - 20:44:09"

Various constants

tztz <- "UTC"  # Apparently all of our time zones are UTC...
s2date <- function(strdate){
  return(as.POSIXct(strdate,tz=tztz))
}
firstday <- s2date("2015-01-01")  # we will count days from the first day in 2015

smcdates <- c("2016-08-17/red/0-30%","2016-09-1/red/30-50%","2016-09-07/red/50-100%")
xbxdates <- c("2016-10-11/purple/0-10%","2016-10-18/purple/10-30%","2016-11-01/purple/30-50%","2016-12-15/purple/50-90%",
              "2016-11-17/blue/content change")
totdates <- c(smcdates,xbxdates)
smcback <- "lightsteelblue1"
xbxback <- "darkseagreen2"
xabback <- "darkseagreen3"
totback <- "wheat"

fpath <- "../TorontoData"
sdate <- s2date("2016-06-01")
mdate <- s2date("2017-01-01")
edate <- s2date("2017-03-07")
verbose <- 2
justdoone <- F # for testing

Misc utility functions

crackdate <- function(datestr){
  sar <- unlist(strsplit(datestr,"/"))
  sdate <- sar[[1]]
  date <- s2date(sdate)
  val <- 0
  sval <- sar[[3]]
  levpart <- sar[[3]]
  pctpresent <- F
  if (grepl("%",levpart)){
    levpart <- gsub("%","",levpart)
    val <- as.numeric(unlist(strsplit(levpart,"-"))[[2]])
    sval <- gsub("-","_",levpart)
    pctpresent <- T
  }
  return(list(date=date,sdate=sdate,val=val,sval=sval,pctpresent=pctpresent))
}

addStepDateToVek <- function(dates,idx,dtvek,vvek){
  cd1 <- crackdate(dates[[idx]])
  if (!cd1$pctpresent){
    # if there is no % don't do anything
    return(vvek)
  }
  dt1 <- cd1$date
  if (idx<length(dates)){
    cd2 <- crackdate(dates[[idx+1]])
    dt2 <- cd2$date
  } else {
    dt2 <- max(dtvek)
  }
  #print("addstepdate")
  #print(dt1)
  #print(dt2)
  val <- cd1$val
  tochg <-  dt1<=dtvek & dtvek<= dt2
  vvek[ tochg ] <- val
  #print(sprintf("changed %d values to %d",sum(tochg),val))
  return(vvek)
}
getStepDates <- function(dates,dtvek){
 vvek <- rep(0,length(dtvek))
 for (i in 1:length(dates)){
   vvek <- addStepDateToVek(dates,i,dtvek,vvek)
 }
 return(vvek)
}
getSmcStepDates <- function(dtvek){
 return(getStepDates(smcdates,dtvek))
}
getXabStepDates <- function(dtvek){
 return(getStepDates(xbxdates,dtvek))
}

Plot functions

addVlinesAndText <- function(vlines,gp){
  if (is.null(vlines)) return(gp) # do nothing in this case
  
  # split the lines and convert to data.frame
  sar <- strsplit(vlines,"/")
  # the following reforms the date strings into a data.frame for geom_vline
  ldf <- data.frame(t(matrix(unlist(sar),length(sar[[1]]),length(sar)))) #tricky
  names(ldf) <- c("dt","clr","lab")
  ldf$dt <- s2date(ldf$dt)
  ldf$ndt <- as.numeric(ldf$dt)
  # add a newline to the front so as to display the text 
  # this keeps the text from writing on top of the vline
  
  ldf$lab <- paste0("\n",ldf$lab) 
  # now actually add the verticle lines and the text 
  gp <- gp + geom_vline(xintercept=ldf$ndt,color=ldf$clr) +
             annotate(geom="text",x=ldf$dt,y=0,label=ldf$lab,color=ldf$clr,hjust=0,angle=90,na.rm=T)
  return(gp)
}
addBackground <- function(backg,gp){
  if (is.null(backg)) return(gp) # do nothing in this case
  gp <- gp + theme(panel.background = element_rect(fill = backg))
  return(gp)
}
overdate <- function(ovdate,defdate){
  # date override
  rv <- defdate
  if (!is.null(ovdate)) { 
    rv <- ovdate
  }
  return(rv)
}
dailyplot <- function(ddf,x,y,mtit="",xlab="date",ylab=NULL,vlines=NULL,backg=NULL,series=NULL,ovsdate=NULL,ovedate=NULL){
  # Single series plot  with monthly breaks on the x-axis
  
  # override dates if needed
  dpsdate <- overdate(ovsdate,sdate)
  dpedate <- overdate(ovedate,edate)
  
  brkctrl <- "1 month"
  dltdays <- difftime(dpedate,dpsdate,"days") 
  if (dltdays<30) brkctrl <- "1 day"
  
  gp <- ggplot(ddf,aes_string(x=x,y=y)) + 
             geom_line(aes_string(color=series),na.rm=T)  +
             xlab(xlab) + ylab(ylab) + ggtitle(mtit) +
             scale_x_datetime("Date",breaks = date_breaks(brkctrl),limits=c(dpsdate,dpedate))

  gp <- addVlinesAndText(vlines,gp)

  gp <- addBackground(backg,gp)

  return(gp)
}
residplot <- function(ddf,x,y,mtit="",xlab="date",ylab=NULL,vlines=NULL,backg=NULL,series=NULL,ovsdate=NULL,ovedate=NULL){
  # Single series plot  with monthly breaks on the x-axis
  
  # override dates if needed
  dpsdate <- overdate(ovsdate,sdate)
  dpedate <- overdate(ovedate,edate)
  
  brkctrl <- "1 month"
  dltdays <- difftime(dpedate,dpsdate,"days") 
  if (dltdays<30) brkctrl <- "1 day"
  
  gp <- ggplot(ddf,aes_string(x=x,y=y)) + 
             geom_point(aes_string(color=series),na.rm=T)  +
             xlab(xlab) + ylab(ylab) + ggtitle(mtit) +
             scale_x_datetime("Date",breaks = date_breaks(brkctrl),limits=c(dpsdate,dpedate)) +
             theme(axis.text.x = element_text(angle = 30, hjust = 1))

  gp <- addVlinesAndText(vlines,gp)

  gp <- addBackground(backg,gp)
  
  hp <- ggplot(ddf) + geom_histogram(aes_string(x=y),bins=30)
  hp <- addBackground(backg,hp)
  
  ghp <- grid.arrange(gp, hp, ncol=2,widths = c(2,1))

  return(ghp)
}

Read in consolidated chat, call, and session volume data

stload <- Sys.time()
tfname <- sprintf("%s/%s",fpath,"colsolidatedTorontoData01.csv")
condf <- read.csv(tfname)

minsessfilt <- 5000
nbef <- nrow(condf)
condf <- condf %>% filter(minsessfilt<actsess)
naft <- nrow(condf)
print(sprintf("Filtered %d of %d hours because sessions less than %d",(nbef-naft),naft,minsessfilt))
## [1] "Filtered 9 of 5815 hours because sessions less than 5000"
condf <- condf %>% mutate( dt = as.POSIXct(dt,tz=tztz) ) %>%
                   mutate( log_winchib = log(winchib) ) %>%
                   mutate( log_wincall = log(wincall) ) %>%
                   mutate( log_xbxchib = log(xbxchib) ) %>%
                   mutate( log_xbxcall = log(xbxcall) ) %>%
                   mutate( rate_winchib = winchib/actsess ) %>%
                   mutate( rate_wincall = wincall/actsess ) %>%
                   mutate( rate_xbxchib = xbxchib/actsess ) %>%
                   mutate( rate_xbxcall = xbxcall/actsess ) 

dcondf <- condf %>% group_by(dnum) %>% summarise(dt=min(dt),
                                                 totchib=sum(totchib),winchib=sum(winchib),xbxchib=sum(xbxchib),
                                                 totcall=sum(totcall),wincall=sum(wincall),xbxcall=sum(xbxcall),
                                                 actsess=sum(actsess),actuser=sum(actuser)
                                                 )

elap <- as.numeric((Sys.time()-stload)[1],units="secs")
print(sprintf("Loading consolidated data took %.1f secs",elap))
## [1] "Loading consolidated data took 0.1 secs"

chats

if (verbose>=2){
  pltdf <- dcondf %>% gather(series,chib,-dt) %>% filter(series %in% c("totchib","winchib","xbxchib"))
  dailyplot(pltdf,"dt","chib",series="series",mtit="Chats In Block",ylab="Sum",vlines=totdates,backg=totback)
}

# calls

if (verbose>=2){
  pltdf <- dcondf %>% gather(series,call,-dt) %>% filter(series %in% c("totcall","wincall","xbxcall"))
  dailyplot(pltdf,"dt","call",series="series",mtit="Calls",ylab="Sum",vlines=totdates,backg=totback)
}

# sessions

if (verbose>=2){
  pltdf <- dcondf %>% gather(series,active,-dt) %>% filter(series %in% c("actsess","actuser"))
  dailyplot(pltdf,"dt","active",series="series",mtit="Sessions",ylab="Sum",vlines=totdates,backg=totback)
}

regression code

results1 <- list()
results2 <- list()
initResults <- function(){
  results1 <<- list()
  results2 <<- list()
}
  
addToResults1 <- function(newresults){
  results1[[length(results1)+1]] <<- newresults
}
addToResults2 <- function(newresults){
  results2[[length(results2)+1]] <<- newresults
}
listtodf <- function(lst){
  nr <- length(lst)    # rows
  if (nr==0) return(NULL)
  nvek <- names(lst[[1]])
  nc <- length(nvek)   # columns
  df <- data.frame(idx=1:nr) # preallocate length
  for (i in 1:nc){
    iname <- nvek[i]
    df[[iname]] <- sapply(lst,`[[`,iname)
  }
  return(df)
}
getResults1 <- function(){
  return(listtodf(results1))
}
getResults2 <- function(){
  return(listtodf(results2))
}

getregdf <- function(df,vname,discdate,befdays,aftdays,dates){
  
  # filter on time
  sregdate <- discdate - days(befdays)
  eregdate <- discdate + days(aftdays)
  df <- df %>% filter( sregdate <= dt & dt <= eregdate) 
  
  df <- df %>% mutate( hour=as.factor(hour(dt)) ) %>%
               mutate( dow=as.factor(wday(dt))) %>%  
               mutate( idx=1:nrow(df) ) 

  # add level variables
  cutidx <- which(df$dt==discdate)

    df <- df %>% mutate( lin1=if_else(idx<cutidx,idx-cutidx,0L)) %>%
                 mutate( lin2=if_else(idx<cutidx,0L,idx-cutidx)) %>%
                 mutate( lin0=1 ) %>%
                 mutate( postchange=if_else(idx<cutidx,0L,1L) )

  return(df)
}
formlist = list(
  formel1="%s ~ hour + dow + lin1 + lin2 + postchange",
  formel2a="%s ~ hour * dow - 1",
  formel2aa="%s ~ hour + dow",
  formel2b="%s ~ lin1 + lin2 + postchange"
#  formel2b="%s ~ lin0 + postchange"
)
doregression <- function(df,vname,formel,chgdate,befdays,aftdays){
  df <- getregdf(df,vname,chgdate,befdays,aftdays,dates)
  formstr <- sprintf(formlist[[formel]],vname)
  form <- as.formula(formstr)
  fit <- glm(form,data=df)
  summary_fit <- summary(fit)
  coeftest_fit <- coeftest(fit, vcov = vcovHC(fit, "HC1"))
  if (verbose>0){
    print(summary_fit)
    print(coeftest_fit)
  }
  df$resid <- resid(fit)
  df$predicted <-  df[[vname]]-df$resid
  # for debugging
  ffit <<- fit
  sfit <<- summary_fit
  cfit <<- coeftest_fit
  return(list(df=df,formstr=formstr,fit=fit,sfit=summary_fit,cfit=coeftest_fit))
}
regress1 <- function(df,vname,area,backg,dates,chgdate,mtit1,befdays=7,aftdays=7,model="model1"){
  crk <- crackdate(chgdate)
  idx <- length(results1)+1
  
  print(sprintf("%d Regression for %s %s on %s",idx,vname,model,chgdate))

  rv <- doregression(df,vname,"formel1",crk$date,befdays,aftdays)
  if (verbose>0){
    df1 <- rv$df %>%  gather_("series",vname,c(vname,"predicted"))
    mtit <- sprintf("%d - %s - %s",idx,mtit1,rv$formstr)
    plt <- dailyplot(df1,"dt",vname,series="series",mtit,ylab="Sum",
                     vlines=dates,backg=backg,ovsdate=min(df1$dt),ovedate=max(df1$dt))
    print(plt)
    mtit <- sprintf("%d - Residuals - %s - %s",idx,mtit1,rv$formstr)
    plt <- residplot(df1,"dt","resid",mtit,ylab="Sum",
                     vlines=dates,backg=backg,ovsdate=min(df1$dt),ovedate=max(df1$dt))
  }
  newresults <- list(var=vname,
                     title=mtit1,
                     chgDate=crk$sdate,
                     chgVal=crk$sval,
                     model=model,
                     formula=rv$formstr,
                     chgCoef.fit=rv$fit$coefficients["postchange"],
                     chgCoef.smry=rv$sfit$coefficients["postchange",1],
                     chgStd.smry=rv$sfit$coefficients["postchange",2],
                     chgTval.smry=rv$sfit$coefficients["postchange",3],
                     chgPval.smry=rv$sfit$coefficients["postchange",4],
                     chgZval.cfit=rv2$cfit["postchange","z value"],
                     chgPval.cfit=rv2$cfit["postchange","Pr(>|z|)"]                     
                     )
  addToResults1(newresults)
}
regress2 <- function(df,vname,area,backg,dates,chgdate,mtit1,befdays=7,aftdays=7,model="model2c"){
  crk <- crackdate(chgdate)
  idx <- length(results1)+1
  
  print(sprintf("%d Regression for %s %s on %s",idx,vname,model,chgdate))

  # step 1
  rv1 <- doregression(df,vname,"formel2aa",crk$date,befdays,aftdays)
  if (verbose>0){
    df1 <- rv1$df %>%  gather_("series",vname,c(vname,"predicted"))
    mtit <- sprintf("%d - %s - %s",idx,mtit1,rv1$formstr)
    plt <- dailyplot(df1,"dt",vname,series="series",mtit,ylab="Sum",
                     vlines=dates,backg=backg,ovsdate=min(df1$dt),ovedate=max(df1$dt))
    print(plt)
    mtit <- sprintf("%d - Step 1 Residuals - %s - %s",idx,mtit1,rv1$formstr)
    plt <- residplot(df1,"dt","resid",mtit,ylab="Sum",
                     vlines=dates,backg=backg,ovsdate=min(df1$dt),ovedate=max(df1$dt))
  }
  newresults <- list(var=vname,
                     title=mtit1,
                     chgDate=crk$sdate,
                     chgVal=crk$sval,
                     model=model,
                     formula=rv1$formstr,
                     aic=rv1$fit$aic,
                     deviance=rv1$fit$dev
                     )
  addToResults1(newresults)
  
  # step 2
  df2 <- rv1$df
  df2$step1residuals <- rv1$fit$residuals
  vname2 <- "step1residuals"
  rv2 <- doregression(df2,vname2,"formel2b",crk$date,7,7)
  if (verbose>0){
    df1 <- rv2$df %>%  gather_("series",vname2,c(vname2,"predicted"))
    mtit <- sprintf("%d - %s - %s",idx,mtit1,rv2$formstr)
    plt <- dailyplot(df1,"dt",vname2,series="series",mtit,ylab="Sum",
                     vlines=dates,backg=backg,ovsdate=min(df1$dt),ovedate=max(df1$dt))
    print(plt)
    mtit <- sprintf("%d - Step 2 Residuals - %s - %s",idx,mtit1,rv2$formstr)
    plt <- residplot(df1,"dt","resid",mtit,ylab="Sum",
                     vlines=dates,backg=backg,ovsdate=min(df1$dt),ovedate=max(df1$dt))
  }

  newresults <- list(var=vname,
                     title=mtit1,
                     chgDate=crk$sdate,
                     chgVal=crk$sval,
                     model=model,
                     formula=rv2$formstr,
                     chgCoef.fit=rv2$fit$coefficients["postchange"],
                     chgCoef.smry=rv2$sfit$coefficients["postchange",1],
                     chgStd.smry=rv2$sfit$coefficients["postchange",2],
                     chgTval.smry=rv2$sfit$coefficients["postchange",3],
                     chgPval.smry=rv2$sfit$coefficients["postchange",4],
                     chgZval.cfit=rv2$cfit["postchange","z value"],
                     chgPval.cfit=rv2$cfit["postchange","Pr(>|z|)"]                     
                     )
  addToResults2(newresults)
}
regressoverdates <- function(model,df,vname,area,title,befdays=7,aftdays=7){
  if (model=="model2" & grepl("^log",vname)) return(NULL)  # model2 is pointless for logs
  if (area=="SMC"){
    dates <- smcdates
    backg <- smcback
  } else {
    dates <- xbxdates
    backg <- xbxback
  }
  atit <- sprintf("%s - %s",area,title)
  for (chgdate in dates){
    if (model=="model1"){
      regress1(condf,vname,area,backg,dates,chgdate,atit,befdays=7,aftdays=7) 
    } else {
      regress2(condf,vname,area,backg,dates,chgdate,atit,befdays=7,aftdays=7) 
    }
  }
}
initResults()

justdoone <- F # for testing
verbose <- 2

m <- "model2"
if (justdoone){
  regressoverdates(m,condf,"winchib","SMC","Hourly Chats")
#  regressoverdates(m,condf,"winchib","Xbox","Hourly Chats")
} else {
  regressoverdates(m,condf,"winchib","SMC","Hourly Chats")
  regressoverdates(m,condf,"log_winchib","SMC","log(Hourly Chats)")
  regressoverdates(m,condf,"rate_winchib","SMC","Hourly Chat Rate per Session")
  
  regressoverdates(m,condf,"wincall","SMC","Hourly Chats")
  regressoverdates(m,condf,"log_wincall","SMC","log(Hourly Chats)")
  regressoverdates(m,condf,"rate_wincall","SMC","Hourly Chat Rate per Session")
  
  regressoverdates(m,condf,"xbxchib","Xbox","Hourly Chats")
  regressoverdates(m,condf,"log_xbxchib","Xbox","log(Hourly Chats)")
  regressoverdates(m,condf,"rate_xbxchib","Xbox","Hourly Chat Rate per Session")
  
  regressoverdates(m,condf,"xbxcall","Xbox","Hourly Chats")
  regressoverdates(m,condf,"log_xbxcall","Xbox","log(Hourly Chats)")
  regressoverdates(m,condf,"rate_xbxcall","Xbox","Hourly Chat Rate per Session")
}
## [1] "1 Regression for winchib model2c on 2016-08-17/red/0-30%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -194.26   -27.77    -1.19    31.93   110.49  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  315.520     13.553  23.280  < 2e-16 ***
## hour1        -35.223     17.281  -2.038 0.042387 *  
## hour2        -76.366     17.281  -4.419 1.38e-05 ***
## hour3       -101.223     17.281  -5.857 1.21e-08 ***
## hour4       -120.937     17.281  -6.998 1.64e-11 ***
## hour5       -159.937     17.281  -9.255  < 2e-16 ***
## hour6       -186.366     17.281 -10.784  < 2e-16 ***
## hour7       -193.509     17.281 -11.197  < 2e-16 ***
## hour8       -195.152     17.281 -11.293  < 2e-16 ***
## hour9       -202.294     17.281 -11.706  < 2e-16 ***
## hour10      -203.580     17.281 -11.780  < 2e-16 ***
## hour11      -155.866     17.281  -9.019  < 2e-16 ***
## hour12       -99.009     17.281  -5.729 2.41e-08 ***
## hour13       -43.437     17.281  -2.514 0.012465 *  
## hour14        38.206     17.281   2.211 0.027788 *  
## hour15        30.706     17.281   1.777 0.076594 .  
## hour16        60.634     17.281   3.509 0.000518 ***
## hour17        88.706     17.281   5.133 5.07e-07 ***
## hour18        96.348     17.281   5.575 5.42e-08 ***
## hour19        40.563     17.281   2.347 0.019552 *  
## hour20        10.634     17.281   0.615 0.538783    
## hour21         4.491     17.281   0.260 0.795126    
## hour22        31.063     17.281   1.797 0.073247 .  
## hour23        -1.437     17.281  -0.083 0.933768    
## dow2          74.250      9.490   7.824 8.31e-14 ***
## dow3          98.000      9.490  10.327  < 2e-16 ***
## dow4         108.357      9.444  11.473  < 2e-16 ***
## dow5          67.563      9.490   7.119 7.72e-12 ***
## dow6          43.062      9.490   4.538 8.17e-06 ***
## dow7          -9.812      9.490  -1.034 0.301955    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 2161.416)
## 
##     Null deviance: 4550168  on 336  degrees of freedom
## Residual deviance:  663555  on 307  degrees of freedom
## AIC: 3574.6
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  315.5203    13.1727  23.9526 < 2.2e-16 ***
## hour1        -35.2231    17.2966  -2.0364  0.041709 *  
## hour2        -76.3659    18.4023  -4.1498 3.328e-05 ***
## hour3       -101.2231    16.3329  -6.1975 5.737e-10 ***
## hour4       -120.9374    14.7398  -8.2048 2.310e-16 ***
## hour5       -159.9374    14.8951 -10.7376 < 2.2e-16 ***
## hour6       -186.3659    18.2887 -10.1902 < 2.2e-16 ***
## hour7       -193.5088    22.3161  -8.6712 < 2.2e-16 ***
## hour8       -195.1516    17.0854 -11.4221 < 2.2e-16 ***
## hour9       -202.2945    13.9144 -14.5385 < 2.2e-16 ***
## hour10      -203.5802    14.2608 -14.2755 < 2.2e-16 ***
## hour11      -155.8659    13.9912 -11.1403 < 2.2e-16 ***
## hour12       -99.0088    12.6895  -7.8024 6.074e-15 ***
## hour13       -43.4374    15.3326  -2.8330  0.004611 ** 
## hour14        38.2055    18.0979   2.1110  0.034768 *  
## hour15        30.7055    20.7421   1.4803  0.138781    
## hour16        60.6341    15.4081   3.9352 8.312e-05 ***
## hour17        88.7055    20.0789   4.4178 9.969e-06 ***
## hour18        96.3484    17.9696   5.3618 8.242e-08 ***
## hour19        40.5626    18.4752   2.1955  0.028126 *  
## hour20        10.6341    17.4406   0.6097  0.542040    
## hour21         4.4912    19.1020   0.2351  0.814117    
## hour22        31.0626    17.6692   1.7580  0.078746 .  
## hour23        -1.4374    13.8910  -0.1035  0.917587    
## dow2          74.2500    10.5547   7.0348 1.995e-12 ***
## dow3          98.0000     8.6167  11.3733 < 2.2e-16 ***
## dow4         108.3567     9.8739  10.9740 < 2.2e-16 ***
## dow5          67.5625     8.5319   7.9188 2.399e-15 ***
## dow6          43.0625     8.4897   5.0723 3.930e-07 ***
## dow7          -9.8125    11.5368  -0.8505  0.395025    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -188.27   -26.72    -0.82    30.42   104.31  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)
## (Intercept) -10.66729    6.87059  -1.553    0.121
## lin1         -0.11394    0.07052  -1.616    0.107
## lin2          0.09631    0.06989   1.378    0.169
## postchange    3.60991    9.65933   0.374    0.709
## 
## (Dispersion parameter for gaussian family taken to be 1964.961)
## 
##     Null deviance: 663555  on 336  degrees of freedom
## Residual deviance: 654332  on 333  degrees of freedom
## AIC: 3517.9
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error z value Pr(>|z|)  
## (Intercept) -10.667289   7.851414 -1.3586  0.17426  
## lin1         -0.113944   0.075287 -1.5135  0.13016  
## lin2          0.096313   0.056630  1.7008  0.08899 .
## postchange    3.609910   9.569332  0.3772  0.70600  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "2 Regression for winchib model2c on 2016-09-1/red/30-50%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -222.590   -25.908     0.954    27.259    83.125  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  309.0282    11.9855  25.784  < 2e-16 ***
## hour1          0.7221    15.2827   0.047 0.962345    
## hour2        -47.2065    15.2827  -3.089 0.002193 ** 
## hour3        -99.8493    15.2827  -6.534 2.67e-10 ***
## hour4       -119.9922    15.2827  -7.852 6.93e-14 ***
## hour5       -151.2065    15.2827  -9.894  < 2e-16 ***
## hour6       -175.8493    15.2827 -11.506  < 2e-16 ***
## hour7       -181.0636    15.2827 -11.848  < 2e-16 ***
## hour8       -183.2779    15.2827 -11.993  < 2e-16 ***
## hour9       -205.4208    15.2827 -13.441  < 2e-16 ***
## hour10      -183.7065    15.2827 -12.021  < 2e-16 ***
## hour11      -153.9208    15.2827 -10.072  < 2e-16 ***
## hour12       -77.1351    15.2827  -5.047 7.69e-07 ***
## hour13       -13.8493    15.2827  -0.906 0.365534    
## hour14        44.2221    15.2827   2.894 0.004081 ** 
## hour15        56.0078    15.2827   3.665 0.000292 ***
## hour16        61.6507    15.2827   4.034 6.93e-05 ***
## hour17        97.0792    15.2827   6.352 7.64e-10 ***
## hour18        97.9364    15.2827   6.408 5.53e-10 ***
## hour19        45.2221    15.2827   2.959 0.003326 ** 
## hour20        32.1507    15.2827   2.104 0.036215 *  
## hour21        25.2935    15.2827   1.655 0.098937 .  
## hour22         8.5792    15.2827   0.561 0.574955    
## hour23        -5.0636    15.2827  -0.331 0.740620    
## dow2          69.7292     8.3923   8.309 3.15e-15 ***
## dow3          72.6250     8.3923   8.654 2.85e-16 ***
## dow4          85.5625     8.3923  10.195  < 2e-16 ***
## dow5          91.0811     8.3521  10.905  < 2e-16 ***
## dow6          54.7708     8.3923   6.526 2.78e-10 ***
## dow7         -11.0208     8.3923  -1.313 0.190094    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1690.34)
## 
##     Null deviance: 4246614  on 336  degrees of freedom
## Residual deviance:  518934  on 307  degrees of freedom
## AIC: 3491.8
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  309.02823   11.25528  27.4563 < 2.2e-16 ***
## hour1          0.72209   13.17371   0.0548 0.9562874    
## hour2        -47.20648   15.64756  -3.0169 0.0025541 ** 
## hour3        -99.84934   13.51356  -7.3888 1.481e-13 ***
## hour4       -119.99219   13.53227  -8.8671 < 2.2e-16 ***
## hour5       -151.20648   13.48824 -11.2102 < 2.2e-16 ***
## hour6       -175.84934   14.09367 -12.4772 < 2.2e-16 ***
## hour7       -181.06362   12.94573 -13.9864 < 2.2e-16 ***
## hour8       -183.27791   12.95992 -14.1419 < 2.2e-16 ***
## hour9       -205.42077   13.69590 -14.9987 < 2.2e-16 ***
## hour10      -183.70648   13.81014 -13.3023 < 2.2e-16 ***
## hour11      -153.92077   12.06867 -12.7537 < 2.2e-16 ***
## hour12       -77.13505   11.33145  -6.8072 9.954e-12 ***
## hour13       -13.84934   12.20614  -1.1346 0.2565345    
## hour14        44.22209   15.22184   2.9052 0.0036705 ** 
## hour15        56.00781   16.24107   3.4485 0.0005637 ***
## hour16        61.65066   13.02913   4.7318 2.226e-06 ***
## hour17        97.07923   15.53487   6.2491 4.128e-10 ***
## hour18        97.93638   14.48562   6.7609 1.371e-11 ***
## hour19        45.22209   14.46154   3.1271 0.0017656 ** 
## hour20        32.15066   18.17263   1.7692 0.0768638 .  
## hour21        25.29352   18.77521   1.3472 0.1779235    
## hour22         8.57923   15.04795   0.5701 0.5685919    
## hour23        -5.06362   21.67395  -0.2336 0.8152745    
## dow2          69.72917    7.79336   8.9473 < 2.2e-16 ***
## dow3          72.62500    9.10873   7.9731 1.547e-15 ***
## dow4          85.56250    7.26739  11.7735 < 2.2e-16 ***
## dow5          91.08106    8.01940  11.3576 < 2.2e-16 ***
## dow6          54.77083    7.21107   7.5954 3.069e-14 ***
## dow7         -11.02083    9.55160  -1.1538 0.2485736    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -220.743   -25.763     2.223    25.617    84.395  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)
## (Intercept)  4.74664    6.08677   0.780    0.436
## lin1         0.07913    0.06247   1.267    0.206
## lin2        -0.06399    0.06192  -1.033    0.302
## postchange   2.55692    8.55736   0.299    0.765
## 
## (Dispersion parameter for gaussian family taken to be 1542.199)
## 
##     Null deviance: 518934  on 336  degrees of freedom
## Residual deviance: 513552  on 333  degrees of freedom
## AIC: 3436.2
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error z value Pr(>|z|)  
## (Intercept)  4.746641   4.549205  1.0434  0.29676  
## lin1         0.079134   0.046540  1.7003  0.08907 .
## lin2        -0.063987   0.069857 -0.9160  0.35968  
## postchange   2.556921   7.725168  0.3310  0.74066  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "3 Regression for winchib model2c on 2016-09-07/red/50-100%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -208.801   -23.964     1.464    26.363    83.771  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  300.369     11.302  26.577  < 2e-16 ***
## hour1          7.071     14.411   0.491 0.624007    
## hour2        -51.143     14.411  -3.549 0.000447 ***
## hour3        -92.858     14.411  -6.444 4.51e-10 ***
## hour4       -115.143     14.411  -7.990 2.75e-14 ***
## hour5       -155.715     14.411 -10.805  < 2e-16 ***
## hour6       -179.429     14.411 -12.451  < 2e-16 ***
## hour7       -188.572     14.411 -13.085  < 2e-16 ***
## hour8       -184.429     14.411 -12.798  < 2e-16 ***
## hour9       -196.858     14.411 -13.660  < 2e-16 ***
## hour10      -180.786     14.411 -12.545  < 2e-16 ***
## hour11      -149.715     14.411 -10.389  < 2e-16 ***
## hour12       -79.858     14.411  -5.541 6.46e-08 ***
## hour13       -12.429     14.411  -0.862 0.389101    
## hour14        56.214     14.411   3.901 0.000118 ***
## hour15        65.071     14.411   4.515 9.02e-06 ***
## hour16        66.642     14.411   4.624 5.54e-06 ***
## hour17       106.357     14.411   7.380 1.49e-12 ***
## hour18        96.357     14.411   6.686 1.08e-10 ***
## hour19        44.928     14.411   3.118 0.001996 ** 
## hour20        40.785     14.411   2.830 0.004959 ** 
## hour21        34.500     14.411   2.394 0.017266 *  
## hour22        23.857     14.411   1.655 0.098851 .  
## hour23        -6.215     14.411  -0.431 0.666589    
## dow2          59.563      7.914   7.527 5.82e-13 ***
## dow3          68.646      7.914   8.674 2.46e-16 ***
## dow4          90.625      7.876  11.507  < 2e-16 ***
## dow5          95.292      7.914  12.042  < 2e-16 ***
## dow6          60.167      7.914   7.603 3.55e-13 ***
## dow7          -2.875      7.914  -0.363 0.716630    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1502.986)
## 
##     Null deviance: 4321974  on 336  degrees of freedom
## Residual deviance:  461417  on 307  degrees of freedom
## AIC: 3452.2
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  300.3694    10.1497  29.5939 < 2.2e-16 ***
## hour1          7.0710    11.8261   0.5979 0.5498945    
## hour2        -51.1432    13.6974  -3.7338 0.0001886 ***
## hour3        -92.8575    11.3703  -8.1667 3.169e-16 ***
## hour4       -115.1432    11.9472  -9.6377 < 2.2e-16 ***
## hour5       -155.7147    12.0199 -12.9548 < 2.2e-16 ***
## hour6       -179.4290    13.2955 -13.4954 < 2.2e-16 ***
## hour7       -188.5718    14.1546 -13.3223 < 2.2e-16 ***
## hour8       -184.4290    11.8329 -15.5861 < 2.2e-16 ***
## hour9       -196.8575    11.4235 -17.2327 < 2.2e-16 ***
## hour10      -180.7861    12.5485 -14.4070 < 2.2e-16 ***
## hour11      -149.7147    11.4620 -13.0619 < 2.2e-16 ***
## hour12       -79.8575    11.6524  -6.8533 7.215e-12 ***
## hour13       -12.4290    11.3165  -1.0983 0.2720721    
## hour14        56.2139    14.6461   3.8381 0.0001240 ***
## hour15        65.0710    13.9063   4.6793 2.879e-06 ***
## hour16        66.6425    12.2842   5.4251 5.793e-08 ***
## hour17       106.3568    14.2934   7.4410 9.993e-14 ***
## hour18        96.3568    14.7365   6.5386 6.208e-11 ***
## hour19        44.9282    13.7526   3.2669 0.0010874 ** 
## hour20        40.7853    16.0758   2.5371 0.0111786 *  
## hour21        34.4996    16.2322   2.1254 0.0335545 *  
## hour22        23.8568    14.8182   1.6100 0.1074058    
## hour23        -6.2147    19.8169  -0.3136 0.7538207    
## dow2          59.5625     6.8048   8.7530 < 2.2e-16 ***
## dow3          68.6458     8.7517   7.8437 4.374e-15 ***
## dow4          90.6251     6.7297  13.4665 < 2.2e-16 ***
## dow5          95.2917     7.8022  12.2134 < 2.2e-16 ***
## dow6          60.1667     6.6020   9.1134 < 2.2e-16 ***
## dow7          -2.8750     8.8986  -0.3231 0.7466315    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -212.757   -25.379     0.875    27.996    79.898  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)
## (Intercept)  3.965423   5.747449   0.690    0.491
## lin1         0.009277   0.058992   0.157    0.875
## lin2        -0.015965   0.058469  -0.273    0.785
## postchange  -5.787096   8.080309  -0.716    0.474
## 
## (Dispersion parameter for gaussian family taken to be 1375.043)
## 
##     Null deviance: 461417  on 336  degrees of freedom
## Residual deviance: 457889  on 333  degrees of freedom
## AIC: 3397.6
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error z value Pr(>|z|)
## (Intercept)  3.9654232  7.9033773  0.5017   0.6159
## lin1         0.0092768  0.0721808  0.1285   0.8977
## lin2        -0.0159647  0.0485488 -0.3288   0.7423
## postchange  -5.7870962  9.3365365 -0.6198   0.5354

## [1] "4 Regression for rate_winchib model2c on 2016-08-17/red/0-30%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -1.728e-03  -1.971e-04  -9.260e-06   2.225e-04   1.197e-03  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.240e-03  1.071e-04  39.595  < 2e-16 ***
## hour1       -3.475e-04  1.366e-04  -2.545   0.0114 *  
## hour2       -5.469e-04  1.366e-04  -4.005 7.78e-05 ***
## hour3       -6.661e-04  1.366e-04  -4.878 1.72e-06 ***
## hour4       -6.951e-04  1.366e-04  -5.090 6.25e-07 ***
## hour5       -1.063e-03  1.366e-04  -7.787 1.06e-13 ***
## hour6       -1.453e-03  1.366e-04 -10.638  < 2e-16 ***
## hour7       -1.564e-03  1.366e-04 -11.455  < 2e-16 ***
## hour8       -1.589e-03  1.366e-04 -11.637  < 2e-16 ***
## hour9       -1.742e-03  1.366e-04 -12.754  < 2e-16 ***
## hour10      -1.928e-03  1.366e-04 -14.121  < 2e-16 ***
## hour11      -1.949e-03  1.366e-04 -14.274  < 2e-16 ***
## hour12      -2.010e-03  1.366e-04 -14.722  < 2e-16 ***
## hour13      -1.899e-03  1.366e-04 -13.911  < 2e-16 ***
## hour14      -1.530e-03  1.366e-04 -11.205  < 2e-16 ***
## hour15      -1.641e-03  1.366e-04 -12.019  < 2e-16 ***
## hour16      -1.268e-03  1.366e-04  -9.285  < 2e-16 ***
## hour17      -1.068e-03  1.366e-04  -7.819 8.62e-14 ***
## hour18      -9.498e-04  1.366e-04  -6.956 2.12e-11 ***
## hour19      -1.200e-03  1.366e-04  -8.790  < 2e-16 ***
## hour20      -1.197e-03  1.366e-04  -8.770  < 2e-16 ***
## hour21      -8.861e-04  1.366e-04  -6.490 3.45e-10 ***
## hour22      -1.456e-04  1.366e-04  -1.066   0.2870    
## hour23      -9.847e-05  1.366e-04  -0.721   0.4714    
## dow2        -9.256e-04  7.498e-05 -12.344  < 2e-16 ***
## dow3        -5.589e-04  7.498e-05  -7.454 9.31e-13 ***
## dow4        -4.348e-04  7.462e-05  -5.827 1.43e-08 ***
## dow5        -8.865e-04  7.498e-05 -11.823  < 2e-16 ***
## dow6        -7.431e-04  7.498e-05  -9.911  < 2e-16 ***
## dow7        -1.459e-04  7.498e-05  -1.945   0.0526 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1.349366e-07)
## 
##     Null deviance: 2.0249e-04  on 336  degrees of freedom
## Residual deviance: 4.1426e-05  on 307  degrees of freedom
## AIC: -4343.9
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  4.2401e-03  1.2502e-04  33.9162 < 2.2e-16 ***
## hour1       -3.4745e-04  1.6008e-04  -2.1705 0.0299685 *  
## hour2       -5.4691e-04  1.5146e-04  -3.6108 0.0003052 ***
## hour3       -6.6612e-04  1.2273e-04  -5.4274 5.719e-08 ***
## hour4       -6.9505e-04  1.4681e-04  -4.7343 2.198e-06 ***
## hour5       -1.0633e-03  1.1986e-04  -8.8708 < 2.2e-16 ***
## hour6       -1.4526e-03  1.3886e-04 -10.4608 < 2.2e-16 ***
## hour7       -1.5641e-03  1.9025e-04  -8.2213 < 2.2e-16 ***
## hour8       -1.5890e-03  1.2861e-04 -12.3551 < 2.2e-16 ***
## hour9       -1.7416e-03  1.1759e-04 -14.8107 < 2.2e-16 ***
## hour10      -1.9281e-03  1.2701e-04 -15.1808 < 2.2e-16 ***
## hour11      -1.9491e-03  1.2802e-04 -15.2254 < 2.2e-16 ***
## hour12      -2.0103e-03  1.3610e-04 -14.7707 < 2.2e-16 ***
## hour13      -1.8995e-03  1.6382e-04 -11.5952 < 2.2e-16 ***
## hour14      -1.5300e-03  1.6271e-04  -9.4036 < 2.2e-16 ***
## hour15      -1.6412e-03  1.2910e-04 -12.7126 < 2.2e-16 ***
## hour16      -1.2678e-03  1.1230e-04 -11.2893 < 2.2e-16 ***
## hour17      -1.0676e-03  1.3550e-04  -7.8790 3.299e-15 ***
## hour18      -9.4978e-04  1.3589e-04  -6.9895 2.758e-12 ***
## hour19      -1.2003e-03  1.1966e-04 -10.0310 < 2.2e-16 ***
## hour20      -1.1975e-03  1.2578e-04  -9.5202 < 2.2e-16 ***
## hour21      -8.8612e-04  1.6095e-04  -5.5058 3.676e-08 ***
## hour22      -1.4563e-04  1.6741e-04  -0.8698 0.3843832    
## hour23      -9.8472e-05  1.4981e-04  -0.6573 0.5109893    
## dow2        -9.2556e-04  9.6402e-05  -9.6010 < 2.2e-16 ***
## dow3        -5.5891e-04  8.0319e-05  -6.9586 3.437e-12 ***
## dow4        -4.3481e-04  8.2618e-05  -5.2629 1.418e-07 ***
## dow5        -8.8654e-04  8.3923e-05 -10.5637 < 2.2e-16 ***
## dow6        -7.4313e-04  8.5198e-05  -8.7223 < 2.2e-16 ***
## dow7        -1.4587e-04  1.0095e-04  -1.4449 0.1484731    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -1.667e-03  -2.078e-04  -1.730e-06   2.158e-04   1.253e-03  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)
## (Intercept) -7.108e-05  5.408e-05  -1.314    0.190
## lin1        -2.458e-07  5.551e-07  -0.443    0.658
## lin2         1.393e-07  5.502e-07   0.253    0.800
## postchange   1.094e-04  7.604e-05   1.439    0.151
## 
## (Dispersion parameter for gaussian family taken to be 1.217596e-07)
## 
##     Null deviance: 4.1426e-05  on 336  degrees of freedom
## Residual deviance: 4.0546e-05  on 333  degrees of freedom
## AIC: -4403.1
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)
## (Intercept) -7.1078e-05  5.9081e-05 -1.2031   0.2290
## lin1        -2.4578e-07  5.0454e-07 -0.4871   0.6262
## lin2         1.3934e-07  4.4023e-07  0.3165   0.7516
## postchange   1.0938e-04  7.3513e-05  1.4880   0.1368

## [1] "5 Regression for rate_winchib model2c on 2016-09-1/red/30-50%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -1.073e-03  -2.293e-04  -2.926e-05   2.023e-04   1.265e-03  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.318e-03  1.148e-04  37.612  < 2e-16 ***
## hour1       -3.574e-05  1.464e-04  -0.244   0.8073    
## hour2       -3.356e-04  1.464e-04  -2.293   0.0225 *  
## hour3       -7.401e-04  1.464e-04  -5.056 7.39e-07 ***
## hour4       -7.891e-04  1.464e-04  -5.390 1.41e-07 ***
## hour5       -1.067e-03  1.464e-04  -7.286 2.72e-12 ***
## hour6       -1.475e-03  1.464e-04 -10.075  < 2e-16 ***
## hour7       -1.602e-03  1.464e-04 -10.941  < 2e-16 ***
## hour8       -1.613e-03  1.464e-04 -11.018  < 2e-16 ***
## hour9       -1.891e-03  1.464e-04 -12.917  < 2e-16 ***
## hour10      -1.828e-03  1.464e-04 -12.487  < 2e-16 ***
## hour11      -1.997e-03  1.464e-04 -13.643  < 2e-16 ***
## hour12      -1.902e-03  1.464e-04 -12.991  < 2e-16 ***
## hour13      -1.745e-03  1.464e-04 -11.919  < 2e-16 ***
## hour14      -1.549e-03  1.464e-04 -10.583  < 2e-16 ***
## hour15      -1.497e-03  1.464e-04 -10.228  < 2e-16 ***
## hour16      -1.267e-03  1.464e-04  -8.653 2.86e-16 ***
## hour17      -9.958e-04  1.464e-04  -6.802 5.40e-11 ***
## hour18      -9.030e-04  1.464e-04  -6.168 2.18e-09 ***
## hour19      -1.133e-03  1.464e-04  -7.742 1.43e-13 ***
## hour20      -1.034e-03  1.464e-04  -7.061 1.11e-11 ***
## hour21      -6.774e-04  1.464e-04  -4.627 5.48e-06 ***
## hour22      -1.355e-04  1.464e-04  -0.926   0.3552    
## hour23       3.593e-05  1.464e-04   0.245   0.8063    
## dow2        -5.710e-04  8.039e-05  -7.102 8.59e-12 ***
## dow3        -8.410e-04  8.039e-05 -10.461  < 2e-16 ***
## dow4        -6.460e-04  8.039e-05  -8.035 2.02e-14 ***
## dow5        -8.429e-04  8.001e-05 -10.536  < 2e-16 ***
## dow6        -7.870e-04  8.039e-05  -9.790  < 2e-16 ***
## dow7        -1.735e-04  8.039e-05  -2.158   0.0317 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1.551134e-07)
## 
##     Null deviance: 2.163e-04  on 336  degrees of freedom
## Residual deviance: 4.762e-05  on 307  degrees of freedom
## AIC: -4296.9
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  4.3183e-03  1.1540e-04  37.4216 < 2.2e-16 ***
## hour1       -3.5740e-05  1.3255e-04  -0.2696 0.7874376    
## hour2       -3.3565e-04  1.3780e-04  -2.4357 0.0148613 *  
## hour3       -7.4013e-04  1.3517e-04  -5.4754 4.364e-08 ***
## hour4       -7.8907e-04  1.5567e-04  -5.0690 3.999e-07 ***
## hour5       -1.0666e-03  1.2958e-04  -8.2315 < 2.2e-16 ***
## hour6       -1.4750e-03  1.1474e-04 -12.8553 < 2.2e-16 ***
## hour7       -1.6018e-03  1.1314e-04 -14.1582 < 2.2e-16 ***
## hour8       -1.6130e-03  1.1585e-04 -13.9227 < 2.2e-16 ***
## hour9       -1.8911e-03  1.0901e-04 -17.3482 < 2.2e-16 ***
## hour10      -1.8281e-03  1.1797e-04 -15.4955 < 2.2e-16 ***
## hour11      -1.9973e-03  1.3338e-04 -14.9754 < 2.2e-16 ***
## hour12      -1.9019e-03  1.7101e-04 -11.1212 < 2.2e-16 ***
## hour13      -1.7450e-03  1.8354e-04  -9.5071 < 2.2e-16 ***
## hour14      -1.5494e-03  1.6432e-04  -9.4294 < 2.2e-16 ***
## hour15      -1.4973e-03  1.5104e-04  -9.9130 < 2.2e-16 ***
## hour16      -1.2668e-03  1.4879e-04  -8.5138 < 2.2e-16 ***
## hour17      -9.9580e-04  1.5163e-04  -6.5674 5.122e-11 ***
## hour18      -9.0296e-04  1.6957e-04  -5.3250 1.010e-07 ***
## hour19      -1.1335e-03  1.3481e-04  -8.4078 < 2.2e-16 ***
## hour20      -1.0337e-03  1.4497e-04  -7.1302 1.003e-12 ***
## hour21      -6.7737e-04  1.8069e-04  -3.7488 0.0001777 ***
## hour22      -1.3555e-04  1.5898e-04  -0.8526 0.3938667    
## hour23       3.5929e-05  1.4252e-04   0.2521 0.8009665    
## dow2        -5.7097e-04  1.1065e-04  -5.1603 2.466e-07 ***
## dow3        -8.4100e-04  7.6303e-05 -11.0219 < 2.2e-16 ***
## dow4        -6.4600e-04  6.5011e-05  -9.9368 < 2.2e-16 ***
## dow5        -8.4295e-04  7.0250e-05 -11.9993 < 2.2e-16 ***
## dow6        -7.8705e-04  7.3057e-05 -10.7730 < 2.2e-16 ***
## dow7        -1.7350e-04  8.3077e-05  -2.0884 0.0367614 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -1.053e-03  -2.279e-04  -1.363e-05   2.088e-04   1.177e-03  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)  
## (Intercept) -8.180e-05  5.719e-05  -1.430   0.1536  
## lin1         2.937e-09  5.870e-07   0.005   0.9960  
## lin2         2.397e-07  5.818e-07   0.412   0.6806  
## postchange   1.432e-04  8.041e-05   1.781   0.0758 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1.361603e-07)
## 
##     Null deviance: 4.7620e-05  on 336  degrees of freedom
## Residual deviance: 4.5341e-05  on 333  degrees of freedom
## AIC: -4365.4
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)  
## (Intercept) -8.1802e-05  4.5008e-05 -1.8175  0.06914 .
## lin1         2.9370e-09  4.0498e-07  0.0073  0.99421  
## lin2         2.3971e-07  5.1433e-07  0.4661  0.64116  
## postchange   1.4323e-04  6.5115e-05  2.1996  0.02783 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "6 Regression for rate_winchib model2c on 2016-09-07/red/50-100%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0010797  -0.0002476  -0.0000230   0.0001990   0.0013280  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.308e-03  1.171e-04  36.798  < 2e-16 ***
## hour1        8.087e-05  1.493e-04   0.542 0.588335    
## hour2       -3.801e-04  1.493e-04  -2.546 0.011377 *  
## hour3       -6.434e-04  1.493e-04  -4.311 2.20e-05 ***
## hour4       -6.939e-04  1.493e-04  -4.649 4.97e-06 ***
## hour5       -1.109e-03  1.493e-04  -7.427 1.11e-12 ***
## hour6       -1.514e-03  1.493e-04 -10.146  < 2e-16 ***
## hour7       -1.666e-03  1.493e-04 -11.164  < 2e-16 ***
## hour8       -1.630e-03  1.493e-04 -10.923  < 2e-16 ***
## hour9       -1.813e-03  1.493e-04 -12.149  < 2e-16 ***
## hour10      -1.778e-03  1.493e-04 -11.912  < 2e-16 ***
## hour11      -1.937e-03  1.493e-04 -12.977  < 2e-16 ***
## hour12      -1.880e-03  1.493e-04 -12.593  < 2e-16 ***
## hour13      -1.738e-03  1.493e-04 -11.641  < 2e-16 ***
## hour14      -1.456e-03  1.493e-04  -9.754  < 2e-16 ***
## hour15      -1.391e-03  1.493e-04  -9.316  < 2e-16 ***
## hour16      -1.186e-03  1.493e-04  -7.947 3.66e-14 ***
## hour17      -9.013e-04  1.493e-04  -6.038 4.49e-09 ***
## hour18      -8.921e-04  1.493e-04  -5.977 6.31e-09 ***
## hour19      -1.113e-03  1.493e-04  -7.455 9.21e-13 ***
## hour20      -9.489e-04  1.493e-04  -6.357 7.42e-10 ***
## hour21      -5.708e-04  1.493e-04  -3.824 0.000159 ***
## hour22       3.409e-06  1.493e-04   0.023 0.981793    
## hour23       2.469e-05  1.493e-04   0.165 0.868723    
## dow2        -6.447e-04  8.197e-05  -7.866 6.31e-14 ***
## dow3        -9.506e-04  8.197e-05 -11.598  < 2e-16 ***
## dow4        -6.840e-04  8.157e-05  -8.385 1.86e-15 ***
## dow5        -7.634e-04  8.197e-05  -9.313  < 2e-16 ***
## dow6        -8.403e-04  8.197e-05 -10.251  < 2e-16 ***
## dow7        -1.646e-04  8.197e-05  -2.008 0.045479 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1.612434e-07)
## 
##     Null deviance: 2.2662e-04  on 336  degrees of freedom
## Residual deviance: 4.9502e-05  on 307  degrees of freedom
## AIC: -4283.9
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  4.3076e-03  1.0979e-04  39.2349 < 2.2e-16 ***
## hour1        8.0874e-05  1.2687e-04   0.6374 0.5238331    
## hour2       -3.8006e-04  1.3434e-04  -2.8292 0.0046665 ** 
## hour3       -6.4342e-04  1.2928e-04  -4.9771 6.455e-07 ***
## hour4       -6.9388e-04  1.3480e-04  -5.1476 2.638e-07 ***
## hour5       -1.1086e-03  1.2733e-04  -8.7064 < 2.2e-16 ***
## hour6       -1.5144e-03  1.1049e-04 -13.7064 < 2.2e-16 ***
## hour7       -1.6664e-03  1.1583e-04 -14.3865 < 2.2e-16 ***
## hour8       -1.6303e-03  1.1508e-04 -14.1674 < 2.2e-16 ***
## hour9       -1.8134e-03  1.0416e-04 -17.4091 < 2.2e-16 ***
## hour10      -1.7780e-03  1.0898e-04 -16.3158 < 2.2e-16 ***
## hour11      -1.9370e-03  1.1913e-04 -16.2598 < 2.2e-16 ***
## hour12      -1.8797e-03  1.6600e-04 -11.3237 < 2.2e-16 ***
## hour13      -1.7376e-03  1.6061e-04 -10.8192 < 2.2e-16 ***
## hour14      -1.4560e-03  1.6429e-04  -8.8620 < 2.2e-16 ***
## hour15      -1.3906e-03  1.5494e-04  -8.9751 < 2.2e-16 ***
## hour16      -1.1862e-03  1.6231e-04  -7.3081 2.709e-13 ***
## hour17      -9.0132e-04  1.6017e-04  -5.6273 1.830e-08 ***
## hour18      -8.9210e-04  1.7569e-04  -5.0777 3.821e-07 ***
## hour19      -1.1128e-03  1.3934e-04  -7.9865 1.389e-15 ***
## hour20      -9.4892e-04  1.4097e-04  -6.7315 1.679e-11 ***
## hour21      -5.7079e-04  1.6863e-04  -3.3848 0.0007124 ***
## hour22       3.4091e-06  1.5989e-04   0.0213 0.9829893    
## hour23       2.4691e-05  1.5822e-04   0.1561 0.8759898    
## dow2        -6.4472e-04  1.1046e-04  -5.8369 5.318e-09 ***
## dow3        -9.5062e-04  7.0346e-05 -13.5135 < 2.2e-16 ***
## dow4        -6.8402e-04  6.1946e-05 -11.0422 < 2.2e-16 ***
## dow5        -7.6337e-04  6.7938e-05 -11.2363 < 2.2e-16 ***
## dow6        -8.4026e-04  7.4674e-05 -11.2524 < 2.2e-16 ***
## dow7        -1.6462e-04  8.4743e-05  -1.9426 0.0520686 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0011584  -0.0002256  -0.0000182   0.0002094   0.0011827  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)   
## (Intercept)  1.897e-04  5.778e-05   3.282  0.00114 **
## lin1         1.480e-06  5.931e-07   2.495  0.01307 * 
## lin2        -1.546e-06  5.878e-07  -2.630  0.00892 **
## postchange  -1.240e-04  8.123e-05  -1.526  0.12785   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1.389687e-07)
## 
##     Null deviance: 4.9502e-05  on 336  degrees of freedom
## Residual deviance: 4.6277e-05  on 333  degrees of freedom
## AIC: -4358.6
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)   
## (Intercept)  1.8965e-04  7.2560e-05  2.6137 0.008956 **
## lin1         1.4798e-06  6.1047e-07  2.4241 0.015348 * 
## lin2        -1.5462e-06  4.8246e-07 -3.2048 0.001352 **
## postchange  -1.2399e-04  8.3741e-05 -1.4807 0.138688   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "7 Regression for wincall model2c on 2016-08-17/red/0-30%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -168.176   -36.998    -3.197    42.690   147.592  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   305.68      18.77  16.290  < 2e-16 ***
## hour1         -65.56      23.93  -2.740  0.00651 ** 
## hour2        -145.06      23.93  -6.062 3.93e-09 ***
## hour3        -209.27      23.93  -8.746  < 2e-16 ***
## hour4        -284.49      23.93 -11.889  < 2e-16 ***
## hour5        -325.56      23.93 -13.606  < 2e-16 ***
## hour6        -355.20      23.93 -14.845  < 2e-16 ***
## hour7        -370.27      23.93 -15.475  < 2e-16 ***
## hour8        -381.84      23.93 -15.958  < 2e-16 ***
## hour9        -383.34      23.93 -16.021  < 2e-16 ***
## hour10       -368.91      23.93 -15.418  < 2e-16 ***
## hour11       -322.06      23.93 -13.460  < 2e-16 ***
## hour12       -200.63      23.93  -8.385 1.86e-15 ***
## hour13        -29.84      23.93  -1.247  0.21327    
## hour14         71.66      23.93   2.995  0.00297 ** 
## hour15        176.30      23.93   7.368 1.61e-12 ***
## hour16        208.94      23.93   8.732  < 2e-16 ***
## hour17        245.01      23.93  10.240  < 2e-16 ***
## hour18        290.87      23.93  12.156  < 2e-16 ***
## hour19        286.73      23.93  11.983  < 2e-16 ***
## hour20        252.23      23.93  10.541  < 2e-16 ***
## hour21        194.59      23.93   8.132 1.05e-14 ***
## hour22        130.51      23.93   5.455 1.01e-07 ***
## hour23         74.87      23.93   3.129  0.00192 ** 
## dow2          147.52      13.14  11.227  < 2e-16 ***
## dow3          172.42      13.14  13.122  < 2e-16 ***
## dow4          162.52      13.08  12.428  < 2e-16 ***
## dow5          169.69      13.14  12.914  < 2e-16 ***
## dow6          146.37      13.14  11.140  < 2e-16 ***
## dow7           62.63      13.14   4.766 2.90e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 4143.579)
## 
##     Null deviance: 22539371  on 336  degrees of freedom
## Residual deviance:  1272079  on 307  degrees of freedom
## AIC: 3793.9
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##             Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  305.680     21.170  14.4395 < 2.2e-16 ***
## hour1        -65.557     21.622  -3.0319  0.002430 ** 
## hour2       -145.057     22.146  -6.5502 5.746e-11 ***
## hour3       -209.272     23.182  -9.0273 < 2.2e-16 ***
## hour4       -284.486     23.397 -12.1591 < 2.2e-16 ***
## hour5       -325.557     22.917 -14.2060 < 2.2e-16 ***
## hour6       -355.200     23.803 -14.9223 < 2.2e-16 ***
## hour7       -370.272     24.963 -14.8329 < 2.2e-16 ***
## hour8       -381.843     23.764 -16.0679 < 2.2e-16 ***
## hour9       -383.343     23.581 -16.2565 < 2.2e-16 ***
## hour10      -368.915     23.229 -15.8815 < 2.2e-16 ***
## hour11      -322.057     21.688 -14.8494 < 2.2e-16 ***
## hour12      -200.629     18.764 -10.6921 < 2.2e-16 ***
## hour13       -29.843     21.715  -1.3743  0.169341    
## hour14        71.657     23.872   3.0017  0.002685 ** 
## hour15       176.300     27.581   6.3922 1.635e-10 ***
## hour16       208.943     27.527   7.5903 3.190e-14 ***
## hour17       245.014     24.410  10.0376 < 2.2e-16 ***
## hour18       290.871     30.209   9.6286 < 2.2e-16 ***
## hour19       286.728     26.614  10.7734 < 2.2e-16 ***
## hour20       252.228     29.297   8.6094 < 2.2e-16 ***
## hour21       194.585     23.955   8.1231 4.544e-16 ***
## hour22       130.514     22.737   5.7401 9.460e-09 ***
## hour23        74.871     19.642   3.8118  0.000138 ***
## dow2         147.521     18.355   8.0369 9.212e-16 ***
## dow3         172.417     16.333  10.5562 < 2.2e-16 ***
## dow4         162.517     16.012  10.1497 < 2.2e-16 ***
## dow5         169.688     16.145  10.5100 < 2.2e-16 ***
## dow6         146.375     15.649   9.3535 < 2.2e-16 ***
## dow7          62.625     19.462   3.2178  0.001292 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -160.035   -37.021     1.096    44.614   155.968  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)
## (Intercept)   9.66297    9.49595   1.018    0.310
## lin1          0.01872    0.09747   0.192    0.848
## lin2         -0.01804    0.09660  -0.187    0.852
## postchange  -16.18079   13.35030  -1.212    0.226
## 
## (Dispersion parameter for gaussian family taken to be 3753.552)
## 
##     Null deviance: 1272079  on 336  degrees of freedom
## Residual deviance: 1249933  on 333  degrees of freedom
## AIC: 3736
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error z value Pr(>|z|)
## (Intercept)   9.662965   9.289869  1.0402   0.2983
## lin1          0.018722   0.076152  0.2459   0.8058
## lin2         -0.018039   0.072759 -0.2479   0.8042
## postchange  -16.180786  10.835204 -1.4934   0.1353

## [1] "8 Regression for wincall model2c on 2016-09-1/red/30-50%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -179.595   -45.163    -3.702    52.298   183.226  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   314.17      20.76  15.134  < 2e-16 ***
## hour1         -80.36      26.47  -3.036 0.002605 ** 
## hour2        -152.93      26.47  -5.777 1.86e-08 ***
## hour3        -220.64      26.47  -8.336 2.62e-15 ***
## hour4        -283.14      26.47 -10.697  < 2e-16 ***
## hour5        -323.36      26.47 -12.216  < 2e-16 ***
## hour6        -349.86      26.47 -13.217  < 2e-16 ***
## hour7        -367.71      26.47 -13.892  < 2e-16 ***
## hour8        -379.14      26.47 -14.324  < 2e-16 ***
## hour9        -374.78      26.47 -14.159  < 2e-16 ***
## hour10       -368.43      26.47 -13.919  < 2e-16 ***
## hour11       -328.43      26.47 -12.408  < 2e-16 ***
## hour12       -224.43      26.47  -8.479 9.72e-16 ***
## hour13        -66.71      26.47  -2.520 0.012231 *  
## hour14         32.64      26.47   1.233 0.218421    
## hour15         97.93      26.47   3.700 0.000256 ***
## hour16        143.57      26.47   5.424 1.18e-07 ***
## hour17        177.14      26.47   6.692 1.04e-10 ***
## hour18        233.57      26.47   8.824  < 2e-16 ***
## hour19        213.43      26.47   8.063 1.68e-14 ***
## hour20        205.64      26.47   7.769 1.20e-13 ***
## hour21        172.43      26.47   6.514 2.99e-10 ***
## hour22         90.36      26.47   3.414 0.000727 ***
## hour23         46.22      26.47   1.746 0.081815 .  
## dow2           95.25      14.54   6.553 2.38e-10 ***
## dow3          150.44      14.54  10.350  < 2e-16 ***
## dow4          149.35      14.54  10.275  < 2e-16 ***
## dow5          150.85      14.47  10.428  < 2e-16 ***
## dow6          137.42      14.54   9.454  < 2e-16 ***
## dow7           48.52      14.54   3.338 0.000947 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 5070.812)
## 
##     Null deviance: 18896875  on 336  degrees of freedom
## Residual deviance:  1556739  on 307  degrees of freedom
## AIC: 3862
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##             Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  314.165     20.292  15.4819 < 2.2e-16 ***
## hour1        -80.356     20.859  -3.8523 0.0001170 ***
## hour2       -152.927     21.969  -6.9611 3.376e-12 ***
## hour3       -220.641     22.215  -9.9321 < 2.2e-16 ***
## hour4       -283.141     22.058 -12.8363 < 2.2e-16 ***
## hour5       -323.356     23.462 -13.7824 < 2.2e-16 ***
## hour6       -349.856     23.494 -14.8914 < 2.2e-16 ***
## hour7       -367.713     22.976 -16.0043 < 2.2e-16 ***
## hour8       -379.141     22.592 -16.7820 < 2.2e-16 ***
## hour9       -374.784     22.239 -16.8524 < 2.2e-16 ***
## hour10      -368.427     22.027 -16.7259 < 2.2e-16 ***
## hour11      -328.427     20.916 -15.7021 < 2.2e-16 ***
## hour12      -224.427     20.059 -11.1881 < 2.2e-16 ***
## hour13       -66.713     21.529  -3.0987 0.0019436 ** 
## hour14        32.644     25.723   1.2691 0.2044198    
## hour15        97.930     28.180   3.4751 0.0005106 ***
## hour16       143.573     29.023   4.9468 7.545e-07 ***
## hour17       177.144     28.474   6.2213 4.931e-10 ***
## hour18       233.573     29.817   7.8335 4.745e-15 ***
## hour19       213.430     33.258   6.4174 1.386e-10 ***
## hour20       205.644     33.929   6.0610 1.353e-09 ***
## hour21       172.430     30.103   5.7281 1.016e-08 ***
## hour22        90.359     28.383   3.1835 0.0014551 ** 
## hour23        46.216     24.265   1.9046 0.0568340 .  
## dow2          95.250     19.914   4.7831 1.726e-06 ***
## dow3         150.437     17.306   8.6929 < 2.2e-16 ***
## dow4         149.354     15.053   9.9221 < 2.2e-16 ***
## dow5         150.854     14.784  10.2040 < 2.2e-16 ***
## dow6         137.417     15.219   9.0294 < 2.2e-16 ***
## dow7          48.521     17.127   2.8331 0.0046103 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -161.724   -49.937    -1.773    44.686   162.767  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)  
## (Intercept)  25.63804   10.22378   2.508   0.0126 *
## lin1          0.09905    0.10494   0.944   0.3459  
## lin2         -0.10081    0.10401  -0.969   0.3331  
## postchange  -34.33582   14.37356  -2.389   0.0175 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 4350.998)
## 
##     Null deviance: 1556739  on 336  degrees of freedom
## Residual deviance: 1448882  on 333  degrees of freedom
## AIC: 3785.8
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error z value Pr(>|z|)   
## (Intercept)  25.638038   9.716927  2.6385 0.008328 **
## lin1          0.099055   0.085114  1.1638 0.244508   
## lin2         -0.100809   0.084187 -1.1974 0.231133   
## postchange  -34.335817  12.184789 -2.8179 0.004834 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "9 Regression for wincall model2c on 2016-09-07/red/50-100%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -341.71   -48.63    -5.27    51.86   169.87  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   304.93      21.83  13.971  < 2e-16 ***
## hour1         -81.56      27.83  -2.931 0.003635 ** 
## hour2        -149.56      27.83  -5.374 1.52e-07 ***
## hour3        -207.63      27.83  -7.461 8.89e-13 ***
## hour4        -263.35      27.83  -9.463  < 2e-16 ***
## hour5        -312.92      27.83 -11.244  < 2e-16 ***
## hour6        -342.13      27.83 -12.294  < 2e-16 ***
## hour7        -355.28      27.83 -12.766  < 2e-16 ***
## hour8        -370.63      27.83 -13.318  < 2e-16 ***
## hour9        -363.99      27.83 -13.079  < 2e-16 ***
## hour10       -354.78      27.83 -12.748  < 2e-16 ***
## hour11       -314.35      27.83 -11.296  < 2e-16 ***
## hour12       -213.35      27.83  -7.666 2.35e-13 ***
## hour13        -52.78      27.83  -1.896 0.058842 .  
## hour14         33.94      27.83   1.220 0.223589    
## hour15         89.30      27.83   3.209 0.001474 ** 
## hour16        146.65      27.83   5.270 2.58e-07 ***
## hour17        173.94      27.83   6.250 1.37e-09 ***
## hour18        225.72      27.83   8.111 1.22e-14 ***
## hour19        214.87      27.83   7.721 1.64e-13 ***
## hour20        183.44      27.83   6.592 1.90e-10 ***
## hour21        144.58      27.83   5.195 3.73e-07 ***
## hour22         65.01      27.83   2.336 0.020135 *  
## hour23         48.15      27.83   1.730 0.084587 .  
## dow2           91.33      15.28   5.976 6.32e-09 ***
## dow3          153.25      15.28  10.028  < 2e-16 ***
## dow4          157.21      15.21  10.336  < 2e-16 ***
## dow5          132.21      15.28   8.651 2.90e-16 ***
## dow6          129.58      15.28   8.479 9.67e-16 ***
## dow7           51.38      15.28   3.362 0.000872 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 5605.065)
## 
##     Null deviance: 17678273  on 336  degrees of freedom
## Residual deviance:  1720755  on 307  degrees of freedom
## AIC: 3895.7
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##             Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  304.926     19.998  15.2476 < 2.2e-16 ***
## hour1        -81.562     20.569  -3.9654 7.328e-05 ***
## hour2       -149.562     21.550  -6.9402 3.916e-12 ***
## hour3       -207.633     21.943  -9.4625 < 2.2e-16 ***
## hour4       -263.348     22.003 -11.9688 < 2.2e-16 ***
## hour5       -312.919     22.835 -13.7034 < 2.2e-16 ***
## hour6       -342.133     22.813 -14.9975 < 2.2e-16 ***
## hour7       -355.276     21.406 -16.5967 < 2.2e-16 ***
## hour8       -370.633     22.297 -16.6228 < 2.2e-16 ***
## hour9       -363.991     22.388 -16.2586 < 2.2e-16 ***
## hour10      -354.776     21.370 -16.6013 < 2.2e-16 ***
## hour11      -314.348     20.548 -15.2981 < 2.2e-16 ***
## hour12      -213.348     19.977 -10.6796 < 2.2e-16 ***
## hour13       -52.776     22.686  -2.3264 0.0199973 *  
## hour14        33.938     26.213   1.2947 0.1954183    
## hour15        89.295     31.239   2.8585 0.0042568 ** 
## hour16       146.652     26.362   5.5631 2.650e-08 ***
## hour17       173.938     25.637   6.7845 1.164e-11 ***
## hour18       225.724     29.779   7.5800 3.457e-14 ***
## hour19       214.867     31.748   6.7679 1.307e-11 ***
## hour20       183.438     33.150   5.5335 3.139e-08 ***
## hour21       144.581     38.568   3.7487 0.0001778 ***
## hour22        65.009     35.844   1.8137 0.0697293 .  
## hour23        48.152     23.483   2.0505 0.0403165 *  
## dow2          91.333     18.408   4.9616 6.991e-07 ***
## dow3         153.250     17.315   8.8506 < 2.2e-16 ***
## dow4         157.206     15.307  10.2700 < 2.2e-16 ***
## dow5         132.208     18.414   7.1797 6.986e-13 ***
## dow6         129.583     14.758   8.7805 < 2.2e-16 ***
## dow7          51.375     16.330   3.1460 0.0016552 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -336.69   -45.48    -6.31    48.90   192.97  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)   
## (Intercept) -30.7293    10.8747  -2.826  0.00500 **
## lin1         -0.2868     0.1116  -2.569  0.01063 * 
## lin2          0.2945     0.1106   2.662  0.00814 **
## postchange   12.4490    15.2886   0.814  0.41608   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 4922.627)
## 
##     Null deviance: 1720755  on 336  degrees of freedom
## Residual deviance: 1639235  on 333  degrees of freedom
## AIC: 3827.4
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error z value Pr(>|z|)   
## (Intercept) -30.72925   12.43926 -2.4703 0.013498 * 
## lin1         -0.28678    0.11082 -2.5877 0.009661 **
## lin2          0.29451    0.10739  2.7425 0.006098 **
## postchange   12.44899   16.21733  0.7676 0.442704   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "10 Regression for rate_wincall model2c on 2016-08-17/red/0-30%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -1.188e-03  -1.775e-04  -2.145e-05   1.489e-04   2.139e-03  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.262e-03  1.017e-04  41.908  < 2e-16 ***
## hour1       -6.334e-04  1.297e-04  -4.885 1.67e-06 ***
## hour2       -1.253e-03  1.297e-04  -9.661  < 2e-16 ***
## hour3       -1.789e-03  1.297e-04 -13.797  < 2e-16 ***
## hour4       -2.525e-03  1.297e-04 -19.471  < 2e-16 ***
## hour5       -2.970e-03  1.297e-04 -22.901  < 2e-16 ***
## hour6       -3.369e-03  1.297e-04 -25.977  < 2e-16 ***
## hour7       -3.551e-03  1.297e-04 -27.385  < 2e-16 ***
## hour8       -3.697e-03  1.297e-04 -28.508  < 2e-16 ***
## hour9       -3.729e-03  1.297e-04 -28.755  < 2e-16 ***
## hour10      -3.622e-03  1.297e-04 -27.935  < 2e-16 ***
## hour11      -3.401e-03  1.297e-04 -26.230  < 2e-16 ***
## hour12      -2.907e-03  1.297e-04 -22.416  < 2e-16 ***
## hour13      -2.182e-03  1.297e-04 -16.823  < 2e-16 ***
## hour14      -1.689e-03  1.297e-04 -13.026  < 2e-16 ***
## hour15      -1.097e-03  1.297e-04  -8.457 1.13e-15 ***
## hour16      -7.187e-04  1.297e-04  -5.542 6.43e-08 ***
## hour17      -3.760e-04  1.297e-04  -2.899 0.004009 ** 
## hour18      -8.846e-05  1.297e-04  -0.682 0.495642    
## hour19       1.088e-04  1.297e-04   0.839 0.401980    
## hour20       2.004e-04  1.297e-04   1.545 0.123324    
## hour21       4.757e-04  1.297e-04   3.668 0.000288 ***
## hour22       7.147e-04  1.297e-04   5.511 7.55e-08 ***
## hour23       6.710e-04  1.297e-04   5.174 4.14e-07 ***
## dow2        -4.650e-04  7.121e-05  -6.530 2.73e-10 ***
## dow3        -4.073e-06  7.121e-05  -0.057 0.954420    
## dow4         2.537e-05  7.087e-05   0.358 0.720631    
## dow5        -1.845e-04  7.121e-05  -2.591 0.010039 *  
## dow6         3.065e-05  7.121e-05   0.430 0.667160    
## dow7         6.396e-04  7.121e-05   8.983  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1.216986e-07)
## 
##     Null deviance: 8.7506e-04  on 336  degrees of freedom
## Residual deviance: 3.7361e-05  on 307  degrees of freedom
## AIC: -4378.7
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  4.2620e-03  1.3422e-04  31.7535 < 2.2e-16 ***
## hour1       -6.3343e-04  1.4973e-04  -4.2305 2.332e-05 ***
## hour2       -1.2528e-03  1.4770e-04  -8.4826 < 2.2e-16 ***
## hour3       -1.7892e-03  1.4364e-04 -12.4563 < 2.2e-16 ***
## hour4       -2.5249e-03  1.4504e-04 -17.4089 < 2.2e-16 ***
## hour5       -2.9697e-03  1.3371e-04 -22.2104 < 2.2e-16 ***
## hour6       -3.3685e-03  1.3660e-04 -24.6595 < 2.2e-16 ***
## hour7       -3.5512e-03  1.4880e-04 -23.8661 < 2.2e-16 ***
## hour8       -3.6967e-03  1.3290e-04 -27.8169 < 2.2e-16 ***
## hour9       -3.7288e-03  1.3657e-04 -27.3026 < 2.2e-16 ***
## hour10      -3.6224e-03  1.3660e-04 -26.5180 < 2.2e-16 ***
## hour11      -3.4014e-03  1.3081e-04 -26.0020 < 2.2e-16 ***
## hour12      -2.9067e-03  1.3556e-04 -21.4429 < 2.2e-16 ***
## hour13      -2.1815e-03  1.6101e-04 -13.5495 < 2.2e-16 ***
## hour14      -1.6891e-03  1.7938e-04  -9.4165 < 2.2e-16 ***
## hour15      -1.0966e-03  2.1288e-04  -5.1512 2.588e-07 ***
## hour16      -7.1871e-04  1.4596e-04  -4.9239 8.483e-07 ***
## hour17      -3.7598e-04  1.4550e-04  -2.5840 0.0097651 ** 
## hour18      -8.8461e-05  1.3871e-04  -0.6377 0.5236471    
## hour19       1.0883e-04  1.4725e-04   0.7391 0.4598419    
## hour20       2.0038e-04  1.3422e-04   1.4929 0.1354718    
## hour21       4.7569e-04  1.4178e-04   3.3551 0.0007932 ***
## hour22       7.1469e-04  1.7215e-04   4.1515 3.303e-05 ***
## hour23       6.7095e-04  1.9872e-04   3.3763 0.0007347 ***
## dow2        -4.6497e-04  7.7105e-05  -6.0303 1.637e-09 ***
## dow3        -4.0735e-06  6.7128e-05  -0.0607 0.9516120    
## dow4         2.5367e-05  6.2220e-05   0.4077 0.6834970    
## dow5        -1.8447e-04  7.2105e-05  -2.5584 0.0105160 *  
## dow6         3.0653e-05  6.5951e-05   0.4648 0.6420848    
## dow7         6.3965e-04  1.0256e-04   6.2367 4.469e-10 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -1.123e-03  -1.883e-04  -1.144e-05   1.543e-04   2.096e-03  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)  
## (Intercept)  1.028e-04  5.124e-05   2.007   0.0456 *
## lin1         7.382e-07  5.259e-07   1.403   0.1614  
## lin2        -7.078e-07  5.213e-07  -1.358   0.1755  
## postchange  -8.359e-05  7.204e-05  -1.160   0.2467  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1.092985e-07)
## 
##     Null deviance: 3.7361e-05  on 336  degrees of freedom
## Residual deviance: 3.6396e-05  on 333  degrees of freedom
## AIC: -4439.5
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)  
## (Intercept)  1.0283e-04  4.8386e-05  2.1252  0.03357 *
## lin1         7.3816e-07  3.7100e-07  1.9897  0.04663 *
## lin2        -7.0777e-07  4.0276e-07 -1.7573  0.07887 .
## postchange  -8.3595e-05  6.1669e-05 -1.3555  0.17524  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "11 Regression for rate_wincall model2c on 2016-09-1/red/30-50%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0017696  -0.0002415   0.0000117   0.0002054   0.0075458  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.395e-03  1.800e-04  24.409  < 2e-16 ***
## hour1       -8.711e-04  2.296e-04  -3.795 0.000178 ***
## hour2       -1.456e-03  2.296e-04  -6.342 8.11e-10 ***
## hour3       -2.048e-03  2.296e-04  -8.923  < 2e-16 ***
## hour4       -2.666e-03  2.296e-04 -11.613  < 2e-16 ***
## hour5       -3.106e-03  2.296e-04 -13.530  < 2e-16 ***
## hour6       -3.489e-03  2.296e-04 -15.196  < 2e-16 ***
## hour7       -3.724e-03  2.296e-04 -16.220  < 2e-16 ***
## hour8       -3.868e-03  2.296e-04 -16.849  < 2e-16 ***
## hour9       -3.829e-03  2.296e-04 -16.680  < 2e-16 ***
## hour10      -3.801e-03  2.296e-04 -16.558  < 2e-16 ***
## hour11      -3.585e-03  2.296e-04 -15.615  < 2e-16 ***
## hour12      -3.170e-03  2.296e-04 -13.809  < 2e-16 ***
## hour13      -2.462e-03  2.296e-04 -10.724  < 2e-16 ***
## hour14      -1.979e-03  2.296e-04  -8.621 3.60e-16 ***
## hour15      -1.607e-03  2.296e-04  -7.000 1.62e-11 ***
## hour16      -1.102e-03  2.296e-04  -4.800 2.48e-06 ***
## hour17      -8.035e-04  2.296e-04  -3.500 0.000534 ***
## hour18      -4.019e-04  2.296e-04  -1.751 0.081003 .  
## hour19      -3.341e-04  2.296e-04  -1.455 0.146613    
## hour20      -4.515e-05  2.296e-04  -0.197 0.844228    
## hour21       3.746e-04  2.296e-04   1.632 0.103764    
## hour22       6.739e-04  2.296e-04   2.935 0.003583 ** 
## hour23       8.981e-04  2.296e-04   3.912 0.000113 ***
## dow2        -3.373e-04  1.261e-04  -2.676 0.007857 ** 
## dow3        -2.612e-05  1.261e-04  -0.207 0.835984    
## dow4        -8.034e-05  1.261e-04  -0.637 0.524430    
## dow5        -3.059e-04  1.255e-04  -2.438 0.015327 *  
## dow6        -6.574e-05  1.261e-04  -0.522 0.602385    
## dow7         5.206e-04  1.261e-04   4.130 4.69e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 3.814184e-07)
## 
##     Null deviance: 0.0009546  on 336  degrees of freedom
## Residual deviance: 0.0001171  on 307  degrees of freedom
## AIC: -3993.7
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  4.3946e-03  1.3916e-04  31.5788 < 2.2e-16 ***
## hour1       -8.7114e-04  1.5998e-04  -5.4452 5.173e-08 ***
## hour2       -1.4559e-03  1.6173e-04  -9.0019 < 2.2e-16 ***
## hour3       -2.0484e-03  1.4657e-04 -13.9759 < 2.2e-16 ***
## hour4       -2.6660e-03  1.4506e-04 -18.3786 < 2.2e-16 ***
## hour5       -3.1061e-03  1.3821e-04 -22.4728 < 2.2e-16 ***
## hour6       -3.4886e-03  1.4244e-04 -24.4906 < 2.2e-16 ***
## hour7       -3.7235e-03  1.4100e-04 -26.4074 < 2.2e-16 ***
## hour8       -3.8679e-03  1.4245e-04 -27.1526 < 2.2e-16 ***
## hour9       -3.8292e-03  1.4960e-04 -25.5955 < 2.2e-16 ***
## hour10      -3.8012e-03  1.4436e-04 -26.3319 < 2.2e-16 ***
## hour11      -3.5847e-03  1.3787e-04 -26.0002 < 2.2e-16 ***
## hour12      -3.1700e-03  1.4703e-04 -21.5603 < 2.2e-16 ***
## hour13      -2.4618e-03  1.5892e-04 -15.4913 < 2.2e-16 ***
## hour14      -1.9790e-03  1.6969e-04 -11.6627 < 2.2e-16 ***
## hour15      -1.6069e-03  1.8139e-04  -8.8586 < 2.2e-16 ***
## hour16      -1.1019e-03  1.7608e-04  -6.2581 3.897e-10 ***
## hour17      -8.0347e-04  1.8197e-04  -4.4155 1.008e-05 ***
## hour18      -4.0190e-04  1.7479e-04  -2.2993 0.0214899 *  
## hour19      -3.3409e-04  1.6845e-04  -1.9833 0.0473343 *  
## hour20      -4.5146e-05  1.4942e-04  -0.3021 0.7625418    
## hour21       3.7459e-04  1.5023e-04   2.4935 0.0126485 *  
## hour22       6.7387e-04  3.4731e-04   1.9403 0.0523482 .  
## hour23       8.9814e-04  6.1574e-04   1.4586 0.1446665    
## dow2        -3.3731e-04  8.9061e-05  -3.7875 0.0001522 ***
## dow3        -2.6122e-05  2.0556e-04  -0.1271 0.8988777    
## dow4        -8.0336e-05  7.0668e-05  -1.1368 0.2556155    
## dow5        -3.0590e-04  8.0895e-05  -3.7814 0.0001559 ***
## dow6        -6.5745e-05  7.0158e-05  -0.9371 0.3487110    
## dow7         5.2060e-04  9.7234e-05   5.3541 8.598e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0017534  -0.0002340   0.0000151   0.0002088   0.0075375  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.650e-05  9.177e-05  -0.180    0.857
## lin1        -4.553e-07  9.420e-07  -0.483    0.629
## lin2         5.106e-07  9.336e-07   0.547    0.585
## postchange  -4.822e-05  1.290e-04  -0.374    0.709
## 
## (Dispersion parameter for gaussian family taken to be 3.505918e-07)
## 
##     Null deviance: 0.00011710  on 336  degrees of freedom
## Residual deviance: 0.00011675  on 333  degrees of freedom
## AIC: -4046.7
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)
## (Intercept) -1.6504e-05  4.4991e-05 -0.3668   0.7137
## lin1        -4.5527e-07  4.5805e-07 -0.9939   0.3203
## lin2         5.1061e-07  1.3365e-06  0.3820   0.7024
## postchange  -4.8223e-05  8.4465e-05 -0.5709   0.5681

## [1] "12 Regression for rate_wincall model2c on 2016-09-07/red/50-100%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0029927  -0.0002288   0.0000111   0.0002253   0.0076133  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.365e-03  1.935e-04  22.565  < 2e-16 ***
## hour1       -8.627e-04  2.467e-04  -3.497 0.000539 ***
## hour2       -1.411e-03  2.467e-04  -5.720 2.52e-08 ***
## hour3       -1.899e-03  2.467e-04  -7.698 1.91e-13 ***
## hour4       -2.431e-03  2.467e-04  -9.855  < 2e-16 ***
## hour5       -3.007e-03  2.467e-04 -12.191  < 2e-16 ***
## hour6       -3.434e-03  2.467e-04 -13.923  < 2e-16 ***
## hour7       -3.631e-03  2.467e-04 -14.719  < 2e-16 ***
## hour8       -3.809e-03  2.467e-04 -15.443  < 2e-16 ***
## hour9       -3.735e-03  2.467e-04 -15.143  < 2e-16 ***
## hour10      -3.679e-03  2.467e-04 -14.915  < 2e-16 ***
## hour11      -3.454e-03  2.467e-04 -14.004  < 2e-16 ***
## hour12      -3.053e-03  2.467e-04 -12.377  < 2e-16 ***
## hour13      -2.332e-03  2.467e-04  -9.453  < 2e-16 ***
## hour14      -1.934e-03  2.467e-04  -7.842 7.40e-14 ***
## hour15      -1.602e-03  2.467e-04  -6.494 3.35e-10 ***
## hour16      -1.061e-03  2.467e-04  -4.300 2.30e-05 ***
## hour17      -7.744e-04  2.467e-04  -3.139 0.001858 ** 
## hour18      -4.266e-04  2.467e-04  -1.730 0.084711 .  
## hour19      -2.684e-04  2.467e-04  -1.088 0.277420    
## hour20      -1.256e-04  2.467e-04  -0.509 0.610972    
## hour21       2.075e-04  2.467e-04   0.841 0.400822    
## hour22       4.566e-04  2.467e-04   1.851 0.065128 .  
## hour23       9.212e-04  2.467e-04   3.734 0.000224 ***
## dow2        -3.856e-04  1.355e-04  -2.847 0.004717 ** 
## dow3        -8.745e-05  1.355e-04  -0.646 0.519037    
## dow4        -9.073e-05  1.348e-04  -0.673 0.501433    
## dow5        -3.807e-04  1.355e-04  -2.811 0.005263 ** 
## dow6        -2.000e-04  1.355e-04  -1.477 0.140830    
## dow7         4.940e-04  1.355e-04   3.647 0.000312 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 4.403789e-07)
## 
##     Null deviance: 0.00090768  on 336  degrees of freedom
## Residual deviance: 0.00013520  on 307  degrees of freedom
## AIC: -3945.3
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  4.3653e-03  1.3535e-04  32.2531 < 2.2e-16 ***
## hour1       -8.6270e-04  1.5421e-04  -5.5942 2.217e-08 ***
## hour2       -1.4111e-03  1.4098e-04 -10.0087 < 2.2e-16 ***
## hour3       -1.8989e-03  1.4155e-04 -13.4153 < 2.2e-16 ***
## hour4       -2.4310e-03  1.4502e-04 -16.7635 < 2.2e-16 ***
## hour5       -3.0071e-03  1.3273e-04 -22.6564 < 2.2e-16 ***
## hour6       -3.4345e-03  1.3157e-04 -26.1043 < 2.2e-16 ***
## hour7       -3.6307e-03  1.3046e-04 -27.8310 < 2.2e-16 ***
## hour8       -3.8095e-03  1.3317e-04 -28.6062 < 2.2e-16 ***
## hour9       -3.7353e-03  1.3796e-04 -27.0756 < 2.2e-16 ***
## hour10      -3.6792e-03  1.3181e-04 -27.9122 < 2.2e-16 ***
## hour11      -3.4543e-03  1.3458e-04 -25.6674 < 2.2e-16 ***
## hour12      -3.0531e-03  1.4441e-04 -21.1421 < 2.2e-16 ***
## hour13      -2.3318e-03  1.5523e-04 -15.0218 < 2.2e-16 ***
## hour14      -1.9344e-03  1.6003e-04 -12.0871 < 2.2e-16 ***
## hour15      -1.6020e-03  1.8753e-04  -8.5429 < 2.2e-16 ***
## hour16      -1.0606e-03  1.4635e-04  -7.2474 4.249e-13 ***
## hour17      -7.7438e-04  1.7810e-04  -4.3479 1.374e-05 ***
## hour18      -4.2664e-04  1.5429e-04  -2.7651  0.005690 ** 
## hour19      -2.6840e-04  1.6865e-04  -1.5915  0.111502    
## hour20      -1.2561e-04  1.6805e-04  -0.7474  0.454797    
## hour21       2.0753e-04  2.3947e-04   0.8666  0.386143    
## hour22       4.5660e-04  4.1117e-04   1.1105  0.266791    
## hour23       9.2115e-04  6.1812e-04   1.4902  0.136163    
## dow2        -3.8559e-04  8.9526e-05  -4.3070 1.655e-05 ***
## dow3        -8.7449e-05  2.0778e-04  -0.4209  0.673847    
## dow4        -9.0731e-05  7.2956e-05  -1.2436  0.213630    
## dow5        -3.8072e-04  1.1898e-04  -3.1997  0.001376 ** 
## dow6        -2.0001e-04  8.1724e-05  -2.4474  0.014391 *  
## dow7         4.9402e-04  9.2947e-05   5.3151 1.066e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0029702  -0.0002166   0.0000196   0.0002217   0.0075429  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)
## (Intercept)  7.080e-05  9.857e-05   0.718    0.473
## lin1         4.241e-07  1.012e-06   0.419    0.675
## lin2        -3.230e-07  1.003e-06  -0.322    0.748
## postchange  -7.842e-05  1.386e-04  -0.566    0.572
## 
## (Dispersion parameter for gaussian family taken to be 4.044257e-07)
## 
##     Null deviance: 0.00013520  on 336  degrees of freedom
## Residual deviance: 0.00013467  on 333  degrees of freedom
## AIC: -3998.6
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)
## (Intercept)  7.0798e-05  2.1841e-04  0.3242   0.7458
## lin1         4.2406e-07  1.9250e-06  0.2203   0.8256
## lin2        -3.2304e-07  5.6371e-07 -0.5731   0.5666
## postchange  -7.8421e-05  2.2949e-04 -0.3417   0.7326

## [1] "13 Regression for xbxchib model2c on 2016-10-11/purple/0-10%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -118.303   -17.402    -3.291    17.663   184.862  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  396.32342   10.99981  36.030  < 2e-16 ***
## hour1         -7.19396   13.96065  -0.515 0.606725    
## hour2        -67.26539   13.96065  -4.818 2.31e-06 ***
## hour3       -115.83682   13.96065  -8.297 3.77e-15 ***
## hour4       -156.26539   13.96065 -11.193  < 2e-16 ***
## hour5       -204.69396   13.96065 -14.662  < 2e-16 ***
## hour6       -250.90825   13.96065 -17.973  < 2e-16 ***
## hour7       -291.97967   13.96065 -20.914  < 2e-16 ***
## hour8       -318.94079   14.23979 -22.398  < 2e-16 ***
## hour9       -316.24849   14.23979 -22.209  < 2e-16 ***
## hour10      -318.86387   14.23979 -22.392  < 2e-16 ***
## hour11      -309.78695   14.23979 -21.755  < 2e-16 ***
## hour12      -260.19396   13.96065 -18.638  < 2e-16 ***
## hour13      -222.62253   13.96065 -15.946  < 2e-16 ***
## hour14      -204.15789   14.55690 -14.025  < 2e-16 ***
## hour15      -157.90789   14.55690 -10.848  < 2e-16 ***
## hour16      -130.76539   13.96065  -9.367  < 2e-16 ***
## hour17      -117.76539   13.96065  -8.436 1.46e-15 ***
## hour18       -89.19396   13.96065  -6.389 6.44e-10 ***
## hour19       -93.05110   13.96065  -6.665 1.28e-10 ***
## hour20       -77.09464   14.23979  -5.414 1.27e-07 ***
## hour21       -28.63310   14.23979  -2.011 0.045252 *  
## hour22         0.09175   13.96065   0.007 0.994760    
## hour23       -28.69396   13.96065  -2.055 0.040720 *  
## dow2         -21.84544    7.75574  -2.817 0.005177 ** 
## dow3         -25.03887    7.71892  -3.244 0.001314 ** 
## dow4         -17.49127    7.75574  -2.255 0.024844 *  
## dow5         -16.86627    7.75574  -2.175 0.030443 *  
## dow6         -29.83829    8.04295  -3.710 0.000247 ***
## dow7          16.17391    7.83121   2.065 0.039761 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1410.542)
## 
##     Null deviance: 4282918  on 326  degrees of freedom
## Residual deviance:  418931  on 297  degrees of freedom
## AIC: 3329.8
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  396.323422   12.745091  31.0962 < 2.2e-16 ***
## hour1         -7.193960   20.469219  -0.3515 0.7252488    
## hour2        -67.265389   14.818457  -4.5393 5.644e-06 ***
## hour3       -115.836818   18.008942  -6.4322 1.258e-10 ***
## hour4       -156.265389   16.258455  -9.6113 < 2.2e-16 ***
## hour5       -204.693960   14.675892 -13.9476 < 2.2e-16 ***
## hour6       -250.908246   13.146463 -19.0856 < 2.2e-16 ***
## hour7       -291.979675   12.939889 -22.5643 < 2.2e-16 ***
## hour8       -318.940794   13.571912 -23.5001 < 2.2e-16 ***
## hour9       -316.248486   13.108979 -24.1246 < 2.2e-16 ***
## hour10      -318.863871   13.162141 -24.2258 < 2.2e-16 ***
## hour11      -309.786948   14.370308 -21.5574 < 2.2e-16 ***
## hour12      -260.193960   13.948114 -18.6544 < 2.2e-16 ***
## hour13      -222.622532   13.857136 -16.0656 < 2.2e-16 ***
## hour14      -204.157891   13.824301 -14.7680 < 2.2e-16 ***
## hour15      -157.907891   17.464036  -9.0419 < 2.2e-16 ***
## hour16      -130.765389   15.007156  -8.7135 < 2.2e-16 ***
## hour17      -117.765389   16.594277  -7.0967 1.277e-12 ***
## hour18       -89.193960   17.106225  -5.2141 1.847e-07 ***
## hour19       -93.051103   15.985055  -5.8211 5.845e-09 ***
## hour20       -77.094640   18.723627  -4.1175 3.830e-05 ***
## hour21       -28.633101   14.910035  -1.9204 0.0548085 .  
## hour22         0.091754   17.596833   0.0052 0.9958397    
## hour23       -28.693960   21.870389  -1.3120 0.1895200    
## dow2         -21.845439    5.440084  -4.0156 5.928e-05 ***
## dow3         -25.038867    7.165975  -3.4941 0.0004756 ***
## dow4         -17.491272    8.129569  -2.1516 0.0314319 *  
## dow5         -16.866272    7.930487  -2.1268 0.0334397 *  
## dow6         -29.838295    6.916071  -4.3143 1.601e-05 ***
## dow7          16.173913    7.756691   2.0852 0.0370551 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -125.384   -19.608    -0.806    16.318   168.037  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  -7.75812    5.30086  -1.464 0.144287    
## lin1         -0.22147    0.05784  -3.829 0.000154 ***
## lin2          0.21183    0.05228   4.052 6.38e-05 ***
## postchange  -19.24286    7.34105  -2.621 0.009175 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1099.41)
## 
##     Null deviance: 418931  on 326  degrees of freedom
## Residual deviance: 355110  on 323  degrees of freedom
## AIC: 3223.8
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error z value  Pr(>|z|)    
## (Intercept)  -7.758121   4.848456 -1.6001  0.109572    
## lin1         -0.221467   0.056248 -3.9373 8.240e-05 ***
## lin2          0.211827   0.045021  4.7051 2.537e-06 ***
## postchange  -19.242855   6.843848 -2.8117  0.004928 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "14 Regression for xbxchib model2c on 2016-10-18/purple/10-30%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -107.310   -17.173     0.727    17.464    76.355  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  404.865      8.687  46.605  < 2e-16 ***
## hour1         -8.598     11.076  -0.776   0.4382    
## hour2        -64.384     11.076  -5.813 1.55e-08 ***
## hour3       -124.455     11.076 -11.237  < 2e-16 ***
## hour4       -156.527     11.076 -14.132  < 2e-16 ***
## hour5       -201.027     11.076 -18.150  < 2e-16 ***
## hour6       -249.598     11.076 -22.535  < 2e-16 ***
## hour7       -286.027     11.076 -25.824  < 2e-16 ***
## hour8       -315.384     11.076 -28.475  < 2e-16 ***
## hour9       -311.170     11.076 -28.094  < 2e-16 ***
## hour10      -315.384     11.076 -28.475  < 2e-16 ***
## hour11      -302.027     11.076 -27.269  < 2e-16 ***
## hour12      -266.455     11.076 -24.057  < 2e-16 ***
## hour13      -228.027     11.076 -20.588  < 2e-16 ***
## hour14      -196.670     11.076 -17.757  < 2e-16 ***
## hour15      -158.241     11.076 -14.287  < 2e-16 ***
## hour16      -136.884     11.076 -12.359  < 2e-16 ***
## hour17      -127.241     11.076 -11.488  < 2e-16 ***
## hour18      -102.384     11.076  -9.244  < 2e-16 ***
## hour19       -86.027     11.076  -7.767 1.23e-13 ***
## hour20       -83.313     11.076  -7.522 6.09e-13 ***
## hour21       -26.655     11.297  -2.359   0.0189 *  
## hour22       -16.116     11.297  -1.427   0.1547    
## hour23       -45.670     11.076  -4.123 4.82e-05 ***
## dow2         -38.708      6.082  -6.364 7.19e-10 ***
## dow3         -33.242      6.053  -5.492 8.39e-08 ***
## dow4         -51.896      6.082  -8.532 6.82e-16 ***
## dow5         -58.917      6.082  -9.687  < 2e-16 ***
## dow6         -42.001      6.153  -6.826 4.71e-11 ***
## dow7           7.896      6.082   1.298   0.1952    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 887.8292)
## 
##     Null deviance: 4102565  on 334  degrees of freedom
## Residual deviance:  270788  on 305  degrees of freedom
## AIC: 3255.5
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  404.8652     9.4339  42.9158 < 2.2e-16 ***
## hour1         -8.5984    13.4166  -0.6409 0.5216056    
## hour2        -64.3841    11.6308  -5.5356 3.101e-08 ***
## hour3       -124.4555    11.4425 -10.8766 < 2.2e-16 ***
## hour4       -156.5269     9.7145 -16.1127 < 2.2e-16 ***
## hour5       -201.0269    10.4772 -19.1871 < 2.2e-16 ***
## hour6       -249.5984    10.3749 -24.0578 < 2.2e-16 ***
## hour7       -286.0269    10.8007 -26.4823 < 2.2e-16 ***
## hour8       -315.3841    10.2377 -30.8063 < 2.2e-16 ***
## hour9       -311.1698    10.5718 -29.4340 < 2.2e-16 ***
## hour10      -315.3841    10.9826 -28.7166 < 2.2e-16 ***
## hour11      -302.0269    11.5475 -26.1552 < 2.2e-16 ***
## hour12      -266.4555     9.8674 -27.0036 < 2.2e-16 ***
## hour13      -228.0269     9.3817 -24.3054 < 2.2e-16 ***
## hour14      -196.6698    10.7433 -18.3062 < 2.2e-16 ***
## hour15      -158.2412    13.4857 -11.7340 < 2.2e-16 ***
## hour16      -136.8841    10.8479 -12.6185 < 2.2e-16 ***
## hour17      -127.2412    10.9970 -11.5705 < 2.2e-16 ***
## hour18      -102.3841    13.9595  -7.3343 2.228e-13 ***
## hour19       -86.0269    13.7330  -6.2643 3.746e-10 ***
## hour20       -83.3126    13.3176  -6.2558 3.954e-10 ***
## hour21       -26.6548    12.0997  -2.2029 0.0275999 *  
## hour22       -16.1164    15.6131  -1.0322 0.3019645    
## hour23       -45.6698    13.6859  -3.3370 0.0008469 ***
## dow2         -38.7083     5.6067  -6.9039 5.059e-12 ***
## dow3         -33.2421     6.8763  -4.8343 1.336e-06 ***
## dow4         -51.8958     5.6033  -9.2616 < 2.2e-16 ***
## dow5         -58.9167     6.7419  -8.7389 < 2.2e-16 ***
## dow6         -42.0009     5.4939  -7.6450 2.090e-14 ***
## dow7           7.8958     6.9361   1.1384 0.2549669    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -101.626   -16.128     0.338    16.634    72.001  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)  
## (Intercept) 10.41230    4.35900   2.389   0.0175 *
## lin1         0.10876    0.04474   2.431   0.0156 *
## lin2        -0.09910    0.04514  -2.195   0.0288 *
## postchange  -3.41632    6.14624  -0.556   0.5787  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 790.9338)
## 
##     Null deviance: 270788  on 334  degrees of freedom
## Residual deviance: 261799  on 331  degrees of freedom
## AIC: 3192.2
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error z value Pr(>|z|)   
## (Intercept) 10.412299   3.972957  2.6208 0.008773 **
## lin1         0.108758   0.044123  2.4649 0.013706 * 
## lin2        -0.099103   0.041402 -2.3937 0.016680 * 
## postchange  -3.416319   5.631930 -0.6066 0.544118   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "15 Regression for xbxchib model2c on 2016-11-01/purple/30-50%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -168.016   -27.546    -1.423    21.464   154.383  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  379.370     13.201  28.738  < 2e-16 ***
## hour1        -15.470     16.832  -0.919 0.358789    
## hour2        -69.398     16.832  -4.123 4.82e-05 ***
## hour3       -110.113     16.832  -6.542 2.54e-10 ***
## hour4       -137.255     16.832  -8.154 9.06e-15 ***
## hour5       -180.684     16.832 -10.734  < 2e-16 ***
## hour6       -228.398     16.832 -13.569  < 2e-16 ***
## hour7       -274.113     16.832 -16.285  < 2e-16 ***
## hour8       -291.398     16.832 -17.312  < 2e-16 ***
## hour9       -295.470     16.832 -17.554  < 2e-16 ***
## hour10      -295.827     16.832 -17.575  < 2e-16 ***
## hour11      -295.613     16.832 -17.562  < 2e-16 ***
## hour12      -266.541     16.832 -15.835  < 2e-16 ***
## hour13      -225.041     16.832 -13.370  < 2e-16 ***
## hour14      -188.184     16.832 -11.180  < 2e-16 ***
## hour15      -156.113     16.832  -9.275  < 2e-16 ***
## hour16      -132.041     16.832  -7.845 7.26e-14 ***
## hour17      -113.541     16.832  -6.745 7.59e-11 ***
## hour18       -91.398     16.832  -5.430 1.15e-07 ***
## hour19       -63.041     16.832  -3.745 0.000215 ***
## hour20       -51.898     16.832  -3.083 0.002233 ** 
## hour21       -16.755     16.832  -0.995 0.320308    
## hour22        -1.541     16.832  -0.092 0.927107    
## hour23         3.459     16.832   0.205 0.837327    
## dow2         -44.729      9.243  -4.839 2.07e-06 ***
## dow3         -19.946      9.199  -2.168 0.030902 *  
## dow4         -37.979      9.243  -4.109 5.11e-05 ***
## dow5         -44.833      9.243  -4.850 1.96e-06 ***
## dow6         -25.958      9.243  -2.808 0.005299 ** 
## dow7          22.646      9.243   2.450 0.014844 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 2050.505)
## 
##     Null deviance: 4484764  on 336  degrees of freedom
## Residual deviance:  629505  on 307  degrees of freedom
## AIC: 3556.9
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  379.3698    17.3044  21.9234 < 2.2e-16 ***
## hour1        -15.4698    22.1381  -0.6988 0.4846878    
## hour2        -69.3983    18.8798  -3.6758 0.0002371 ***
## hour3       -110.1126    21.5665  -5.1057 3.295e-07 ***
## hour4       -137.2555    20.6587  -6.6440 3.054e-11 ***
## hour5       -180.6840    18.6297  -9.6987 < 2.2e-16 ***
## hour6       -228.3983    19.6419 -11.6281 < 2.2e-16 ***
## hour7       -274.1126    18.7893 -14.5887 < 2.2e-16 ***
## hour8       -291.3983    16.5383 -17.6197 < 2.2e-16 ***
## hour9       -295.4698    16.8762 -17.5081 < 2.2e-16 ***
## hour10      -295.8269    16.9506 -17.4523 < 2.2e-16 ***
## hour11      -295.6126    16.8262 -17.5686 < 2.2e-16 ***
## hour12      -266.5412    16.9715 -15.7052 < 2.2e-16 ***
## hour13      -225.0412    17.4508 -12.8958 < 2.2e-16 ***
## hour14      -188.1840    17.9929 -10.4588 < 2.2e-16 ***
## hour15      -156.1126    20.5183  -7.6085 2.774e-14 ***
## hour16      -132.0412    19.5686  -6.7476 1.503e-11 ***
## hour17      -113.5412    20.4838  -5.5430 2.974e-08 ***
## hour18       -91.3983    21.8259  -4.1876 2.819e-05 ***
## hour19       -63.0412    20.3345  -3.1002 0.0019338 ** 
## hour20       -51.8983    19.7461  -2.6283 0.0085818 ** 
## hour21       -16.7555    22.0977  -0.7582 0.4483050    
## hour22        -1.5412    21.7656  -0.0708 0.9435506    
## hour23         3.4588    23.2553   0.1487 0.8817646    
## dow2         -44.7292     9.2791  -4.8204 1.433e-06 ***
## dow3         -19.9463     9.3000  -2.1448 0.0319721 *  
## dow4         -37.9792     8.0482  -4.7190 2.370e-06 ***
## dow5         -44.8333     8.3480  -5.3705 7.851e-08 ***
## dow6         -25.9583     9.9087  -2.6198 0.0087994 ** 
## dow7          22.6458    12.7580   1.7750 0.0758933 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -147.356   -22.000     1.483    19.239   128.361  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -44.33846    5.72354  -7.747 1.15e-13 ***
## lin1         -0.32887    0.05875  -5.598 4.54e-08 ***
## lin2          0.31901    0.05823   5.479 8.46e-08 ***
## postchange   33.99293    8.04669   4.224 3.10e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1363.625)
## 
##     Null deviance: 629505  on 336  degrees of freedom
## Residual deviance: 454087  on 333  degrees of freedom
## AIC: 3394.8
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error z value  Pr(>|z|)    
## (Intercept) -44.338456   5.143806 -8.6198 < 2.2e-16 ***
## lin1         -0.328865   0.050407 -6.5242 6.835e-11 ***
## lin2          0.319011   0.061668  5.1731 2.303e-07 ***
## postchange   33.992935   7.077758  4.8028 1.565e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "16 Regression for xbxchib model2c on 2016-12-15/purple/50-90%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -120.541   -22.457    -1.859    18.364   143.102  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  383.3987    11.7490  32.632  < 2e-16 ***
## hour1         -1.8788    14.9812  -0.125 0.900280    
## hour2        -31.8074    14.9812  -2.123 0.034540 *  
## hour3        -53.2359    14.9812  -3.554 0.000440 ***
## hour4       -100.3074    14.9812  -6.696 1.02e-10 ***
## hour5       -119.4502    14.9812  -7.973 3.07e-14 ***
## hour6       -165.2359    14.9812 -11.030  < 2e-16 ***
## hour7       -224.3788    14.9812 -14.977  < 2e-16 ***
## hour8       -261.0217    14.9812 -17.423  < 2e-16 ***
## hour9       -285.3788    14.9812 -19.049  < 2e-16 ***
## hour10      -286.1645    14.9812 -19.102  < 2e-16 ***
## hour11      -289.6645    14.9812 -19.335  < 2e-16 ***
## hour12      -275.9502    14.9812 -18.420  < 2e-16 ***
## hour13      -242.6645    14.9812 -16.198  < 2e-16 ***
## hour14      -206.6645    14.9812 -13.795  < 2e-16 ***
## hour15      -162.9502    14.9812 -10.877  < 2e-16 ***
## hour16      -128.4502    14.9812  -8.574 4.99e-16 ***
## hour17      -107.2359    14.9812  -7.158 6.06e-12 ***
## hour18       -84.0931    14.9812  -5.613 4.44e-08 ***
## hour19       -70.6645    14.9812  -4.717 3.64e-06 ***
## hour20       -39.8788    14.9812  -2.662 0.008179 ** 
## hour21       -13.5217    14.9812  -0.903 0.367457    
## hour22        10.2641    14.9812   0.685 0.493778    
## hour23        -0.2359    14.9812  -0.016 0.987444    
## dow2         -38.5625     8.2267  -4.687 4.17e-06 ***
## dow3         -33.8750     8.2267  -4.118 4.92e-05 ***
## dow4         -30.2083     8.2267  -3.672 0.000284 ***
## dow5         -76.7020     8.1873  -9.368  < 2e-16 ***
## dow6         -72.9792     8.2267  -8.871  < 2e-16 ***
## dow7           4.6875     8.2267   0.570 0.569238    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1624.301)
## 
##     Null deviance: 4386095  on 336  degrees of freedom
## Residual deviance:  498660  on 307  degrees of freedom
## AIC: 3478.3
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  383.39874   11.99472  31.9640 < 2.2e-16 ***
## hour1         -1.87881   19.54090  -0.0961 0.9234035    
## hour2        -31.80738   13.25677  -2.3993 0.0164251 *  
## hour3        -53.23595   15.46308  -3.4428 0.0005758 ***
## hour4       -100.30738   12.49816  -8.0258 1.009e-15 ***
## hour5       -119.45023   11.62435 -10.2759 < 2.2e-16 ***
## hour6       -165.23595   13.26378 -12.4577 < 2.2e-16 ***
## hour7       -224.37881   13.17943 -17.0249 < 2.2e-16 ***
## hour8       -261.02166   12.42414 -21.0092 < 2.2e-16 ***
## hour9       -285.37881   12.68090 -22.5046 < 2.2e-16 ***
## hour10      -286.16452   12.34804 -23.1749 < 2.2e-16 ***
## hour11      -289.66452   12.15535 -23.8302 < 2.2e-16 ***
## hour12      -275.95023   12.45153 -22.1619 < 2.2e-16 ***
## hour13      -242.66452   12.09972 -20.0554 < 2.2e-16 ***
## hour14      -206.66452   12.70927 -16.2609 < 2.2e-16 ***
## hour15      -162.95023   15.31636 -10.6390 < 2.2e-16 ***
## hour16      -128.45023   17.13425  -7.4967 6.545e-14 ***
## hour17      -107.23595   14.20984  -7.5466 4.468e-14 ***
## hour18       -84.09309   15.48329  -5.4312 5.597e-08 ***
## hour19       -70.66452   17.96728  -3.9330 8.391e-05 ***
## hour20       -39.87881   13.52922  -2.9476 0.0032025 ** 
## hour21       -13.52166   19.62020  -0.6892 0.4907160    
## hour22        10.26405   20.84236   0.4925 0.6223934    
## hour23        -0.23595   15.32370  -0.0154 0.9877150    
## dow2         -38.56250    8.98185  -4.2934 1.760e-05 ***
## dow3         -33.87500    7.76971  -4.3599 1.301e-05 ***
## dow4         -30.20833    6.59244  -4.5823 4.600e-06 ***
## dow5         -76.70202    7.46948 -10.2687 < 2.2e-16 ***
## dow6         -72.97917    9.23560  -7.9019 2.746e-15 ***
## dow7           4.68750    9.56043   0.4903 0.6239200    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -139.343   -22.564    -2.382    18.881   145.234  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -11.30298    5.42139  -2.085   0.0378 *  
## lin1          0.05922    0.05565   1.064   0.2880    
## lin2         -0.04392    0.05515  -0.796   0.4264    
## postchange   31.20286    7.62190   4.094 5.33e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1223.452)
## 
##     Null deviance: 498660  on 336  degrees of freedom
## Residual deviance: 407410  on 333  degrees of freedom
## AIC: 3358.2
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error z value  Pr(>|z|)    
## (Intercept) -11.302978   4.768488 -2.3703   0.01777 *  
## lin1          0.059225   0.045954  1.2888   0.19748    
## lin2         -0.043916   0.048835 -0.8993   0.36851    
## postchange   31.202856   7.454532  4.1858 2.842e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "17 Regression for xbxchib model2c on 2016-11-17/blue/content change"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -231.311   -23.464    -0.621    21.871   127.461  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   428.19      14.06  30.451  < 2e-16 ***
## hour1          12.54      16.97   0.739  0.46037    
## hour2         -30.23      16.97  -1.781  0.07592 .  
## hour3         -72.84      16.97  -4.293 2.42e-05 ***
## hour4        -131.46      16.97  -7.748 1.67e-13 ***
## hour5        -157.84      16.97  -9.303  < 2e-16 ***
## hour6        -220.46      16.97 -12.993  < 2e-16 ***
## hour7        -262.69      16.97 -15.482  < 2e-16 ***
## hour8        -279.92      16.97 -16.497  < 2e-16 ***
## hour9        -301.53      16.97 -17.771  < 2e-16 ***
## hour10       -307.69      16.97 -18.134  < 2e-16 ***
## hour11       -300.61      16.97 -17.717  < 2e-16 ***
## hour12       -291.92      16.97 -17.204  < 2e-16 ***
## hour13       -247.61      16.97 -14.593  < 2e-16 ***
## hour14       -206.84      16.97 -12.190  < 2e-16 ***
## hour15       -174.76      16.97 -10.300  < 2e-16 ***
## hour16       -138.94      16.66  -8.339 3.27e-15 ***
## hour17       -102.15      16.99  -6.014 5.56e-09 ***
## hour18        -98.69      16.99  -5.810 1.67e-08 ***
## hour19        -82.84      16.99  -4.877 1.79e-06 ***
## hour20        -52.84      16.99  -3.111  0.00205 ** 
## hour21        -28.46      16.99  -1.675  0.09494 .  
## hour22        -10.08      16.99  -0.593  0.55355    
## hour23        -14.61      16.99  -0.860  0.39033    
## dow2          -55.94      10.09  -5.544 6.74e-08 ***
## dow3          -69.55      10.09  -6.893 3.54e-11 ***
## dow4          -55.17      10.09  -5.468 9.98e-08 ***
## dow5          -93.13      10.05  -9.262  < 2e-16 ***
## dow6          -64.44      10.09  -6.387 6.90e-10 ***
## dow7          -32.04      10.50  -3.051  0.00249 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1939.499)
## 
##     Null deviance: 4439367  on 313  degrees of freedom
## Residual deviance:  550818  on 284  degrees of freedom
## AIC: 3298.6
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##             Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  428.190     13.947  30.7014 < 2.2e-16 ***
## hour1         12.543     16.649   0.7534  0.451228    
## hour2        -30.226     14.858  -2.0343  0.041924 *  
## hour3        -72.841     17.096  -4.2609 2.036e-05 ***
## hour4       -131.457     20.683  -6.3557 2.075e-10 ***
## hour5       -157.841     23.207  -6.8015 1.035e-11 ***
## hour6       -220.457     16.358 -13.4770 < 2.2e-16 ***
## hour7       -262.688     12.435 -21.1241 < 2.2e-16 ***
## hour8       -279.918     12.238 -22.8735 < 2.2e-16 ***
## hour9       -301.534     12.604 -23.9233 < 2.2e-16 ***
## hour10      -307.688     12.044 -25.5460 < 2.2e-16 ***
## hour11      -300.611     12.079 -24.8867 < 2.2e-16 ***
## hour12      -291.918     12.177 -23.9733 < 2.2e-16 ***
## hour13      -247.611     11.475 -21.5776 < 2.2e-16 ***
## hour14      -206.841     14.268 -14.4973 < 2.2e-16 ***
## hour15      -174.765     15.405 -11.3446 < 2.2e-16 ***
## hour16      -138.938     17.442  -7.9656 1.644e-15 ***
## hour17      -102.152     18.044  -5.6612 1.503e-08 ***
## hour18       -98.691     15.750  -6.2662 3.700e-10 ***
## hour19       -82.844     18.657  -4.4405 8.976e-06 ***
## hour20       -52.844     18.463  -2.8622  0.004208 ** 
## hour21       -28.460     17.101  -1.6642  0.096071 .  
## hour22       -10.075     18.466  -0.5456  0.585346    
## hour23       -14.614     17.275  -0.8460  0.397580    
## dow2         -55.941     10.536  -5.3096 1.098e-07 ***
## dow3         -69.546      9.840  -7.0677 1.575e-12 ***
## dow4         -55.171     10.511  -5.2488 1.531e-07 ***
## dow5         -93.130     10.470  -8.8947 < 2.2e-16 ***
## dow6         -64.442     10.518  -6.1266 8.975e-10 ***
## dow7         -32.038     15.632  -2.0494  0.040419 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -227.20   -23.52    -0.84    23.69   123.22  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)  
## (Intercept) -12.30730    6.95783  -1.769   0.0779 .
## lin1         -0.08908    0.08269  -1.077   0.2822  
## lin2          0.05705    0.06571   0.868   0.3860  
## postchange   12.49535    9.44230   1.323   0.1867  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1736.824)
## 
##     Null deviance: 550818  on 313  degrees of freedom
## Residual deviance: 538415  on 310  degrees of freedom
## AIC: 3239.4
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error z value Pr(>|z|)  
## (Intercept) -12.307302   5.150465 -2.3896  0.01687 *
## lin1         -0.089077   0.066256 -1.3444  0.17881  
## lin2          0.057051   0.048140  1.1851  0.23597  
## postchange   12.495352   7.547330  1.6556  0.09780 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "18 Regression for rate_xbxchib model2c on 2016-10-11/purple/0-10%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -1.288e-03  -2.448e-04  -2.723e-05   1.952e-04   1.728e-03  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.664e-03  1.310e-04  35.592  < 2e-16 ***
## hour1       -1.766e-05  1.663e-04  -0.106  0.91552    
## hour2       -4.374e-04  1.663e-04  -2.630  0.00898 ** 
## hour3       -8.037e-04  1.663e-04  -4.833 2.16e-06 ***
## hour4       -1.086e-03  1.663e-04  -6.527 2.89e-10 ***
## hour5       -1.573e-03  1.663e-04  -9.460  < 2e-16 ***
## hour6       -2.219e-03  1.663e-04 -13.341  < 2e-16 ***
## hour7       -2.761e-03  1.663e-04 -16.603  < 2e-16 ***
## hour8       -3.086e-03  1.696e-04 -18.194  < 2e-16 ***
## hour9       -3.077e-03  1.696e-04 -18.142  < 2e-16 ***
## hour10      -3.161e-03  1.696e-04 -18.638  < 2e-16 ***
## hour11      -3.237e-03  1.696e-04 -19.083  < 2e-16 ***
## hour12      -3.020e-03  1.663e-04 -18.157  < 2e-16 ***
## hour13      -2.898e-03  1.663e-04 -17.427  < 2e-16 ***
## hour14      -2.835e-03  1.734e-04 -16.348  < 2e-16 ***
## hour15      -2.573e-03  1.734e-04 -14.840  < 2e-16 ***
## hour16      -2.347e-03  1.663e-04 -14.114  < 2e-16 ***
## hour17      -2.191e-03  1.663e-04 -13.176  < 2e-16 ***
## hour18      -1.958e-03  1.663e-04 -11.771  < 2e-16 ***
## hour19      -1.900e-03  1.663e-04 -11.423  < 2e-16 ***
## hour20      -1.642e-03  1.696e-04  -9.682  < 2e-16 ***
## hour21      -9.660e-04  1.696e-04  -5.695 2.97e-08 ***
## hour22      -2.802e-04  1.663e-04  -1.685  0.09310 .  
## hour23      -3.230e-04  1.663e-04  -1.942  0.05305 .  
## dow2        -1.164e-03  9.239e-05 -12.595  < 2e-16 ***
## dow3        -1.167e-03  9.195e-05 -12.687  < 2e-16 ***
## dow4        -1.095e-03  9.239e-05 -11.848  < 2e-16 ***
## dow5        -1.170e-03  9.239e-05 -12.663  < 2e-16 ***
## dow6        -1.085e-03  9.581e-05 -11.325  < 2e-16 ***
## dow7         1.579e-04  9.329e-05   1.693  0.09149 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 2.001542e-07)
## 
##     Null deviance: 5.4694e-04  on 326  degrees of freedom
## Residual deviance: 5.9446e-05  on 297  degrees of freedom
## AIC: -4085.2
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  4.6637e-03  1.7250e-04  27.0355 < 2.2e-16 ***
## hour1       -1.7656e-05  2.3940e-04  -0.0737 0.9412097    
## hour2       -4.3739e-04  1.8209e-04  -2.4021 0.0163021 *  
## hour3       -8.0373e-04  2.2032e-04  -3.6480 0.0002643 ***
## hour4       -1.0855e-03  2.1283e-04  -5.1005 3.388e-07 ***
## hour5       -1.5732e-03  1.9457e-04  -8.0855 6.191e-16 ***
## hour6       -2.2186e-03  1.6640e-04 -13.3331 < 2.2e-16 ***
## hour7       -2.7611e-03  1.7028e-04 -16.2150 < 2.2e-16 ***
## hour8       -3.0862e-03  1.8496e-04 -16.6852 < 2.2e-16 ***
## hour9       -3.0774e-03  1.8992e-04 -16.2032 < 2.2e-16 ***
## hour10      -3.1614e-03  1.9097e-04 -16.5550 < 2.2e-16 ***
## hour11      -3.2370e-03  1.9105e-04 -16.9432 < 2.2e-16 ***
## hour12      -3.0196e-03  1.6495e-04 -18.3059 < 2.2e-16 ***
## hour13      -2.8980e-03  1.6134e-04 -17.9622 < 2.2e-16 ***
## hour14      -2.8349e-03  1.6161e-04 -17.5414 < 2.2e-16 ***
## hour15      -2.5733e-03  1.7978e-04 -14.3135 < 2.2e-16 ***
## hour16      -2.3473e-03  1.6761e-04 -14.0040 < 2.2e-16 ***
## hour17      -2.1912e-03  1.8852e-04 -11.6234 < 2.2e-16 ***
## hour18      -1.9575e-03  2.0635e-04  -9.4861 < 2.2e-16 ***
## hour19      -1.8996e-03  1.9540e-04  -9.7216 < 2.2e-16 ***
## hour20      -1.6423e-03  2.1583e-04  -7.6090 2.762e-14 ***
## hour21      -9.6599e-04  2.2417e-04  -4.3092 1.639e-05 ***
## hour22      -2.8017e-04  2.3037e-04  -1.2162 0.2239144    
## hour23      -3.2300e-04  2.5354e-04  -1.2739 0.2026813    
## dow2        -1.1636e-03  8.0568e-05 -14.4422 < 2.2e-16 ***
## dow3        -1.1666e-03  8.4902e-05 -13.7404 < 2.2e-16 ***
## dow4        -1.0946e-03  9.4599e-05 -11.5714 < 2.2e-16 ***
## dow5        -1.1699e-03  9.6937e-05 -12.0686 < 2.2e-16 ***
## dow6        -1.0850e-03  9.3122e-05 -11.6518 < 2.2e-16 ***
## dow7         1.5794e-04  1.1841e-04   1.3339 0.1822497    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -1.391e-03  -2.507e-04  -3.394e-05   1.834e-04   1.572e-03  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)  
## (Intercept)  2.178e-05  6.545e-05   0.333   0.7395  
## lin1        -1.214e-06  7.141e-07  -1.700   0.0901 .
## lin2         1.158e-06  6.456e-07   1.793   0.0738 .
## postchange  -2.296e-04  9.065e-05  -2.533   0.0118 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1.676253e-07)
## 
##     Null deviance: 5.9446e-05  on 326  degrees of freedom
## Residual deviance: 5.4143e-05  on 323  degrees of freedom
## AIC: -4167.7
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)   
## (Intercept)  2.1783e-05  6.4600e-05  0.3372 0.735966   
## lin1        -1.2142e-06  6.1241e-07 -1.9826 0.047413 * 
## lin2         1.1578e-06  5.6435e-07  2.0515 0.040219 * 
## postchange  -2.2964e-04  8.3629e-05 -2.7460 0.006033 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "19 Regression for rate_xbxchib model2c on 2016-10-18/purple/10-30%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##       Min         1Q     Median         3Q        Max  
## -0.005581  -0.000388   0.000199   0.000471   0.045578  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.469e-03  1.071e-03   4.173 3.92e-05 ***
## hour1       -2.871e-05  1.365e-03  -0.021   0.9832    
## hour2       -4.663e-04  1.365e-03  -0.342   0.7329    
## hour3       -9.245e-04  1.365e-03  -0.677   0.4988    
## hour4       -1.145e-03  1.365e-03  -0.839   0.4024    
## hour5       -1.601e-03  1.365e-03  -1.172   0.2419    
## hour6       -2.277e-03  1.365e-03  -1.668   0.0964 .  
## hour7       -2.763e-03  1.365e-03  -2.024   0.0439 *  
## hour8       -3.101e-03  1.365e-03  -2.271   0.0238 *  
## hour9       -3.083e-03  1.365e-03  -2.258   0.0247 *  
## hour10      -3.190e-03  1.365e-03  -2.336   0.0201 *  
## hour11      -3.209e-03  1.365e-03  -2.350   0.0194 *  
## hour12      -3.087e-03  1.365e-03  -2.261   0.0245 *  
## hour13      -2.939e-03  1.365e-03  -2.153   0.0321 *  
## hour14      -2.751e-03  1.365e-03  -2.015   0.0448 *  
## hour15      -2.526e-03  1.365e-03  -1.850   0.0653 .  
## hour16      -2.376e-03  1.365e-03  -1.741   0.0828 .  
## hour17      -2.288e-03  1.365e-03  -1.675   0.0949 .  
## hour18      -1.602e-03  1.365e-03  -1.173   0.2416    
## hour19       1.239e-03  1.365e-03   0.907   0.3650    
## hour20       1.934e-03  1.365e-03   1.417   0.1575    
## hour21      -8.546e-04  1.393e-03  -0.614   0.5399    
## hour22      -3.301e-04  1.393e-03  -0.237   0.8128    
## hour23      -5.239e-06  1.365e-03  -0.004   0.9969    
## dow2        -1.402e-03  7.497e-04  -1.870   0.0624 .  
## dow3        -1.232e-03  7.462e-04  -1.651   0.0997 .  
## dow4        -1.407e-03  7.497e-04  -1.877   0.0615 .  
## dow5        -1.536e-03  7.497e-04  -2.049   0.0414 *  
## dow6         1.098e-03  7.584e-04   1.448   0.1487    
## dow7         2.829e-05  7.497e-04   0.038   0.9699    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1.349078e-05)
## 
##     Null deviance: 0.0051101  on 334  degrees of freedom
## Residual deviance: 0.0041147  on 305  degrees of freedom
## AIC: -2775.3
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value  Pr(>|z|)    
## (Intercept)  4.4692e-03  3.5185e-04 12.7022 < 2.2e-16 ***
## hour1       -2.8710e-05  3.8414e-04 -0.0747 0.9404233    
## hour2       -4.6635e-04  3.6160e-04 -1.2897 0.1971557    
## hour3       -9.2448e-04  3.8573e-04 -2.3967 0.0165440 *  
## hour4       -1.1449e-03  3.5801e-04 -3.1981 0.0013834 ** 
## hour5       -1.6007e-03  3.6510e-04 -4.3843 1.164e-05 ***
## hour6       -2.2772e-03  3.5537e-04 -6.4079 1.476e-10 ***
## hour7       -2.7629e-03  3.5210e-04 -7.8467 4.271e-15 ***
## hour8       -3.1010e-03  3.4660e-04 -8.9469 < 2.2e-16 ***
## hour9       -3.0826e-03  3.6272e-04 -8.4986 < 2.2e-16 ***
## hour10      -3.1897e-03  3.6325e-04 -8.7809 < 2.2e-16 ***
## hour11      -3.2090e-03  3.6764e-04 -8.7288 < 2.2e-16 ***
## hour12      -3.0866e-03  3.6051e-04 -8.5619 < 2.2e-16 ***
## hour13      -2.9391e-03  3.5644e-04 -8.2457 < 2.2e-16 ***
## hour14      -2.7513e-03  3.6675e-04 -7.5019 6.290e-14 ***
## hour15      -2.5257e-03  3.7754e-04 -6.6900 2.231e-11 ***
## hour16      -2.3763e-03  3.6480e-04 -6.5140 7.315e-11 ***
## hour17      -2.2875e-03  3.6261e-04 -6.3085 2.818e-10 ***
## hour18      -1.6018e-03  4.5584e-04 -3.5141 0.0004413 ***
## hour19       1.2386e-03  2.9870e-03  0.4147 0.6783823    
## hour20       1.9345e-03  3.5470e-03  0.5454 0.5854976    
## hour21      -8.5460e-04  3.7441e-04 -2.2825 0.0224578 *  
## hour22      -3.3015e-04  3.8050e-04 -0.8677 0.3855826    
## hour23      -5.2395e-06  5.5741e-04 -0.0094 0.9925002    
## dow2        -1.4021e-03  2.0537e-04 -6.8273 8.651e-12 ***
## dow3        -1.2323e-03  2.1964e-04 -5.6104 2.018e-08 ***
## dow4        -1.4070e-03  2.0596e-04 -6.8317 8.391e-12 ***
## dow5        -1.5359e-03  2.1108e-04 -7.2761 3.435e-13 ***
## dow6         1.0980e-03  1.4114e-03  0.7779 0.4366110    
## dow7         2.8290e-05  1.9323e-04  0.1464 0.8836031    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##       Min         1Q     Median         3Q        Max  
## -0.005215  -0.000566   0.000094   0.000600   0.045208  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)
## (Intercept) -3.998e-04  5.435e-04  -0.736    0.462
## lin1        -4.466e-07  5.579e-06  -0.080    0.936
## lin2         6.440e-07  5.629e-06   0.114    0.909
## postchange   7.107e-04  7.664e-04   0.927    0.354
## 
## (Dispersion parameter for gaussian family taken to be 1.229686e-05)
## 
##     Null deviance: 0.0041147  on 334  degrees of freedom
## Residual deviance: 0.0040703  on 331  degrees of freedom
## AIC: -2830.9
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)  
## (Intercept) -3.9984e-04  1.8457e-04 -2.1663  0.03029 *
## lin1        -4.4663e-07  1.7689e-06 -0.2525  0.80066  
## lin2         6.4401e-07  2.1633e-06  0.2977  0.76593  
## postchange   7.1066e-04  3.5190e-04  2.0195  0.04344 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "20 Regression for rate_xbxchib model2c on 2016-11-01/purple/30-50%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -2.320e-03  -3.017e-04  -3.385e-05   2.250e-04   2.158e-03  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.0047393  0.0001715  27.633  < 2e-16 ***
## hour1       -0.0001227  0.0002187  -0.561 0.575159    
## hour2       -0.0005298  0.0002187  -2.423 0.015989 *  
## hour3       -0.0007750  0.0002187  -3.544 0.000456 ***
## hour4       -0.0009341  0.0002187  -4.271 2.59e-05 ***
## hour5       -0.0013643  0.0002187  -6.239 1.46e-09 ***
## hour6       -0.0020765  0.0002187  -9.495  < 2e-16 ***
## hour7       -0.0027135  0.0002187 -12.408  < 2e-16 ***
## hour8       -0.0029465  0.0002187 -13.474  < 2e-16 ***
## hour9       -0.0030528  0.0002187 -13.960  < 2e-16 ***
## hour10      -0.0031161  0.0002187 -14.249  < 2e-16 ***
## hour11      -0.0032365  0.0002187 -14.799  < 2e-16 ***
## hour12      -0.0031576  0.0002187 -14.439  < 2e-16 ***
## hour13      -0.0029841  0.0002187 -13.645  < 2e-16 ***
## hour14      -0.0028065  0.0002187 -12.833  < 2e-16 ***
## hour15      -0.0025958  0.0002187 -11.870  < 2e-16 ***
## hour16      -0.0024052  0.0002187 -10.998  < 2e-16 ***
## hour17      -0.0022098  0.0002187 -10.105  < 2e-16 ***
## hour18      -0.0020428  0.0002187  -9.341  < 2e-16 ***
## hour19      -0.0017844  0.0002187  -8.160 8.74e-15 ***
## hour20      -0.0015745  0.0002187  -7.200 4.68e-12 ***
## hour21      -0.0010302  0.0002187  -4.711 3.74e-06 ***
## hour22      -0.0005108  0.0002187  -2.336 0.020145 *  
## hour23      -0.0001203  0.0002187  -0.550 0.582541    
## dow2        -0.0014134  0.0001201 -11.770  < 2e-16 ***
## dow3        -0.0012047  0.0001195 -10.080  < 2e-16 ***
## dow4        -0.0013058  0.0001201 -10.873  < 2e-16 ***
## dow5        -0.0014544  0.0001201 -12.111  < 2e-16 ***
## dow6        -0.0011582  0.0001201  -9.645  < 2e-16 ***
## dow7         0.0001855  0.0001201   1.545 0.123478    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 3.461259e-07)
## 
##     Null deviance: 0.00063658  on 336  degrees of freedom
## Residual deviance: 0.00010626  on 307  degrees of freedom
## AIC: -4026.4
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  0.00473934  0.00026966  17.5752 < 2.2e-16 ***
## hour1       -0.00012270  0.00032331  -0.3795 0.7043105    
## hour2       -0.00052979  0.00027131  -1.9528 0.0508491 .  
## hour3       -0.00077499  0.00030260  -2.5611 0.0104336 *  
## hour4       -0.00093413  0.00029842  -3.1303 0.0017465 ** 
## hour5       -0.00136433  0.00028624  -4.7663 1.876e-06 ***
## hour6       -0.00207648  0.00028274  -7.3442 2.070e-13 ***
## hour7       -0.00271354  0.00027133 -10.0008 < 2.2e-16 ***
## hour8       -0.00294653  0.00025489 -11.5600 < 2.2e-16 ***
## hour9       -0.00305283  0.00026772 -11.4031 < 2.2e-16 ***
## hour10      -0.00311607  0.00026966 -11.5555 < 2.2e-16 ***
## hour11      -0.00323645  0.00026236 -12.3358 < 2.2e-16 ***
## hour12      -0.00315764  0.00025084 -12.5882 < 2.2e-16 ***
## hour13      -0.00298408  0.00024271 -12.2949 < 2.2e-16 ***
## hour14      -0.00280646  0.00024243 -11.5765 < 2.2e-16 ***
## hour15      -0.00259577  0.00025710 -10.0962 < 2.2e-16 ***
## hour16      -0.00240519  0.00025579  -9.4031 < 2.2e-16 ***
## hour17      -0.00220982  0.00026885  -8.2195 < 2.2e-16 ***
## hour18      -0.00204283  0.00029010  -7.0419 1.896e-12 ***
## hour19      -0.00178441  0.00028272  -6.3116 2.761e-10 ***
## hour20      -0.00157447  0.00028194  -5.5845 2.344e-08 ***
## hour21      -0.00103021  0.00030976  -3.3259 0.0008815 ***
## hour22      -0.00051082  0.00030035  -1.7007 0.0889975 .  
## hour23      -0.00012034  0.00030929  -0.3891 0.6972256    
## dow2        -0.00141342  0.00013151 -10.7477 < 2.2e-16 ***
## dow3        -0.00120469  0.00012848  -9.3763 < 2.2e-16 ***
## dow4        -0.00130579  0.00012064 -10.8241 < 2.2e-16 ***
## dow5        -0.00145445  0.00012334 -11.7925 < 2.2e-16 ***
## dow6        -0.00115824  0.00013839  -8.3691 < 2.2e-16 ***
## dow7         0.00018549  0.00018566   0.9991 0.3177631    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -2.064e-03  -2.622e-04  -1.408e-05   2.567e-04   1.838e-03  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -5.057e-04  7.689e-05  -6.577 1.86e-10 ***
## lin1        -3.461e-06  7.892e-07  -4.385 1.56e-05 ***
## lin2         3.283e-06  7.822e-07   4.197 3.47e-05 ***
## postchange   4.420e-04  1.081e-04   4.089 5.45e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 2.460993e-07)
## 
##     Null deviance: 1.0626e-04  on 336  degrees of freedom
## Residual deviance: 8.1951e-05  on 333  degrees of freedom
## AIC: -4166
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value  Pr(>|z|)    
## (Intercept) -5.0574e-04  6.6690e-05 -7.5833 3.368e-14 ***
## lin1        -3.4608e-06  6.1225e-07 -5.6526 1.580e-08 ***
## lin2         3.2833e-06  8.0996e-07  4.0537 5.042e-05 ***
## postchange   4.4197e-04  8.6305e-05  5.1210 3.039e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "21 Regression for rate_xbxchib model2c on 2016-12-15/purple/50-90%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0014191  -0.0003134   0.0000006   0.0002497   0.0018290  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  4.776e-03  1.522e-04  31.377  < 2e-16 ***
## hour1       -5.298e-05  1.941e-04  -0.273   0.7851    
## hour2       -2.476e-04  1.941e-04  -1.276   0.2030    
## hour3       -3.702e-04  1.941e-04  -1.907   0.0574 .  
## hour4       -7.759e-04  1.941e-04  -3.997 8.02e-05 ***
## hour5       -8.339e-04  1.941e-04  -4.296 2.33e-05 ***
## hour6       -1.374e-03  1.941e-04  -7.077 1.00e-11 ***
## hour7       -2.096e-03  1.941e-04 -10.798  < 2e-16 ***
## hour8       -2.580e-03  1.941e-04 -13.294  < 2e-16 ***
## hour9       -2.912e-03  1.941e-04 -15.001  < 2e-16 ***
## hour10      -2.936e-03  1.941e-04 -15.125  < 2e-16 ***
## hour11      -3.042e-03  1.941e-04 -15.670  < 2e-16 ***
## hour12      -3.046e-03  1.941e-04 -15.695  < 2e-16 ***
## hour13      -2.950e-03  1.941e-04 -15.196  < 2e-16 ***
## hour14      -2.737e-03  1.941e-04 -14.100  < 2e-16 ***
## hour15      -2.525e-03  1.941e-04 -13.008  < 2e-16 ***
## hour16      -2.357e-03  1.941e-04 -12.141  < 2e-16 ***
## hour17      -2.123e-03  1.941e-04 -10.936  < 2e-16 ***
## hour18      -1.887e-03  1.941e-04  -9.725  < 2e-16 ***
## hour19      -1.758e-03  1.941e-04  -9.055  < 2e-16 ***
## hour20      -1.455e-03  1.941e-04  -7.498 7.02e-13 ***
## hour21      -1.097e-03  1.941e-04  -5.653 3.60e-08 ***
## hour22      -5.584e-04  1.941e-04  -2.877   0.0043 ** 
## hour23      -2.014e-04  1.941e-04  -1.038   0.3002    
## dow2        -1.404e-03  1.066e-04 -13.175  < 2e-16 ***
## dow3        -1.262e-03  1.066e-04 -11.842  < 2e-16 ***
## dow4        -1.173e-03  1.066e-04 -11.006  < 2e-16 ***
## dow5        -1.733e-03  1.061e-04 -16.339  < 2e-16 ***
## dow6        -1.585e-03  1.066e-04 -14.868  < 2e-16 ***
## dow7        -1.213e-04  1.066e-04  -1.139   0.2558    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 2.726478e-07)
## 
##     Null deviance: 5.9683e-04  on 336  degrees of freedom
## Residual deviance: 8.3703e-05  on 307  degrees of freedom
## AIC: -4106.8
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  4.7762e-03  2.1398e-04  22.3208 < 2.2e-16 ***
## hour1       -5.2975e-05  3.0361e-04  -0.1745 0.8614866    
## hour2       -2.4761e-04  2.3359e-04  -1.0600 0.2891475    
## hour3       -3.7023e-04  2.5573e-04  -1.4478 0.1476862    
## hour4       -7.7589e-04  2.1779e-04  -3.5625 0.0003673 ***
## hour5       -8.3391e-04  2.0759e-04  -4.0172 5.890e-05 ***
## hour6       -1.3737e-03  2.1358e-04  -6.4315 1.263e-10 ***
## hour7       -2.0958e-03  2.1938e-04  -9.5530 < 2.2e-16 ***
## hour8       -2.5803e-03  2.1584e-04 -11.9547 < 2.2e-16 ***
## hour9       -2.9117e-03  2.2238e-04 -13.0930 < 2.2e-16 ***
## hour10      -2.9357e-03  2.2351e-04 -13.1343 < 2.2e-16 ***
## hour11      -3.0415e-03  2.2371e-04 -13.5959 < 2.2e-16 ***
## hour12      -3.0464e-03  2.1508e-04 -14.1642 < 2.2e-16 ***
## hour13      -2.9495e-03  2.0273e-04 -14.5489 < 2.2e-16 ***
## hour14      -2.7368e-03  2.0118e-04 -13.6038 < 2.2e-16 ***
## hour15      -2.5248e-03  2.2031e-04 -11.4602 < 2.2e-16 ***
## hour16      -2.3565e-03  2.0398e-04 -11.5530 < 2.2e-16 ***
## hour17      -2.1225e-03  2.0793e-04 -10.2080 < 2.2e-16 ***
## hour18      -1.8875e-03  2.2307e-04  -8.4615 < 2.2e-16 ***
## hour19      -1.7575e-03  2.4311e-04  -7.2294 4.853e-13 ***
## hour20      -1.4552e-03  2.2299e-04  -6.5261 6.749e-11 ***
## hour21      -1.0973e-03  2.5231e-04  -4.3490 1.367e-05 ***
## hour22      -5.5838e-04  2.9551e-04  -1.8896 0.0588177 .  
## hour23      -2.0141e-04  2.6005e-04  -0.7745 0.4386250    
## dow2        -1.4042e-03  1.1772e-04 -11.9290 < 2.2e-16 ***
## dow3        -1.2622e-03  1.0976e-04 -11.4999 < 2.2e-16 ***
## dow4        -1.1730e-03  1.0403e-04 -11.2762 < 2.2e-16 ***
## dow5        -1.7332e-03  1.0852e-04 -15.9706 < 2.2e-16 ***
## dow6        -1.5847e-03  1.2175e-04 -13.0159 < 2.2e-16 ***
## dow7        -1.2135e-04  1.5032e-04  -0.8073 0.4195023    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -1.561e-03  -2.612e-04  -4.514e-05   2.555e-04   1.592e-03  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -2.447e-04  6.795e-05  -3.602 0.000364 ***
## lin1        -2.868e-08  6.974e-07  -0.041 0.967228    
## lin2         2.844e-07  6.913e-07   0.411 0.681044    
## postchange   4.617e-04  9.553e-05   4.833 2.05e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 1.921934e-07)
## 
##     Null deviance: 8.3703e-05  on 336  degrees of freedom
## Residual deviance: 6.4000e-05  on 333  degrees of freedom
## AIC: -4249.3
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value  Pr(>|z|)    
## (Intercept) -2.4473e-04  4.6249e-05 -5.2916 1.212e-07 ***
## lin1        -2.8675e-08  5.1578e-07 -0.0556    0.9557    
## lin2         2.8438e-07  5.6387e-07  0.5043    0.6140    
## postchange   4.6172e-04  8.4525e-05  5.4625 4.694e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "22 Regression for rate_xbxchib model2c on 2016-11-17/blue/content change"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0034923  -0.0003257   0.0000144   0.0002926   0.0029251  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.0054945  0.0001958  28.065  < 2e-16 ***
## hour1       -0.0001577  0.0002362  -0.668  0.50492    
## hour2       -0.0004189  0.0002362  -1.773  0.07726 .  
## hour3       -0.0007599  0.0002362  -3.217  0.00145 ** 
## hour4       -0.0012746  0.0002362  -5.395 1.44e-07 ***
## hour5       -0.0014421  0.0002362  -6.104 3.38e-09 ***
## hour6       -0.0022199  0.0002362  -9.397  < 2e-16 ***
## hour7       -0.0027421  0.0002362 -11.607  < 2e-16 ***
## hour8       -0.0029578  0.0002362 -12.520  < 2e-16 ***
## hour9       -0.0032628  0.0002362 -13.812  < 2e-16 ***
## hour10      -0.0033575  0.0002362 -14.212  < 2e-16 ***
## hour11      -0.0033185  0.0002362 -14.047  < 2e-16 ***
## hour12      -0.0034321  0.0002362 -14.528  < 2e-16 ***
## hour13      -0.0032682  0.0002362 -13.834  < 2e-16 ***
## hour14      -0.0030477  0.0002362 -12.901  < 2e-16 ***
## hour15      -0.0029108  0.0002362 -12.322  < 2e-16 ***
## hour16      -0.0026499  0.0002320 -11.423  < 2e-16 ***
## hour17      -0.0023476  0.0002365  -9.927  < 2e-16 ***
## hour18      -0.0022794  0.0002365  -9.638  < 2e-16 ***
## hour19      -0.0021265  0.0002365  -8.992  < 2e-16 ***
## hour20      -0.0018070  0.0002365  -7.641 3.34e-13 ***
## hour21      -0.0014539  0.0002365  -6.148 2.66e-09 ***
## hour22      -0.0009610  0.0002365  -4.064 6.26e-05 ***
## hour23      -0.0004190  0.0002365  -1.772  0.07754 .  
## dow2        -0.0016460  0.0001405 -11.717  < 2e-16 ***
## dow3        -0.0016973  0.0001405 -12.082  < 2e-16 ***
## dow4        -0.0014893  0.0001405 -10.601  < 2e-16 ***
## dow5        -0.0019589  0.0001400 -13.993  < 2e-16 ***
## dow6        -0.0015085  0.0001405 -10.738  < 2e-16 ***
## dow7        -0.0004522  0.0001462  -3.094  0.00218 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 3.759669e-07)
## 
##     Null deviance: 0.00060926  on 313  degrees of freedom
## Residual deviance: 0.00010677  on 284  degrees of freedom
## AIC: -3723.7
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  0.00549447  0.00030651  17.9261 < 2.2e-16 ***
## hour1       -0.00015772  0.00031850  -0.4952 0.6204635    
## hour2       -0.00041890  0.00030947  -1.3536 0.1758650    
## hour3       -0.00075988  0.00032912  -2.3088 0.0209524 *  
## hour4       -0.00127459  0.00036768  -3.4665 0.0005272 ***
## hour5       -0.00144207  0.00041703  -3.4580 0.0005443 ***
## hour6       -0.00221992  0.00034415  -6.4505 1.115e-10 ***
## hour7       -0.00274205  0.00030237  -9.0687 < 2.2e-16 ***
## hour8       -0.00295779  0.00029297 -10.0960 < 2.2e-16 ***
## hour9       -0.00326281  0.00030653 -10.6443 < 2.2e-16 ***
## hour10      -0.00335745  0.00030539 -10.9941 < 2.2e-16 ***
## hour11      -0.00331845  0.00029749 -11.1547 < 2.2e-16 ***
## hour12      -0.00343207  0.00030241 -11.3491 < 2.2e-16 ***
## hour13      -0.00326821  0.00028788 -11.3528 < 2.2e-16 ***
## hour14      -0.00304769  0.00029379 -10.3738 < 2.2e-16 ***
## hour15      -0.00291084  0.00029835  -9.7564 < 2.2e-16 ***
## hour16      -0.00264989  0.00030568  -8.6689 < 2.2e-16 ***
## hour17      -0.00234756  0.00031031  -7.5653 3.869e-14 ***
## hour18      -0.00227940  0.00030435  -7.4894 6.916e-14 ***
## hour19      -0.00212653  0.00032651  -6.5129 7.370e-11 ***
## hour20      -0.00180696  0.00034107  -5.2979 1.171e-07 ***
## hour21      -0.00145390  0.00034234  -4.2470 2.167e-05 ***
## hour22      -0.00096104  0.00035132  -2.7355 0.0062288 ** 
## hour23      -0.00041896  0.00037273  -1.1240 0.2609944    
## dow2        -0.00164598  0.00014489 -11.3602 < 2.2e-16 ***
## dow3        -0.00169729  0.00013497 -12.5755 < 2.2e-16 ***
## dow4        -0.00148927  0.00013922 -10.6973 < 2.2e-16 ***
## dow5        -0.00195887  0.00014070 -13.9223 < 2.2e-16 ***
## dow6        -0.00150852  0.00014714 -10.2523 < 2.2e-16 ***
## dow7        -0.00045221  0.00023765  -1.9029 0.0570549 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0034444  -0.0003097   0.0000039   0.0002961   0.0028819  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.064e-04  9.748e-05  -1.092    0.276
## lin1        -6.360e-07  1.158e-06  -0.549    0.583
## lin2         2.316e-07  9.206e-07   0.252    0.802
## postchange   1.385e-04  1.323e-04   1.047    0.296
## 
## (Dispersion parameter for gaussian family taken to be 3.409033e-07)
## 
##     Null deviance: 0.00010677  on 313  degrees of freedom
## Residual deviance: 0.00010568  on 310  degrees of freedom
## AIC: -3778.9
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)  
## (Intercept) -1.0644e-04  5.2936e-05 -2.0107  0.04436 *
## lin1        -6.3599e-07  7.4290e-07 -0.8561  0.39195  
## lin2         2.3164e-07  6.3496e-07  0.3648  0.71525  
## postchange   1.3847e-04  9.8049e-05  1.4122  0.15788  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "23 Regression for xbxcall model2c on 2016-10-11/purple/0-10%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -861.18  -100.15   -12.90    60.18  1365.53  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1697.833     61.051  27.810  < 2e-16 ***
## hour1        -287.317     77.485  -3.708 0.000249 ***
## hour2        -718.674     77.485  -9.275  < 2e-16 ***
## hour3        -876.674     77.485 -11.314  < 2e-16 ***
## hour4       -1179.531     77.485 -15.223  < 2e-16 ***
## hour5       -1440.460     77.485 -18.590  < 2e-16 ***
## hour6       -1564.960     77.485 -20.197  < 2e-16 ***
## hour7       -1619.317     77.485 -20.899  < 2e-16 ***
## hour8       -1623.478     79.034 -20.542  < 2e-16 ***
## hour9       -1638.632     79.034 -20.733  < 2e-16 ***
## hour10      -1587.863     79.034 -20.091  < 2e-16 ***
## hour11      -1514.401     79.034 -19.161  < 2e-16 ***
## hour12      -1289.174     77.485 -16.638  < 2e-16 ***
## hour13       -898.745     77.485 -11.599  < 2e-16 ***
## hour14       -681.049     80.794  -8.429 1.52e-15 ***
## hour15       -450.383     80.794  -5.574 5.58e-08 ***
## hour16       -258.388     77.485  -3.335 0.000962 ***
## hour17       -218.031     77.485  -2.814 0.005222 ** 
## hour18       -196.960     77.485  -2.542 0.011533 *  
## hour19         16.469     77.485   0.213 0.831827    
## hour20        209.753     79.034   2.654 0.008384 ** 
## hour21        171.137     79.034   2.165 0.031155 *  
## hour22        279.398     77.485   3.606 0.000365 ***
## hour23        159.398     77.485   2.057 0.040545 *  
## dow2           -7.588     43.046  -0.176 0.860193    
## dow3           18.733     42.842   0.437 0.662234    
## dow4          -20.630     43.046  -0.479 0.632112    
## dow5          -33.692     43.046  -0.783 0.434424    
## dow6           48.388     44.640   1.084 0.279266    
## dow7          116.174     43.465   2.673 0.007937 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 43451.58)
## 
##     Null deviance: 165708111  on 326  degrees of freedom
## Residual deviance:  12905121  on 297  degrees of freedom
## AIC: 4450.7
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  1697.8332    71.4728  23.7550 < 2.2e-16 ***
## hour1        -287.3167   113.8586  -2.5235  0.011621 *  
## hour2        -718.6738    75.6349  -9.5019 < 2.2e-16 ***
## hour3        -876.6738   130.0471  -6.7412 1.571e-11 ***
## hour4       -1179.5310    83.4696 -14.1313 < 2.2e-16 ***
## hour5       -1440.4595    67.3790 -21.3785 < 2.2e-16 ***
## hour6       -1564.9595    67.1309 -23.3121 < 2.2e-16 ***
## hour7       -1619.3167    67.6285 -23.9443 < 2.2e-16 ***
## hour8       -1623.4779    75.8873 -21.3933 < 2.2e-16 ***
## hour9       -1638.6317    74.9192 -21.8720 < 2.2e-16 ***
## hour10      -1587.8625    83.3692 -19.0462 < 2.2e-16 ***
## hour11      -1514.4010    71.3007 -21.2396 < 2.2e-16 ***
## hour12      -1289.1738    69.7182 -18.4912 < 2.2e-16 ***
## hour13       -898.7452    68.8102 -13.0612 < 2.2e-16 ***
## hour14       -681.0494    83.4096  -8.1651 3.211e-16 ***
## hour15       -450.3828    89.2483  -5.0464 4.502e-07 ***
## hour16       -258.3881    94.5524  -2.7328  0.006281 ** 
## hour17       -218.0310    86.2698  -2.5273  0.011494 *  
## hour18       -196.9595   105.5561  -1.8659  0.062052 .  
## hour19         16.4690    82.5781   0.1994  0.841922    
## hour20        209.7529   100.2490   2.0923  0.036410 *  
## hour21        171.1375    84.4228   2.0271  0.042647 *  
## hour22        279.3976    91.5569   3.0516  0.002276 ** 
## hour23        159.3976    91.0781   1.7501  0.080098 .  
## dow2           -7.5882    36.8962  -0.2057  0.837053    
## dow3           18.7333    36.2287   0.5171  0.605096    
## dow4          -20.6299    44.1112  -0.4677  0.640014    
## dow5          -33.6924    49.0215  -0.6873  0.491895    
## dow6           48.3877    40.9510   1.1816  0.237364    
## dow7          116.1739    37.2083   3.1223  0.001795 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -863.55   -96.59   -11.32    58.88  1362.13  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)
## (Intercept)  -4.0118    31.9303  -0.126    0.900
## lin1         -0.0693     0.3484  -0.199    0.842
## lin2          0.2124     0.3149   0.674    0.501
## postchange  -15.2302    44.2196  -0.344    0.731
## 
## (Dispersion parameter for gaussian family taken to be 39890.74)
## 
##     Null deviance: 12905121  on 326  degrees of freedom
## Residual deviance: 12884710  on 323  degrees of freedom
## AIC: 4398.2
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error z value Pr(>|z|)
## (Intercept)  -4.011753  27.703958 -0.1448   0.8849
## lin1         -0.069299   0.352631 -0.1965   0.8442
## lin2          0.212403   0.280329  0.7577   0.4486
## postchange  -15.230185  37.540201 -0.4057   0.6850

## [1] "24 Regression for xbxcall model2c on 2016-10-18/purple/10-30%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -452.39   -98.06    -6.70    64.47   770.39  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1699.04      47.52  35.752  < 2e-16 ***
## hour1        -324.08      60.59  -5.349 1.74e-07 ***
## hour2        -692.65      60.59 -11.432  < 2e-16 ***
## hour3        -968.08      60.59 -15.978  < 2e-16 ***
## hour4       -1213.73      60.59 -20.032  < 2e-16 ***
## hour5       -1430.94      60.59 -23.617  < 2e-16 ***
## hour6       -1554.94      60.59 -25.663  < 2e-16 ***
## hour7       -1618.23      60.59 -26.708  < 2e-16 ***
## hour8       -1624.51      60.59 -26.811  < 2e-16 ***
## hour9       -1635.51      60.59 -26.993  < 2e-16 ***
## hour10      -1586.87      60.59 -26.190  < 2e-16 ***
## hour11      -1523.80      60.59 -25.149  < 2e-16 ***
## hour12      -1278.23      60.59 -21.096  < 2e-16 ***
## hour13       -912.65      60.59 -15.063  < 2e-16 ***
## hour14       -670.80      60.59 -11.071  < 2e-16 ***
## hour15       -447.15      60.59  -7.380 1.51e-12 ***
## hour16       -319.23      60.59  -5.269 2.60e-07 ***
## hour17       -226.94      60.59  -3.745 0.000215 ***
## hour18       -128.51      60.59  -2.121 0.034729 *  
## hour19         69.92      60.59   1.154 0.249430    
## hour20        273.85      60.59   4.520 8.87e-06 ***
## hour21        258.74      61.80   4.187 3.71e-05 ***
## hour22        231.43      61.80   3.745 0.000216 ***
## hour23        174.35      60.59   2.877 0.004292 ** 
## dow2          -25.83      33.27  -0.776 0.438103    
## dow3           68.80      33.11   2.078 0.038584 *  
## dow4          -63.81      33.27  -1.918 0.056060 .  
## dow5          -76.19      33.27  -2.290 0.022715 *  
## dow6           34.72      33.66   1.032 0.303062    
## dow7          153.08      33.27   4.601 6.18e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 26569.29)
## 
##     Null deviance: 169592222  on 334  degrees of freedom
## Residual deviance:   8103633  on 305  degrees of freedom
## AIC: 4394.1
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  1699.045     55.226  30.7651 < 2.2e-16 ***
## hour1        -324.083     75.953  -4.2669 1.982e-05 ***
## hour2        -692.654     69.359  -9.9866 < 2.2e-16 ***
## hour3        -968.083     63.640 -15.2118 < 2.2e-16 ***
## hour4       -1213.726     56.334 -21.5451 < 2.2e-16 ***
## hour5       -1430.940     52.302 -27.3592 < 2.2e-16 ***
## hour6       -1554.940     52.136 -29.8247 < 2.2e-16 ***
## hour7       -1618.226     52.974 -30.5477 < 2.2e-16 ***
## hour8       -1624.511     60.443 -26.8769 < 2.2e-16 ***
## hour9       -1635.511     59.088 -27.6791 < 2.2e-16 ***
## hour10      -1586.869     67.302 -23.5782 < 2.2e-16 ***
## hour11      -1523.797     56.450 -26.9939 < 2.2e-16 ***
## hour12      -1278.226     55.626 -22.9789 < 2.2e-16 ***
## hour13       -912.654     55.048 -16.5793 < 2.2e-16 ***
## hour14       -670.797     67.064 -10.0024 < 2.2e-16 ***
## hour15       -447.154     65.300  -6.8477 7.507e-12 ***
## hour16       -319.226     73.849  -4.3227 1.541e-05 ***
## hour17       -226.940     72.497  -3.1303 0.0017460 ** 
## hour18       -128.511     73.915  -1.7386 0.0820985 .  
## hour19         69.917     64.835   1.0784 0.2808615    
## hour20        273.846     83.848   3.2660 0.0010909 ** 
## hour21        258.739     70.179   3.6868 0.0002271 ***
## hour22        231.431     71.603   3.2322 0.0012286 ** 
## hour23        174.346     81.418   2.1414 0.0322441 *  
## dow2          -25.833     35.746  -0.7227 0.4698635    
## dow3           68.795     35.059   1.9623 0.0497327 *  
## dow4          -63.812     31.859  -2.0030 0.0451788 *  
## dow5          -76.187     30.940  -2.4625 0.0137987 *  
## dow6           34.723     37.754   0.9197 0.3577250    
## dow7          153.083     38.020   4.0264 5.663e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -446.84   -97.67   -10.71    66.80   771.80  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)
## (Intercept)  10.2933    24.2221   0.425    0.671
## lin1          0.1541     0.2486   0.620    0.536
## lin2         -0.1429     0.2509  -0.570    0.569
## postchange    4.3103    34.1534   0.126    0.900
## 
## (Dispersion parameter for gaussian family taken to be 24422.43)
## 
##     Null deviance: 8103633  on 334  degrees of freedom
## Residual deviance: 8083823  on 331  degrees of freedom
## AIC: 4341.3
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##             Estimate Std. Error z value Pr(>|z|)
## (Intercept) 10.29335   26.85334  0.3833   0.7015
## lin1         0.15407    0.26569  0.5799   0.5620
## lin2        -0.14291    0.23851 -0.5992   0.5490
## postchange   4.31033   33.68727  0.1280   0.8982

## [1] "25 Regression for xbxcall model2c on 2016-11-01/purple/30-50%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -603.41   -78.85   -12.28    77.59   665.03  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1617.37      53.74  30.097  < 2e-16 ***
## hour1        -309.97      68.52  -4.524 8.70e-06 ***
## hour2        -655.75      68.52  -9.570  < 2e-16 ***
## hour3        -944.61      68.52 -13.786  < 2e-16 ***
## hour4       -1146.11      68.52 -16.726  < 2e-16 ***
## hour5       -1338.54      68.52 -19.535  < 2e-16 ***
## hour6       -1463.47      68.52 -21.358  < 2e-16 ***
## hour7       -1573.25      68.52 -22.960  < 2e-16 ***
## hour8       -1605.90      68.52 -23.436  < 2e-16 ***
## hour9       -1606.83      68.52 -23.450  < 2e-16 ***
## hour10      -1572.47      68.52 -22.949  < 2e-16 ***
## hour11      -1491.04      68.52 -21.760  < 2e-16 ***
## hour12      -1256.40      68.52 -18.336  < 2e-16 ***
## hour13       -930.33      68.52 -13.577  < 2e-16 ***
## hour14       -666.75      68.52  -9.731  < 2e-16 ***
## hour15       -473.90      68.52  -6.916 2.70e-11 ***
## hour16       -368.47      68.52  -5.377 1.50e-07 ***
## hour17       -294.33      68.52  -4.295 2.34e-05 ***
## hour18       -250.18      68.52  -3.651 0.000307 ***
## hour19        -81.18      68.52  -1.185 0.237021    
## hour20         47.25      68.52   0.690 0.491029    
## hour21        184.82      68.52   2.697 0.007379 ** 
## hour22        189.17      68.52   2.761 0.006113 ** 
## hour23         83.96      68.52   1.225 0.221399    
## dow2          -67.92      37.63  -1.805 0.072061 .  
## dow3           70.07      37.45   1.871 0.062260 .  
## dow4           54.52      37.63   1.449 0.148373    
## dow5           14.31      37.63   0.380 0.703934    
## dow6          122.65      37.63   3.259 0.001242 ** 
## dow7           85.08      37.63   2.261 0.024448 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 33980.55)
## 
##     Null deviance: 147749277  on 336  degrees of freedom
## Residual deviance:  10432028  on 307  degrees of freedom
## AIC: 4503
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  1617.366     65.745  24.6004 < 2.2e-16 ***
## hour1        -309.969     92.432  -3.3535 0.0007980 ***
## hour2        -655.754     82.065  -7.9906 1.342e-15 ***
## hour3        -944.611     76.353 -12.3717 < 2.2e-16 ***
## hour4       -1146.111     71.027 -16.1362 < 2.2e-16 ***
## hour5       -1338.540     71.339 -18.7631 < 2.2e-16 ***
## hour6       -1463.469     79.871 -18.3230 < 2.2e-16 ***
## hour7       -1573.254     58.732 -26.7871 < 2.2e-16 ***
## hour8       -1605.897     59.548 -26.9680 < 2.2e-16 ***
## hour9       -1606.826     58.980 -27.2437 < 2.2e-16 ***
## hour10      -1572.469     59.048 -26.6302 < 2.2e-16 ***
## hour11      -1491.040     58.968 -25.2858 < 2.2e-16 ***
## hour12      -1256.397     61.217 -20.5236 < 2.2e-16 ***
## hour13       -930.326     62.465 -14.8936 < 2.2e-16 ***
## hour14       -666.754     74.565  -8.9419 < 2.2e-16 ***
## hour15       -473.897     83.108  -5.7022 1.183e-08 ***
## hour16       -368.469     76.189  -4.8362 1.323e-06 ***
## hour17       -294.326     73.880  -3.9838 6.781e-05 ***
## hour18       -250.183     72.224  -3.4640 0.0005323 ***
## hour19        -81.183     70.410  -1.1530 0.2489081    
## hour20         47.246     78.749   0.6000 0.5485386    
## hour21        184.817     88.391   2.0909 0.0365357 *  
## hour22        189.174     95.900   1.9726 0.0485381 *  
## hour23         83.960    101.589   0.8265 0.4085405    
## dow2          -67.917     38.586  -1.7601 0.0783856 .  
## dow3           70.074     40.499   1.7303 0.0835838 .  
## dow4           54.521     33.525   1.6263 0.1038924    
## dow5           14.313     35.425   0.4040 0.6861992    
## dow6          122.646     41.690   2.9419 0.0032625 ** 
## dow7           85.083     48.948   1.7382 0.0821717 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -493.63   -97.00   -13.55    78.44   652.48  
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -162.8244    23.6198  -6.894 2.74e-11 ***
## lin1          -1.0825     0.2424  -4.465 1.10e-05 ***
## lin2           1.1482     0.2403   4.778 2.65e-06 ***
## postchange   137.3053    33.2070   4.135 4.50e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 23223.09)
## 
##     Null deviance: 10432028  on 336  degrees of freedom
## Residual deviance:  7733288  on 333  degrees of freedom
## AIC: 4350.2
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##               Estimate Std. Error z value  Pr(>|z|)    
## (Intercept) -162.82440   28.89766 -5.6345 1.755e-08 ***
## lin1          -1.08254    0.28506 -3.7976 0.0001461 ***
## lin2           1.14818    0.22462  5.1116 3.194e-07 ***
## postchange   137.30527   34.96035  3.9275 8.585e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "26 Regression for xbxcall model2c on 2016-12-15/purple/50-90%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -578.75   -98.27   -22.55   100.40  1343.75  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1801.24      56.85  31.682  < 2e-16 ***
## hour1        -181.83      72.49  -2.508  0.01265 *  
## hour2        -445.12      72.49  -6.140 2.55e-09 ***
## hour3        -696.83      72.49  -9.612  < 2e-16 ***
## hour4       -1037.69      72.49 -14.314  < 2e-16 ***
## hour5       -1295.62      72.49 -17.872  < 2e-16 ***
## hour6       -1497.97      72.49 -20.663  < 2e-16 ***
## hour7       -1654.90      72.49 -22.828  < 2e-16 ***
## hour8       -1748.76      72.49 -24.123  < 2e-16 ***
## hour9       -1779.90      72.49 -24.552  < 2e-16 ***
## hour10      -1785.97      72.49 -24.636  < 2e-16 ***
## hour11      -1757.12      72.49 -24.238  < 2e-16 ***
## hour12      -1656.76      72.49 -22.854  < 2e-16 ***
## hour13      -1392.33      72.49 -19.206  < 2e-16 ***
## hour14      -1015.62      72.49 -14.010  < 2e-16 ***
## hour15       -725.33      72.49 -10.005  < 2e-16 ***
## hour16       -457.62      72.49  -6.312 9.60e-10 ***
## hour17       -374.19      72.49  -5.162 4.41e-07 ***
## hour18       -307.26      72.49  -4.238 2.98e-05 ***
## hour19       -234.83      72.49  -3.239  0.00133 ** 
## hour20        -73.62      72.49  -1.015  0.31067    
## hour21         83.81      72.49   1.156  0.24853    
## hour22        183.81      72.49   2.536  0.01172 *  
## hour23        121.38      72.49   1.674  0.09507 .  
## dow2           63.62      39.81   1.598  0.11102    
## dow3          119.00      39.81   2.989  0.00302 ** 
## dow4          115.96      39.81   2.913  0.00384 ** 
## dow5         -103.88      39.62  -2.622  0.00918 ** 
## dow6          -66.31      39.81  -1.666  0.09678 .  
## dow7          118.25      39.81   2.970  0.00321 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 38034.76)
## 
##     Null deviance: 176190777  on 336  degrees of freedom
## Residual deviance:  11676672  on 307  degrees of freedom
## AIC: 4541
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  1801.240     63.677  28.2870 < 2.2e-16 ***
## hour1        -181.832     84.418  -2.1539 0.0312444 *  
## hour2        -445.117     86.001  -5.1757 2.271e-07 ***
## hour3        -696.832     85.981  -8.1045 5.296e-16 ***
## hour4       -1037.689     76.315 -13.5974 < 2.2e-16 ***
## hour5       -1295.617     67.386 -19.2267 < 2.2e-16 ***
## hour6       -1497.974     74.270 -20.1692 < 2.2e-16 ***
## hour7       -1654.903     71.648 -23.0976 < 2.2e-16 ***
## hour8       -1748.760     65.077 -26.8721 < 2.2e-16 ***
## hour9       -1779.903     65.718 -27.0841 < 2.2e-16 ***
## hour10      -1785.974     65.205 -27.3903 < 2.2e-16 ***
## hour11      -1757.117     65.786 -26.7095 < 2.2e-16 ***
## hour12      -1656.760     65.842 -25.1628 < 2.2e-16 ***
## hour13      -1392.332     66.844 -20.8296 < 2.2e-16 ***
## hour14      -1015.617     64.643 -15.7111 < 2.2e-16 ***
## hour15       -725.332     81.489  -8.9010 < 2.2e-16 ***
## hour16       -457.617    128.871  -3.5510 0.0003838 ***
## hour17       -374.189     88.187  -4.2432 2.204e-05 ***
## hour18       -307.260     81.620  -3.7645 0.0001669 ***
## hour19       -234.832     81.578  -2.8786 0.0039941 ** 
## hour20        -73.617     75.917  -0.9697 0.3321941    
## hour21         83.811     82.930   1.0106 0.3121967    
## hour22        183.811     94.087   1.9536 0.0507442 .  
## hour23        121.383     88.992   1.3640 0.1725749    
## dow2           63.625     49.467   1.2862 0.1983666    
## dow3          119.000     30.831   3.8598 0.0001135 ***
## dow4          115.958     32.085   3.6141 0.0003013 ***
## dow5         -103.882     31.040  -3.3467 0.0008177 ***
## dow6          -66.312     34.266  -1.9352 0.0529615 .  
## dow7          118.250     36.792   3.2140 0.0013091 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -541.05  -107.40   -11.41    85.29  1373.12  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)  
## (Intercept)   3.9691    27.7824   0.143   0.8865  
## lin1          0.5953     0.2852   2.088   0.0376 *
## lin2         -0.5424     0.2826  -1.919   0.0558 .
## postchange   87.6500    39.0592   2.244   0.0255 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 32129.66)
## 
##     Null deviance: 11676672  on 336  degrees of freedom
## Residual deviance: 10699177  on 333  degrees of freedom
## AIC: 4459.6
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##             Estimate Std. Error z value Pr(>|z|)   
## (Intercept)  3.96907   32.44882  0.1223 0.902647   
## lin1         0.59528    0.27311  2.1797 0.029283 * 
## lin2        -0.54239    0.20134 -2.6939 0.007062 **
## postchange  87.65001   38.68570  2.2657 0.023470 * 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "27 Regression for xbxcall model2c on 2016-11-17/blue/content change"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -964.12   -67.08     3.09    77.22   416.74  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  1648.86      46.15  35.730  < 2e-16 ***
## hour1        -143.18      55.68  -2.571 0.010640 *  
## hour2        -430.42      55.68  -7.730 1.88e-13 ***
## hour3        -765.88      55.68 -13.754  < 2e-16 ***
## hour4       -1012.34      55.68 -18.180  < 2e-16 ***
## hour5       -1247.95      55.68 -22.411  < 2e-16 ***
## hour6       -1461.80      55.68 -26.251  < 2e-16 ***
## hour7       -1566.34      55.68 -28.129  < 2e-16 ***
## hour8       -1616.57      55.68 -29.031  < 2e-16 ***
## hour9       -1632.88      55.68 -29.324  < 2e-16 ***
## hour10      -1616.26      55.68 -29.025  < 2e-16 ***
## hour11      -1609.57      55.68 -28.905  < 2e-16 ***
## hour12      -1477.72      55.68 -26.537  < 2e-16 ***
## hour13      -1226.95      55.68 -22.034  < 2e-16 ***
## hour14       -846.95      55.68 -15.210  < 2e-16 ***
## hour15       -632.11      55.68 -11.352  < 2e-16 ***
## hour16       -557.73      54.68 -10.200  < 2e-16 ***
## hour17       -347.05      55.75  -6.226 1.72e-09 ***
## hour18       -303.05      55.75  -5.436 1.17e-07 ***
## hour19       -217.51      55.75  -3.902 0.000119 ***
## hour20       -126.90      55.75  -2.276 0.023568 *  
## hour21         34.03      55.75   0.610 0.542085    
## hour22        121.95      55.75   2.188 0.029511 *  
## hour23         79.03      55.75   1.418 0.157390    
## dow2           16.71      33.11   0.505 0.614140    
## dow3           32.48      33.11   0.981 0.327425    
## dow4           84.59      33.11   2.555 0.011155 *  
## dow5          -51.26      33.00  -1.553 0.121443    
## dow6           25.86      33.11   0.781 0.435493    
## dow7          150.75      34.46   4.375 1.71e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 20889.16)
## 
##     Null deviance: 128843284  on 313  degrees of freedom
## Residual deviance:   5932522  on 284  degrees of freedom
## AIC: 4044.9
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error  z value  Pr(>|z|)    
## (Intercept)  1648.856     52.043  31.6824 < 2.2e-16 ***
## hour1        -143.185     64.379  -2.2241  0.026143 *  
## hour2        -430.416     59.602  -7.2215 5.144e-13 ***
## hour3        -765.877     53.600 -14.2888 < 2.2e-16 ***
## hour4       -1012.339     53.336 -18.9804 < 2.2e-16 ***
## hour5       -1247.954     45.777 -27.2617 < 2.2e-16 ***
## hour6       -1461.800     38.304 -38.1633 < 2.2e-16 ***
## hour7       -1566.339     39.090 -40.0700 < 2.2e-16 ***
## hour8       -1616.569     40.297 -40.1160 < 2.2e-16 ***
## hour9       -1632.877     39.179 -41.6776 < 2.2e-16 ***
## hour10      -1616.262     43.574 -37.0923 < 2.2e-16 ***
## hour11      -1609.569     40.802 -39.4483 < 2.2e-16 ***
## hour12      -1477.723     41.778 -35.3710 < 2.2e-16 ***
## hour13      -1226.954     42.427 -28.9192 < 2.2e-16 ***
## hour14       -846.954     40.772 -20.7727 < 2.2e-16 ***
## hour15       -632.108     53.830 -11.7426 < 2.2e-16 ***
## hour16       -557.733     90.878  -6.1371 8.402e-10 ***
## hour17       -347.050     56.737  -6.1168 9.547e-10 ***
## hour18       -303.050     65.342  -4.6379 3.519e-06 ***
## hour19       -217.512     62.304  -3.4911  0.000481 ***
## hour20       -126.896     52.855  -2.4008  0.016358 *  
## hour21         34.027     55.764   0.6102  0.541734    
## hour22        121.950     57.278   2.1291  0.033248 *  
## hour23         79.027     63.229   1.2499  0.211354    
## dow2           16.713     47.678   0.3505  0.725932    
## dow3           32.484     42.523   0.7639  0.444924    
## dow4           84.588     45.953   1.8407  0.065660 .  
## dow5          -51.258     42.477  -1.2067  0.227540    
## dow6           25.859     42.052   0.6149  0.538600    
## dow7          150.749     47.431   3.1783  0.001482 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -956.99   -68.83     2.16    71.34   462.23  
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)   
## (Intercept) -60.9644    22.4769  -2.712  0.00705 **
## lin1         -0.6729     0.2671  -2.519  0.01227 * 
## lin2          0.6325     0.2123   2.980  0.00311 **
## postchange   17.9969    30.5029   0.590  0.55562   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 18125.15)
## 
##     Null deviance: 5932522  on 313  degrees of freedom
## Residual deviance: 5618797  on 310  degrees of freedom
## AIC: 3975.9
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##              Estimate Std. Error z value  Pr(>|z|)    
## (Intercept) -60.96438   21.72018 -2.8068 0.0050035 ** 
## lin1         -0.67288    0.21189 -3.1756 0.0014953 ** 
## lin2          0.63250    0.17284  3.6594 0.0002528 ***
## postchange   17.99692   27.18816  0.6619 0.5080099    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "28 Regression for rate_xbxcall model2c on 2016-10-11/purple/0-10%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0053966  -0.0016096  -0.0000894   0.0010578   0.0146357  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.0198963  0.0007224  27.543  < 2e-16 ***
## hour1       -0.0026124  0.0009168  -2.849 0.004687 ** 
## hour2       -0.0063722  0.0009168  -6.950 2.32e-11 ***
## hour3       -0.0076267  0.0009168  -8.319 3.26e-15 ***
## hour4       -0.0106353  0.0009168 -11.600  < 2e-16 ***
## hour5       -0.0136972  0.0009168 -14.940  < 2e-16 ***
## hour6       -0.0153908  0.0009168 -16.787  < 2e-16 ***
## hour7       -0.0160953  0.0009168 -17.556  < 2e-16 ***
## hour8       -0.0162440  0.0009351 -17.370  < 2e-16 ***
## hour9       -0.0164751  0.0009351 -17.618  < 2e-16 ***
## hour10      -0.0160166  0.0009351 -17.127  < 2e-16 ***
## hour11      -0.0156722  0.0009351 -16.759  < 2e-16 ***
## hour12      -0.0145563  0.0009168 -15.877  < 2e-16 ***
## hour13      -0.0125659  0.0009168 -13.706  < 2e-16 ***
## hour14      -0.0114745  0.0009560 -12.003  < 2e-16 ***
## hour15      -0.0103074  0.0009560 -10.782  < 2e-16 ***
## hour16      -0.0082573  0.0009168  -9.007  < 2e-16 ***
## hour17      -0.0079140  0.0009168  -8.632 3.73e-16 ***
## hour18      -0.0074314  0.0009168  -8.106 1.38e-14 ***
## hour19      -0.0057353  0.0009168  -6.256 1.38e-09 ***
## hour20      -0.0036583  0.0009351  -3.912 0.000114 ***
## hour21      -0.0022194  0.0009351  -2.373 0.018264 *  
## hour22       0.0011921  0.0009168   1.300 0.194525    
## hour23       0.0014482  0.0009168   1.580 0.115269    
## dow2        -0.0043895  0.0005093  -8.618 4.11e-16 ***
## dow3        -0.0037254  0.0005069  -7.349 1.94e-12 ***
## dow4        -0.0039083  0.0005093  -7.673 2.42e-13 ***
## dow5        -0.0043553  0.0005093  -8.551 6.56e-16 ***
## dow6        -0.0032073  0.0005282  -6.072 3.85e-09 ***
## dow7         0.0013028  0.0005143   2.533 0.011815 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 6.083323e-06)
## 
##     Null deviance: 0.0143119  on 326  degrees of freedom
## Residual deviance: 0.0018067  on 297  degrees of freedom
## AIC: -2968.7
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  0.01989626  0.00079396  25.0596 < 2.2e-16 ***
## hour1       -0.00261237  0.00106705  -2.4482 0.0143563 *  
## hour2       -0.00637220  0.00073681  -8.6484 < 2.2e-16 ***
## hour3       -0.00762666  0.00135245  -5.6391 1.709e-08 ***
## hour4       -0.01063527  0.00090498 -11.7520 < 2.2e-16 ***
## hour5       -0.01369717  0.00074798 -18.3123 < 2.2e-16 ***
## hour6       -0.01539081  0.00079483 -19.3636 < 2.2e-16 ***
## hour7       -0.01609532  0.00083067 -19.3764 < 2.2e-16 ***
## hour8       -0.01624398  0.00098669 -16.4631 < 2.2e-16 ***
## hour9       -0.01647511  0.00101429 -16.2430 < 2.2e-16 ***
## hour10      -0.01601665  0.00109228 -14.6635 < 2.2e-16 ***
## hour11      -0.01567219  0.00087339 -17.9442 < 2.2e-16 ***
## hour12      -0.01455628  0.00073653 -19.7633 < 2.2e-16 ***
## hour13      -0.01256593  0.00068813 -18.2609 < 2.2e-16 ***
## hour14      -0.01147449  0.00086232 -13.3065 < 2.2e-16 ***
## hour15      -0.01030735  0.00088401 -11.6598 < 2.2e-16 ***
## hour16      -0.00825734  0.00103205  -8.0009 1.235e-15 ***
## hour17      -0.00791400  0.00095148  -8.3176 < 2.2e-16 ***
## hour18      -0.00743140  0.00107811  -6.8930 5.463e-12 ***
## hour19      -0.00573534  0.00097647  -5.8735 4.266e-09 ***
## hour20      -0.00365832  0.00099188  -3.6883 0.0002258 ***
## hour21      -0.00221943  0.00087844  -2.5266 0.0115185 *  
## hour22       0.00119209  0.00094539   1.2610 0.2073257    
## hour23       0.00144818  0.00099860   1.4502 0.1469987    
## dow2        -0.00438947  0.00044544  -9.8543 < 2.2e-16 ***
## dow3        -0.00372537  0.00050931  -7.3145 2.583e-13 ***
## dow4        -0.00390827  0.00052515  -7.4422 9.905e-14 ***
## dow5        -0.00435532  0.00060737  -7.1708 7.455e-13 ***
## dow6        -0.00320727  0.00049879  -6.4301 1.275e-10 ***
## dow7         0.00130283  0.00063525   2.0509 0.0402774 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0055921  -0.0015300  -0.0001391   0.0009793   0.0144543  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)
## (Intercept)  2.817e-04  3.767e-04   0.748    0.455
## lin1         9.376e-07  4.110e-06   0.228    0.820
## lin2         2.768e-07  3.715e-06   0.074    0.941
## postchange  -4.986e-04  5.217e-04  -0.956    0.340
## 
## (Dispersion parameter for gaussian family taken to be 5.552039e-06)
## 
##     Null deviance: 0.0018067  on 326  degrees of freedom
## Residual deviance: 0.0017933  on 323  degrees of freedom
## AIC: -3023.2
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)
## (Intercept)  2.8169e-04  3.0597e-04  0.9207   0.3572
## lin1         9.3762e-07  3.4811e-06  0.2693   0.7877
## lin2         2.7678e-07  3.4521e-06  0.0802   0.9361
## postchange  -4.9861e-04  4.4158e-04 -1.1291   0.2588

## [1] "29 Regression for rate_xbxcall model2c on 2016-10-18/purple/10-30%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -0.03951  -0.00171   0.00165   0.00362   0.34959  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)  
## (Intercept)  0.0178100  0.0082458   2.160   0.0316 *
## hour1       -0.0030695  0.0105131  -0.292   0.7705  
## hour2       -0.0064113  0.0105131  -0.610   0.5424  
## hour3       -0.0088660  0.0105131  -0.843   0.3997  
## hour4       -0.0113337  0.0105131  -1.078   0.2819  
## hour5       -0.0138939  0.0105131  -1.322   0.1873  
## hour6       -0.0155715  0.0105131  -1.481   0.1396  
## hour7       -0.0164120  0.0105131  -1.561   0.1195  
## hour8       -0.0165172  0.0105131  -1.571   0.1172  
## hour9       -0.0166805  0.0105131  -1.587   0.1136  
## hour10      -0.0162493  0.0105131  -1.546   0.1232  
## hour11      -0.0160193  0.0105131  -1.524   0.1286  
## hour12      -0.0147646  0.0105131  -1.404   0.1612  
## hour13      -0.0129473  0.0105131  -1.232   0.2191  
## hour14      -0.0114366  0.0105131  -1.088   0.2775  
## hour15      -0.0102999  0.0105131  -0.980   0.3280  
## hour16      -0.0090027  0.0105131  -0.856   0.3925  
## hour17      -0.0083274  0.0105131  -0.792   0.4289  
## hour18      -0.0038448  0.0105131  -0.366   0.7148  
## hour19       0.0182298  0.0105131   1.734   0.0839 .
## hour20       0.0243723  0.0105131   2.318   0.0211 *
## hour21      -0.0007475  0.0107230  -0.070   0.9445  
## hour22       0.0013488  0.0107230   0.126   0.9000  
## hour23       0.0042471  0.0105131   0.404   0.6865  
## dow2        -0.0049582  0.0057732  -0.859   0.3911  
## dow3        -0.0033451  0.0057455  -0.582   0.5609  
## dow4        -0.0044536  0.0057732  -0.771   0.4411  
## dow5        -0.0049375  0.0057732  -0.855   0.3931  
## dow6         0.0143698  0.0058401   2.461   0.0144 *
## dow7         0.0016000  0.0057732   0.277   0.7819  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 0.0007999039)
## 
##     Null deviance: 0.29392  on 334  degrees of freedom
## Residual deviance: 0.24397  on 305  degrees of freedom
## AIC: -1407.6
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value  Pr(>|z|)    
## (Intercept)  0.01781000  0.00255193  6.9790 2.972e-12 ***
## hour1       -0.00306950  0.00272012 -1.1284 0.2591335    
## hour2       -0.00641134  0.00264097 -2.4276 0.0151971 *  
## hour3       -0.00886602  0.00260649 -3.4015 0.0006701 ***
## hour4       -0.01133368  0.00254824 -4.4476 8.682e-06 ***
## hour5       -0.01389387  0.00256510 -5.4165 6.078e-08 ***
## hour6       -0.01557155  0.00256184 -6.0783 1.215e-09 ***
## hour7       -0.01641196  0.00258665 -6.3449 2.226e-10 ***
## hour8       -0.01651725  0.00264979 -6.2334 4.564e-10 ***
## hour9       -0.01668052  0.00263849 -6.3220 2.582e-10 ***
## hour10      -0.01624929  0.00266207 -6.1040 1.034e-09 ***
## hour11      -0.01601926  0.00260503 -6.1494 7.780e-10 ***
## hour12      -0.01476460  0.00259079 -5.6989 1.206e-08 ***
## hour13      -0.01294729  0.00261153 -4.9577 7.132e-07 ***
## hour14      -0.01143658  0.00270842 -4.2226 2.415e-05 ***
## hour15      -0.01029994  0.00273076 -3.7718 0.0001621 ***
## hour16      -0.00900274  0.00276335 -3.2579 0.0011224 ** 
## hour17      -0.00832736  0.00271726 -3.0646 0.0021795 ** 
## hour18      -0.00384476  0.00357603 -1.0751 0.2823088    
## hour19       0.01822977  0.02334965  0.7807 0.4349613    
## hour20       0.02437226  0.02717636  0.8968 0.3698158    
## hour21      -0.00074748  0.00250015 -0.2990 0.7649615    
## hour22       0.00134876  0.00244715  0.5512 0.5815277    
## hour23       0.00424714  0.00331043  1.2830 0.1995074    
## dow2        -0.00495818  0.00157841 -3.1412 0.0016823 ** 
## dow3        -0.00334512  0.00160845 -2.0797 0.0375508 *  
## dow4        -0.00445357  0.00160241 -2.7793 0.0054477 ** 
## dow5        -0.00493752  0.00162280 -3.0426 0.0023456 ** 
## dow6         0.01436979  0.01088511  1.3201 0.1867908    
## dow7         0.00159996  0.00147110  1.0876 0.2767744    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##      Min        1Q    Median        3Q       Max  
## -0.03670  -0.00332   0.00100   0.00434   0.34674  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)
## (Intercept) -3.421e-03  4.186e-03  -0.817    0.414
## lin1        -8.023e-06  4.296e-05  -0.187    0.852
## lin2         9.295e-06  4.335e-05   0.214    0.830
## postchange   5.410e-03  5.902e-03   0.917    0.360
## 
## (Dispersion parameter for gaussian family taken to be 0.000729231)
## 
##     Null deviance: 0.24397  on 334  degrees of freedom
## Residual deviance: 0.24138  on 331  degrees of freedom
## AIC: -1463.2
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)  
## (Intercept) -3.4212e-03  1.3919e-03 -2.4579  0.01398 *
## lin1        -8.0227e-06  1.3274e-05 -0.6044  0.54557  
## lin2         9.2948e-06  1.6591e-05  0.5602  0.57531  
## postchange   5.4095e-03  2.6877e-03  2.0127  0.04414 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "30 Regression for rate_xbxcall model2c on 2016-11-01/purple/30-50%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0061684  -0.0014791  -0.0001982   0.0010462   0.0079468  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.0198739  0.0006911  28.757  < 2e-16 ***
## hour1       -0.0029718  0.0008812  -3.372 0.000841 ***
## hour2       -0.0060178  0.0008812  -6.829 4.59e-11 ***
## hour3       -0.0086865  0.0008812  -9.857  < 2e-16 ***
## hour4       -0.0107781  0.0008812 -12.231  < 2e-16 ***
## hour5       -0.0131094  0.0008812 -14.876  < 2e-16 ***
## hour6       -0.0149713  0.0008812 -16.989  < 2e-16 ***
## hour7       -0.0163731  0.0008812 -18.580  < 2e-16 ***
## hour8       -0.0168128  0.0008812 -19.079  < 2e-16 ***
## hour9       -0.0168620  0.0008812 -19.135  < 2e-16 ***
## hour10      -0.0165549  0.0008812 -18.786  < 2e-16 ***
## hour11      -0.0160029  0.0008812 -18.160  < 2e-16 ***
## hour12      -0.0147806  0.0008812 -16.773  < 2e-16 ***
## hour13      -0.0131320  0.0008812 -14.902  < 2e-16 ***
## hour14      -0.0114931  0.0008812 -13.042  < 2e-16 ***
## hour15      -0.0103683  0.0008812 -11.766  < 2e-16 ***
## hour16      -0.0095339  0.0008812 -10.819  < 2e-16 ***
## hour17      -0.0088819  0.0008812 -10.079  < 2e-16 ***
## hour18      -0.0084916  0.0008812  -9.636  < 2e-16 ***
## hour19      -0.0070754  0.0008812  -8.029 2.11e-14 ***
## hour20      -0.0055380  0.0008812  -6.285 1.13e-09 ***
## hour21      -0.0028365  0.0008812  -3.219 0.001425 ** 
## hour22      -0.0006412  0.0008812  -0.728 0.467371    
## hour23       0.0001140  0.0008812   0.129 0.897161    
## dow2        -0.0046623  0.0004839  -9.635  < 2e-16 ***
## dow3        -0.0032814  0.0004816  -6.814 5.03e-11 ***
## dow4        -0.0032289  0.0004839  -6.672 1.17e-10 ***
## dow5        -0.0039765  0.0004839  -8.218 5.89e-15 ***
## dow6        -0.0025867  0.0004839  -5.345 1.76e-07 ***
## dow7         0.0009608  0.0004839   1.986 0.047970 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 5.620059e-06)
## 
##     Null deviance: 0.0131476  on 336  degrees of freedom
## Residual deviance: 0.0017254  on 307  degrees of freedom
## AIC: -3087.1
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  0.01987393  0.00085667  23.1990 < 2.2e-16 ***
## hour1       -0.00297183  0.00115783  -2.5667  0.010267 *  
## hour2       -0.00601783  0.00101762  -5.9136 3.346e-09 ***
## hour3       -0.00868650  0.00093132  -9.3271 < 2.2e-16 ***
## hour4       -0.01077807  0.00094984 -11.3472 < 2.2e-16 ***
## hour5       -0.01310939  0.00094842 -13.8224 < 2.2e-16 ***
## hour6       -0.01497131  0.00103441 -14.4733 < 2.2e-16 ***
## hour7       -0.01637314  0.00085331 -19.1877 < 2.2e-16 ***
## hour8       -0.01681279  0.00087187 -19.2837 < 2.2e-16 ***
## hour9       -0.01686200  0.00088052 -19.1500 < 2.2e-16 ***
## hour10      -0.01655492  0.00090004 -18.3936 < 2.2e-16 ***
## hour11      -0.01600286  0.00086389 -18.5241 < 2.2e-16 ***
## hour12      -0.01478056  0.00079266 -18.6467 < 2.2e-16 ***
## hour13      -0.01313197  0.00074998 -17.5097 < 2.2e-16 ***
## hour14      -0.01149312  0.00092187 -12.4671 < 2.2e-16 ***
## hour15      -0.01036831  0.00104154  -9.9548 < 2.2e-16 ***
## hour16      -0.00953391  0.00095184 -10.0162 < 2.2e-16 ***
## hour17      -0.00888189  0.00098320  -9.0337 < 2.2e-16 ***
## hour18      -0.00849161  0.00094587  -8.9775 < 2.2e-16 ***
## hour19      -0.00707539  0.00089167  -7.9350 2.105e-15 ***
## hour20      -0.00553805  0.00085181  -6.5015 7.951e-11 ***
## hour21      -0.00283649  0.00100929  -2.8104  0.004948 ** 
## hour22      -0.00064123  0.00103074  -0.6221  0.533870    
## hour23       0.00011399  0.00116125   0.0982  0.921803    
## dow2        -0.00466226  0.00054765  -8.5131 < 2.2e-16 ***
## dow3        -0.00328144  0.00055312  -5.9326 2.982e-09 ***
## dow4        -0.00322887  0.00049296  -6.5499 5.756e-11 ***
## dow5        -0.00397654  0.00051682  -7.6942 1.424e-14 ***
## dow6        -0.00258671  0.00056758  -4.5575 5.177e-06 ***
## dow7         0.00096084  0.00072181   1.3312  0.183135    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0045761  -0.0013318  -0.0002365   0.0010365   0.0077795  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.913e-03  3.134e-04  -6.105 2.86e-09 ***
## lin1        -1.284e-05  3.217e-06  -3.993 8.05e-05 ***
## lin2         1.307e-05  3.189e-06   4.099 5.22e-05 ***
## postchange   1.639e-03  4.407e-04   3.719 0.000235 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 4.089313e-06)
## 
##     Null deviance: 0.0017254  on 336  degrees of freedom
## Residual deviance: 0.0013617  on 333  degrees of freedom
## AIC: -3218.9
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value  Pr(>|z|)    
## (Intercept) -1.9134e-03  3.2754e-04 -5.8418 5.165e-09 ***
## lin1        -1.2844e-05  3.2503e-06 -3.9517 7.758e-05 ***
## lin2         1.3070e-05  3.0495e-06  4.2861 1.819e-05 ***
## postchange   1.6387e-03  3.8798e-04  4.2237 2.404e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "31 Regression for rate_xbxcall model2c on 2016-12-15/purple/50-90%"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0052973  -0.0016974  -0.0001152   0.0013397   0.0079829  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  2.194e-02  6.834e-04  32.097  < 2e-16 ***
## hour1       -2.110e-03  8.714e-04  -2.422   0.0160 *  
## hour2       -4.207e-03  8.714e-04  -4.827 2.18e-06 ***
## hour3       -6.544e-03  8.714e-04  -7.510 6.49e-13 ***
## hour4       -9.914e-03  8.714e-04 -11.377  < 2e-16 ***
## hour5       -1.255e-02  8.714e-04 -14.407  < 2e-16 ***
## hour6       -1.509e-02  8.714e-04 -17.314  < 2e-16 ***
## hour7       -1.707e-02  8.714e-04 -19.586  < 2e-16 ***
## hour8       -1.824e-02  8.714e-04 -20.929  < 2e-16 ***
## hour9       -1.869e-02  8.714e-04 -21.443  < 2e-16 ***
## hour10      -1.877e-02  8.714e-04 -21.535  < 2e-16 ***
## hour11      -1.852e-02  8.714e-04 -21.253  < 2e-16 ***
## hour12      -1.778e-02  8.714e-04 -20.399  < 2e-16 ***
## hour13      -1.642e-02  8.714e-04 -18.844  < 2e-16 ***
## hour14      -1.433e-02  8.714e-04 -16.447  < 2e-16 ***
## hour15      -1.253e-02  8.714e-04 -14.374  < 2e-16 ***
## hour16      -1.112e-02  8.714e-04 -12.758  < 2e-16 ***
## hour17      -9.814e-03  8.714e-04 -11.262  < 2e-16 ***
## hour18      -8.973e-03  8.714e-04 -10.297  < 2e-16 ***
## hour19      -8.341e-03  8.714e-04  -9.571  < 2e-16 ***
## hour20      -6.748e-03  8.714e-04  -7.743 1.42e-13 ***
## hour21      -4.677e-03  8.714e-04  -5.366 1.58e-07 ***
## hour22      -2.011e-03  8.714e-04  -2.308   0.0217 *  
## hour23       4.491e-05  8.714e-04   0.052   0.9589    
## dow2        -4.131e-03  4.785e-04  -8.633 3.30e-16 ***
## dow3        -2.698e-03  4.785e-04  -5.637 3.92e-08 ***
## dow4        -2.362e-03  4.785e-04  -4.937 1.31e-06 ***
## dow5        -5.016e-03  4.763e-04 -10.532  < 2e-16 ***
## dow6        -4.178e-03  4.785e-04  -8.731  < 2e-16 ***
## dow7         8.303e-04  4.785e-04   1.735   0.0837 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 5.49608e-06)
## 
##     Null deviance: 0.0158925  on 336  degrees of freedom
## Residual deviance: 0.0016873  on 307  degrees of freedom
## AIC: -3094.6
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  2.1936e-02  9.5717e-04  22.9180 < 2.2e-16 ***
## hour1       -2.1104e-03  1.0703e-03  -1.9717 0.0486385 *  
## hour2       -4.2069e-03  1.0976e-03  -3.8330 0.0001266 ***
## hour3       -6.5444e-03  1.1486e-03  -5.6975 1.216e-08 ***
## hour4       -9.9141e-03  1.0517e-03  -9.4271 < 2.2e-16 ***
## hour5       -1.2555e-02  9.4455e-04 -13.2920 < 2.2e-16 ***
## hour6       -1.5088e-02  1.0585e-03 -14.2551 < 2.2e-16 ***
## hour7       -1.7068e-02  1.0667e-03 -16.0011 < 2.2e-16 ***
## hour8       -1.8239e-02  9.9621e-04 -18.3081 < 2.2e-16 ***
## hour9       -1.8686e-02  1.0280e-03 -18.1774 < 2.2e-16 ***
## hour10      -1.8766e-02  1.0252e-03 -18.3058 < 2.2e-16 ***
## hour11      -1.8521e-02  1.0444e-03 -17.7332 < 2.2e-16 ***
## hour12      -1.7777e-02  9.9884e-04 -17.7973 < 2.2e-16 ***
## hour13      -1.6421e-02  9.5059e-04 -17.2747 < 2.2e-16 ***
## hour14      -1.4333e-02  9.1692e-04 -15.6312 < 2.2e-16 ***
## hour15      -1.2526e-02  1.0624e-03 -11.7900 < 2.2e-16 ***
## hour16      -1.1117e-02  1.0964e-03 -10.1399 < 2.2e-16 ***
## hour17      -9.8139e-03  1.1596e-03  -8.4631 < 2.2e-16 ***
## hour18      -8.9730e-03  1.1839e-03  -7.5794 3.472e-14 ***
## hour19      -8.3409e-03  1.1995e-03  -6.9536 3.561e-12 ***
## hour20      -6.7480e-03  1.1326e-03  -5.9578 2.557e-09 ***
## hour21      -4.6766e-03  1.1484e-03  -4.0721 4.659e-05 ***
## hour22      -2.0109e-03  1.1868e-03  -1.6944 0.0901868 .  
## hour23       4.4907e-05  1.1803e-03   0.0380 0.9696489    
## dow2        -4.1313e-03  5.0729e-04  -8.1437 3.832e-16 ***
## dow3        -2.6976e-03  4.7052e-04  -5.7333 9.849e-09 ***
## dow4        -2.3624e-03  4.6986e-04  -5.0278 4.961e-07 ***
## dow5        -5.0157e-03  4.7700e-04 -10.5151 < 2.2e-16 ***
## dow6        -4.1781e-03  4.9316e-04  -8.4721 < 2.2e-16 ***
## dow7         8.3027e-04  6.2440e-04   1.3297 0.1836131    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0043648  -0.0015705  -0.0002936   0.0010743   0.0071270  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -7.515e-04  3.209e-04  -2.342 0.019787 *  
## lin1         1.508e-06  3.294e-06   0.458 0.647342    
## lin2        -4.953e-07  3.265e-06  -0.152 0.879503    
## postchange   1.667e-03  4.512e-04   3.694 0.000258 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 4.286827e-06)
## 
##     Null deviance: 0.0016873  on 336  degrees of freedom
## Residual deviance: 0.0014275  on 333  degrees of freedom
## AIC: -3203
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value  Pr(>|z|)    
## (Intercept) -7.5147e-04  2.7367e-04 -2.7459  0.006034 ** 
## lin1         1.5082e-06  2.8289e-06  0.5331  0.593948    
## lin2        -4.9530e-07  2.5440e-06 -0.1947  0.845634    
## postchange   1.6668e-03  3.9956e-04  4.1716 3.025e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## [1] "32 Regression for rate_xbxcall model2c on 2016-11-17/blue/content change"
## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0081327  -0.0010820  -0.0000354   0.0010022   0.0166581  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.0214892  0.0006888  31.199  < 2e-16 ***
## hour1       -0.0027911  0.0008311  -3.358 0.000892 ***
## hour2       -0.0050140  0.0008311  -6.033 5.01e-09 ***
## hour3       -0.0082452  0.0008311  -9.920  < 2e-16 ***
## hour4       -0.0105464  0.0008311 -12.689  < 2e-16 ***
## hour5       -0.0130700  0.0008311 -15.726  < 2e-16 ***
## hour6       -0.0158277  0.0008311 -19.044  < 2e-16 ***
## hour7       -0.0171759  0.0008311 -20.666  < 2e-16 ***
## hour8       -0.0178586  0.0008311 -21.487  < 2e-16 ***
## hour9       -0.0180537  0.0008311 -21.722  < 2e-16 ***
## hour10      -0.0177976  0.0008311 -21.414  < 2e-16 ***
## hour11      -0.0179295  0.0008311 -21.572  < 2e-16 ***
## hour12      -0.0169602  0.0008311 -20.406  < 2e-16 ***
## hour13      -0.0157293  0.0008311 -18.925  < 2e-16 ***
## hour14      -0.0136625  0.0008311 -16.438  < 2e-16 ***
## hour15      -0.0123838  0.0008311 -14.900  < 2e-16 ***
## hour16      -0.0120646  0.0008161 -14.783  < 2e-16 ***
## hour17      -0.0100940  0.0008320 -12.132  < 2e-16 ***
## hour18      -0.0094978  0.0008320 -11.415  < 2e-16 ***
## hour19      -0.0088177  0.0008320 -10.598  < 2e-16 ***
## hour20      -0.0077850  0.0008320  -9.357  < 2e-16 ***
## hour21      -0.0057563  0.0008320  -6.918 3.03e-11 ***
## hour22      -0.0031857  0.0008320  -3.829 0.000158 ***
## hour23      -0.0004680  0.0008320  -0.562 0.574237    
## dow2        -0.0045357  0.0004942  -9.177  < 2e-16 ***
## dow3        -0.0036820  0.0004942  -7.450 1.13e-12 ***
## dow4        -0.0027765  0.0004942  -5.618 4.61e-08 ***
## dow5        -0.0047412  0.0004925  -9.627  < 2e-16 ***
## dow6        -0.0030225  0.0004942  -6.116 3.18e-09 ***
## dow7         0.0013770  0.0005143   2.677 0.007852 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 4.653619e-06)
## 
##     Null deviance: 0.0127004  on 313  degrees of freedom
## Residual deviance: 0.0013216  on 284  degrees of freedom
## AIC: -2933.7
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error  z value  Pr(>|z|)    
## (Intercept)  0.02148918  0.00140301  15.3165 < 2.2e-16 ***
## hour1       -0.00279110  0.00144507  -1.9315 0.0534259 .  
## hour2       -0.00501397  0.00142793  -3.5114 0.0004458 ***
## hour3       -0.00824518  0.00140762  -5.8575 4.698e-09 ***
## hour4       -0.01054638  0.00141648  -7.4455 9.658e-14 ***
## hour5       -0.01307004  0.00140670  -9.2912 < 2.2e-16 ***
## hour6       -0.01582774  0.00140435 -11.2705 < 2.2e-16 ***
## hour7       -0.01717587  0.00142833 -12.0251 < 2.2e-16 ***
## hour8       -0.01785857  0.00145854 -12.2441 < 2.2e-16 ***
## hour9       -0.01805372  0.00144534 -12.4910 < 2.2e-16 ***
## hour10      -0.01779757  0.00145353 -12.2444 < 2.2e-16 ***
## hour11      -0.01792948  0.00148614 -12.0645 < 2.2e-16 ***
## hour12      -0.01696015  0.00144816 -11.7115 < 2.2e-16 ***
## hour13      -0.01572926  0.00140367 -11.2058 < 2.2e-16 ***
## hour14      -0.01366255  0.00138166  -9.8885 < 2.2e-16 ***
## hour15      -0.01238375  0.00143614  -8.6229 < 2.2e-16 ***
## hour16      -0.01206457  0.00156355  -7.7161 1.199e-14 ***
## hour17      -0.01009398  0.00144610  -6.9801 2.949e-12 ***
## hour18      -0.00949779  0.00147189  -6.4528 1.098e-10 ***
## hour19      -0.00881772  0.00148433  -5.9405 2.841e-09 ***
## hour20      -0.00778501  0.00142937  -5.4465 5.138e-08 ***
## hour21      -0.00575629  0.00139931  -4.1137 3.894e-05 ***
## hour22      -0.00318566  0.00138152  -2.3059 0.0211156 *  
## hour23      -0.00046800  0.00179114  -0.2613 0.7938732    
## dow2        -0.00453571  0.00055853  -8.1208 4.633e-16 ***
## dow3        -0.00368200  0.00054098  -6.8062 1.002e-11 ***
## dow4        -0.00277645  0.00058087  -4.7798 1.755e-06 ***
## dow5        -0.00474119  0.00057250  -8.2815 < 2.2e-16 ***
## dow6        -0.00302249  0.00063696  -4.7452 2.083e-06 ***
## dow7         0.00137696  0.00077073   1.7866 0.0740089 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

## 
## Call:
## glm(formula = form, data = df)
## 
## Deviance Residuals: 
##        Min          1Q      Median          3Q         Max  
## -0.0081227  -0.0011831   0.0000357   0.0010280   0.0168225  
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)  
## (Intercept) -5.645e-04  3.411e-04  -1.655   0.0989 .
## lin1        -6.932e-06  4.054e-06  -1.710   0.0883 .
## lin2         5.962e-06  3.222e-06   1.850   0.0652 .
## postchange   1.140e-04  4.629e-04   0.246   0.8057  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for gaussian family taken to be 4.174845e-06)
## 
##     Null deviance: 0.0013216  on 313  degrees of freedom
## Residual deviance: 0.0012942  on 310  degrees of freedom
## AIC: -2992.3
## 
## Number of Fisher Scoring iterations: 2
## 
## 
## z test of coefficients:
## 
##                Estimate  Std. Error z value Pr(>|z|)  
## (Intercept) -5.6455e-04  2.3515e-04 -2.4008  0.01636 *
## lin1        -6.9321e-06  2.6939e-06 -2.5733  0.01007 *
## lin2         5.9617e-06  2.8448e-06  2.0957  0.03611 *
## postchange   1.1396e-04  4.4457e-04  0.2563  0.79770  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Results

Step 1 Results

idx title model formula chgDate chgVal aic deviance
1 SMC - Hourly Chats model2c winchib ~ hour + dow 2016-08-17 0_30 3574.605 6.635546e+05
2 SMC - Hourly Chats model2c winchib ~ hour + dow 2016-09-1 30_50 3491.759 5.189343e+05
3 SMC - Hourly Chats model2c winchib ~ hour + dow 2016-09-07 50_100 3452.170 4.614167e+05
4 SMC - Hourly Chat Rate per Session model2c rate_winchib ~ hour + dow 2016-08-17 0_30 -4343.877 4.140000e-05
5 SMC - Hourly Chat Rate per Session model2c rate_winchib ~ hour + dow 2016-09-1 30_50 -4296.916 4.760000e-05
6 SMC - Hourly Chat Rate per Session model2c rate_winchib ~ hour + dow 2016-09-07 50_100 -4283.854 4.950000e-05
7 SMC - Hourly Chats model2c wincall ~ hour + dow 2016-08-17 0_30 3793.924 1.272079e+06
8 SMC - Hourly Chats model2c wincall ~ hour + dow 2016-09-1 30_50 3861.978 1.556739e+06
9 SMC - Hourly Chats model2c wincall ~ hour + dow 2016-09-07 50_100 3895.735 1.720755e+06
10 SMC - Hourly Chat Rate per Session model2c rate_wincall ~ hour + dow 2016-08-17 0_30 -4378.675 3.740000e-05
11 SMC - Hourly Chat Rate per Session model2c rate_wincall ~ hour + dow 2016-09-1 30_50 -3993.703 1.171000e-04
12 SMC - Hourly Chat Rate per Session model2c rate_wincall ~ hour + dow 2016-09-07 50_100 -3945.263 1.352000e-04
13 Xbox - Hourly Chats model2c xbxchib ~ hour + dow 2016-10-11 0_10 3329.835 4.189309e+05
14 Xbox - Hourly Chats model2c xbxchib ~ hour + dow 2016-10-18 10_30 3255.501 2.707879e+05
15 Xbox - Hourly Chats model2c xbxchib ~ hour + dow 2016-11-01 30_50 3556.853 6.295051e+05
16 Xbox - Hourly Chats model2c xbxchib ~ hour + dow 2016-12-15 50_90 3478.329 4.986604e+05
17 Xbox - Hourly Chats model2c xbxchib ~ hour + dow 2016-11-17 content change 3298.600 5.508178e+05
18 Xbox - Hourly Chat Rate per Session model2c rate_xbxchib ~ hour + dow 2016-10-11 0_10 -4085.187 5.940000e-05
19 Xbox - Hourly Chat Rate per Session model2c rate_xbxchib ~ hour + dow 2016-10-18 10_30 -2775.264 4.114700e-03
20 Xbox - Hourly Chat Rate per Session model2c rate_xbxchib ~ hour + dow 2016-11-01 30_50 -4026.424 1.063000e-04
21 Xbox - Hourly Chat Rate per Session model2c rate_xbxchib ~ hour + dow 2016-12-15 50_90 -4106.839 8.370000e-05
22 Xbox - Hourly Chat Rate per Session model2c rate_xbxchib ~ hour + dow 2016-11-17 content change -3723.680 1.068000e-04
23 Xbox - Hourly Chats model2c xbxcall ~ hour + dow 2016-10-11 0_10 4450.684 1.290512e+07
24 Xbox - Hourly Chats model2c xbxcall ~ hour + dow 2016-10-18 10_30 4394.076 8.103633e+06
25 Xbox - Hourly Chats model2c xbxcall ~ hour + dow 2016-11-01 30_50 4503.048 1.043203e+07
26 Xbox - Hourly Chats model2c xbxcall ~ hour + dow 2016-12-15 50_90 4541.033 1.167667e+07
27 Xbox - Hourly Chats model2c xbxcall ~ hour + dow 2016-11-17 content change 4044.915 5.932522e+06
28 Xbox - Hourly Chat Rate per Session model2c rate_xbxcall ~ hour + dow 2016-10-11 0_10 -2968.738 1.806700e-03
29 Xbox - Hourly Chat Rate per Session model2c rate_xbxcall ~ hour + dow 2016-10-18 10_30 -1407.632 2.439707e-01
30 Xbox - Hourly Chat Rate per Session model2c rate_xbxcall ~ hour + dow 2016-11-01 30_50 -3087.105 1.725400e-03
31 Xbox - Hourly Chat Rate per Session model2c rate_xbxcall ~ hour + dow 2016-12-15 50_90 -3094.623 1.687300e-03
32 Xbox - Hourly Chat Rate per Session model2c rate_xbxcall ~ hour + dow 2016-11-17 content change -2933.688 1.321600e-03

Step 2 Results

idx title model formula chgDate chgVal chgCoef.fit chgCoef.smry chgStd.smry chgZval.cfit chgPval.cfit s
1 SMC - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-08-17 0_30 3.6099102 3.6099102 9.6593264 0.3772374 0.7059972 |
2 SMC - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-09-1 30_50 2.5569213 2.5569213 8.5573647 0.3309858 0.7406552 |
3 SMC - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-09-07 50_100 -5.7870962 -5.7870962 8.0803093 -0.6198333 0.5353675 |
4 SMC - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-08-17 0_30 0.0001094 0.0001094 0.0000760 1.4879623 0.1367608 |
5 SMC - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-09-1 30_50 0.0001432 0.0001432 0.0000804 2.1996495 0.0278318 | *
6 SMC - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-09-07 50_100 -0.0001240 -0.0001240 0.0000812 -1.4806961 0.1386876 |
7 SMC - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-08-17 0_30 -16.1807861 -16.1807861 13.3503031 -1.4933531 0.1353448 |
8 SMC - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-09-1 30_50 -34.3358166 -34.3358166 14.3735602 -2.8179245 0.0048335 | **
9 SMC - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-09-07 50_100 12.4489877 12.4489877 15.2886221 0.7676351 0.4427040 |
10 SMC - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-08-17 0_30 -0.0000836 -0.0000836 0.0000720 -1.3555447 0.1752441 |
11 SMC - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-09-1 30_50 -0.0000482 -0.0000482 0.0001290 -0.5709237 0.5680514 |
12 SMC - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-09-07 50_100 -0.0000784 -0.0000784 0.0001386 -0.3417116 0.7325680 |
13 Xbox - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-10-11 0_10 -19.2428551 -19.2428551 7.3410513 -2.8117013 0.0049280 | **
14 Xbox - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-10-18 10_30 -3.4163190 -3.4163190 6.1462413 -0.6065983 0.5441176 |
15 Xbox - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-11-01 30_50 33.9929346 33.9929346 8.0466906 4.8027829 0.0000016 | ***
16 Xbox - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-12-15 50_90 31.2028562 31.2028562 7.6219032 4.1857564 0.0000284 | ***
17 Xbox - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-11-17 content change 12.4953521 12.4953521 9.4422984 1.6555991 0.0978030 | .
18 Xbox - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-10-11 0_10 -0.0002296 -0.0002296 0.0000906 -2.7460017 0.0060326 | **
19 Xbox - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-10-18 10_30 0.0007107 0.0007107 0.0007664 2.0194630 0.0434391 | *
20 Xbox - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-11-01 30_50 0.0004420 0.0004420 0.0001081 5.1210113 0.0000003 | ***
21 Xbox - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-12-15 50_90 0.0004617 0.0004617 0.0000955 5.4625347 0.0000000 | ***
22 Xbox - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-11-17 content change 0.0001385 0.0001385 0.0001323 1.4122389 0.1578796 |
23 Xbox - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-10-11 0_10 -15.2301854 -15.2301854 44.2195628 -0.4057033 0.6849606 |
24 Xbox - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-10-18 10_30 4.3103307 4.3103307 34.1533979 0.1279513 0.8981875 |
25 Xbox - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-11-01 30_50 137.3052654 137.3052654 33.2070178 3.9274570 0.0000858 | ***
26 Xbox - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-12-15 50_90 87.6500100 87.6500100 39.0591661 2.2656952 0.0234700 | *
27 Xbox - Hourly Chats model2c step1residuals ~ lin1 + lin2 + postchange 2016-11-17 content change 17.9969247 17.9969247 30.5028676 0.6619397 0.5080099 |
28 Xbox - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-10-11 0_10 -0.0004986 -0.0004986 0.0005217 -1.1291371 0.2588400 |
29 Xbox - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-10-18 10_30 0.0054095 0.0054095 0.0059016 2.0127194 0.0441442 | *
30 Xbox - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-11-01 30_50 0.0016387 0.0016387 0.0004407 4.2236696 0.0000240 | ***
31 Xbox - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-12-15 50_90 0.0016668 0.0016668 0.0004512 4.1715945 0.0000302 | ***
32 Xbox - Hourly Chat Rate per Session model2c step1residuals ~ lin1 + lin2 + postchange 2016-11-17 content change 0.0001140 0.0001140 0.0004629 0.2563250 0.7976999 |

Done

## [1] "Version 0.6 created on 27 Mar 2017 - 20:44:09 took 83.1 secs"